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A filament is mounted on the high voltage end of a Sloan 
radiofrequency resonance transformer and electrons are 
accelerated out of the vacuum tank through a thin window 
at a maximum energy of about 850 kv. About fifty elements 
have been bombarded with these electrons and examined 
with a Geiger counter for induced radioactivity. No posi- 
tive results were obtained. This indicates that the yield is 


Search for Radioactivity Induced by 800-Kilovolt Electrons 


J. J. Livincoop AND ARTHUR H. SNELL,* Radiation Laboratory University of California 
(Received September 17, 1935.) 






less than 1 activation per 10" electrons of approximately 
750 kv, provided that the half-lives have values between 
some seconds and one or two hours and that the products 
of disintegration are electrons of over 200 kv energy. The 
collision cross sections for such reactions have upper limits 
ranging from 10°-** cm? for hydrogen to 10-8 cm? for 
uranium. 








INTRODUCTION 


N the course of development of the Sloan 
radiofrequency high voltage resonance trans- 
former for the purpose of generating x-rays and 
for the acceleration of positive ions, it seemed 
worth while to convert the apparatus into a 
Lenard tube, so that electrons could be directed 
outward into the room where they might be 
conveniently used. High voltage Lenard tubes 
have been studied in some detail by Coolidge! 
but, as far as we are aware, cathode rays have 
never been used in an endeavor to assail the 
nucleus. Theoretically one would expect them 
to have little or no effect, since the de Broglie 
wave-length of an electron traveling with ap- 
proximately the velocity of light is still several 
powers of ten larger than the nucleus. Further- 
more, since the Fermi theory of f-emission 
postulates the simultaneous emission of a neu- 
trino, one would expect that the reverse process 
* 1851 Exhibition Scholar. 


1 Coolidge, Science 62, 441 (1925); J. Frank. Inst. 202, 
693, 722 (1926). Gen. Elec. Rev. 35, 413 (1932). 





2D. H. Sloan, Phys. Rev. 43, 213 (1933); 47, 62 (1935). 
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would require the double entry of an electron 
and a neutrino into the nucleus. If, on the other 
hand, it were found experimentally possible to 
make a nucleus become a #$-emitter by first 
driving only an electron into it, the Fermi 
postulate would seem to require some modifi- 
cation. 


APPARATUS 


The principles of a radiofrequency high voltage 
generator have been described by its designer? 
and its performance as an x-ray tube has been 
published.* It is therefore necessary here only to 
recall that it consists of a powerful short wave 
oscillator feeding into a single resonant circuit 
which consists of a single loop primary and a 
large water-cooled secondary coil both suspended 
in an evacuated metal tank. The secondary coil 
acts as a quarter wave-length antenna; its 
voltage node is at its upper end where it is 
attached to the grounded tank, and its voltage 


Stone, Livingston, Sloan and Chaffee, Radiology 24, 
153 (1935); 24, 298 (1935). 
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Fic. 1. High voltage radiofrequency resonance transformer, 
adapted for the acceleration of electrons. 


loop is at the lower freely hanging end which 
executes high frequency oscillations, of about 44 
meters wave-length, with an amplitude of some 
hundreds of kilovolts. When used as an x-ray 
tube, an anode is mounted on this high voltage 
end of the coil and a filament is placed at 
ground potential in the wall of the tank. Sur- 
rounding the filament is a grid several thousand 
volts negative with respect to it, so as to prevent 
the emission of electrons until the anode has 
neared the peak of its positive voltage oscilla- 
tions, thus permitting an approximation to d.c. 


End of Co:/ : 
Flexible Wire Septum 
Guage Fase / ca J 
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operation. To convert the apparatus into an 
electron accelerator of the Lenard type it was 
necessary to supply a source of electrons at the 
end of the high voltage coil and to mount a 
suitable window in the wall of the tank to allow 
the electrons to escape into the air of the room. 

An early attempt at an electron source was 
simply an eighth-inch tungsten rod sharpened to 
a point and mounted on the end of the coil; 
as the voltage rises to a critical value cold 
emission sets in because of the intense electric 
field, and the source thus supplies its own bias 
to prevent emission at low voltage. This device, 
however, was unstable with regard to quantity 
and quality of emission, the sharpness of the 
point being quite critical, and the vacuum had 
to be lost for every adjustment. 

Consequently advantage was taken of the con- 
struction of the main coil in order to mount a 
hot filament at its end. The coil is made of 7/8- 
inch outside diameter copper tubing within 
which is a copper tube of 7/16-inch diameter, 
the latter being covered with a layer of rubber 
tape to prevent heat transfer from the outgoing 
to the incoming water; the rubber is protected 
against erosion by a layer of thin copper tape. 
The two concentric pipes are therefore insulated 
from each other and serve admirably as lead 
wires to conduct current to the filament from a 
transformer at ground potential. The filament 
assembly is shown in Fig. 2; a copper glass 
seal and a tungsten glass seal are used to bring 
a connection from the inner tube out to the 
vacuum, the 0.080-inch tungsten rod being con- 
nected to the central pipe by a flexible wire for 
ease in assembly. The glass is protected from high 
frequency dielectric hysteresis and from cathode 
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Fic. 2. Details of the filament assembly at the high voltage end of the transformer. The longi- 
tudinal hard rubber strips which align the septum have been omitted in the sketch for greater 
clarity. The direction of the water flow is indicated by arrows. 























RADIOACTIVITY 


rays by the surrounding copper tube and by a 
small shield at the base of the filament. A thin 
cylindrical metal septum, held in place by longi- 
tudinal hard rubber strips (not shown in the 
figure), directs the water down to the end of the 
assembly, thus preventing stagnation and conse- 
quent boiling. By moving the rounded solid 
copper tip forward or backward the space charge 
around the filament can be considerably altered 
so that emission during the low parts of the 
voltage waves can be reduced. The filament is 
20 cm from the tank wall. 

This arrangement, thougn satisfactory, is not 
as good nor as flexible as would be an inde- 
pendently biased grid whose potential is con- 
trolled by a high voltage wire passing down the 
inside of the coil, but there are difficulties in 
obtaining a sufficiently well insulated wire that 
is small enough not seriously to hinder the flow 
of water. An alternative possibility would be to 
build into the filament unit a midget trans- 
former, rectifier tube and filter condenser to 
supply the potential for a grid, driven through 
a circuit consisting of the main coil and a low 
voltage wire inside it. 

The window for the electrons consists in. a 
copper foil 0.001 inch thick soldered to a copper 
plate 2 inches in diameter and 3/16 inch thick 
and perforated by many 1/8-inch holes spaced 
5/32 inch between centers. A water-cooling pipe 
is silver soldered to the periphery of this disk, 
which in turn is soldered to a larger brass plate 
attached with sealing wax to an aperture in the 
tank wall. The electron beam strikes the window 
over an area about 1 inch in diameter and 
currents of 300 microamperes (average, read on 
a d.c. meter) have been passed through without 
damage. Windows with larger holes covered with 
copper, nickel or molybdenum foils melted 
through in a few seconds, even though addi- 
tionally cooled with jets of air or water. Scatter- 
ing of the electrons in the foil and air is sufficient 
to bend the beam backwards so that a hemi- 
spherical glow can be seen, in the darkened 
room, of about six feet radius. All sealing-wax 
joints and rubber tubes in the immediate vicinity 
of the window must be shielded with metal, else 
a very rapid deterioration sets in due to the 
bombardment. 

Knowledge of the energy of the electrons has 
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been gained by rather rough absorption measure- 
ments. A plane collecting plate was mounted on 
an insulated stand about 2 inches from the 
window and connected to ground through a 
microammeter ; absorbing foils of aluminum were 
placed between it and the window. To prevent 
the registration of electrons scattered around the 
foils, the collecting electrode was shielded with a 
diaphragm consisting of a large sheet of 1/4-inch 
lead with a 1-inch hole cut in it. Correlation 
between range and energy was made through 
the data of Varder,‘ the maximum range of 
0.048-inch aluminum indicating an energy of 
845 kv. To obtain this potential, the average 
power dissipated as heat in the coil, determined 
from the temperature rise of the cooling water, 
was 22.5 kw. The anode power supply (60-cycle 
a.c. at 9000 volts) totalled 70 kw, including the 
energy lost in the control resistances. The heating 
current for the oscillator filaments required 15 
kw additional. Absorption curves were taken at 
a number of potentials applied to the anodes of 
the oscillator tubes and the heat dissipated in 
the coil was noted in each case. Comparison 
between the voltage developed (electron range 
plus window correction) and the power in the 
coil shows that these two quantities are quite 
well related by the equation 


Kilovoltsccoit) = 183 (Kilowattscecoit )?. 


This form of equation was to be expected from 
elementary considerations, although the experi- 
mental value of the constant, of course, depends 
upon the inductance and capacity of the par- 
ticular coil and vacuum tank used, as well as 
on the frequency and wave shape. 


SEARCH FOR NUCLEAR DISINTEGRATION 


Because of the large intensity of the incidental 
x-rays, it was feasible at this time to search only 
for delayed radioactivity of targets bombarded 
by the electron beam after the exposure had 
been made and the high voltage shut off. The 
points of unprotected Geiger counters (con- 
nected to an amplifier and scale of two thyratron 
recorder) were quickly rendered unstable and 
useless by scattered electrons, even though no 
potential was applied to the point during bom- 


* Varder, Phil. Mag. 29, 726 (1915). 
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TABLE I, Elements bombarded and examined for radioactivity. 











Z Form Usep | Z 








Form Usep 
1 H' H;PO, 134 Se Metal 
1 D? Ds3PO, 135 Br KBr 
3 Li Metal 38 Sr SrCl. 
4 Be Metal 140 Zr Metal 
5 B Pyrex, borax '41 Cb Metal 
6 C © solid 142 Mo Metal 
7 N  Ca(NOQs)2 '47 Ag = Metal 
8 O Various compounds! 48 Cd Metal 
9 F Ba Fe. /50 Sn Metal 
11 Na NaCl 51 Sb Metal 
12 Mg Metal 52 Te Metal 
13 Al Metal 53 I KI 
14 Si Glass, sand 56 Ba BaF. 
15 P HsPO, 58 Ce  Ce(NOs)s 
16 S Solid 60 Nd Nd (NQs)3 
17 Cl NaCl, KCl 73 Ta Metal 
19 K- KCI, KI 74 sW Metal 
20 Ca Cal ), Ca(N¢ )3)2 78 Pt Metal 
24 Cr KoCr,O; 79. Au Metal 
25 Mn MnCl, 80 Hg Hg—Cu amalgam 
26 Fe Metal 81 TI Metal 
27 Co CoCir 82 Pb Metal 
28 Ni Metal 83 = Bi Metal 
29 Cu Metal 90 Th Th(NOs), 
30 Zn Metal | 92 U U(NOs3)3 
33 As Aso ds 


| 








bardment. Consequently, lead shielded counters 
were used and a system of levers and control 
strings arranged such that the operator could 
move the target into the electron beam and 
close the opening to the counter or, after the 
exposure, open the counter and swing the target 
into position in front of it. This shift could be 
made within two seconds, from the operator's 
station in a lead shielded booth across the room. 
Some of the targets were in the form of metallic 
disks, others as powder held to a base plate with 
a thin layer of sealing wax. All of the latter 
type and many of the former had to be water 
cooled. The targets were bombarded for three 
minutes and examined with the counter for five 
minutes, this procedure being repeated several 
times for each. The electrons of all velocities 
totalled about 100 microamperes and the usual 
peak energy was somewhat over 800 kv. 

Table I lists the elements bombarded and 
examined. In no case was there any evidence of 
induced radioactivity. 

At a very conservative estimate, based on the 
absorption curves, at least 1 microampere of the 
bombarding electrons had an energy of 800 kv 
or over. Counts numbering three times the back- 
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ground would certainly have been recognizable 
so that, taking into account the solid angle sub- 
tended by the counter, a yield of 1 disintegration 
per second would have been significant, provided 
that the half-lives were not shorter than about a 
second nor longer than perhaps an hour, for 
the bombardments of each element totalled 
about 10 minutes, and 2 seconds were lost before 
counting began. If it be assumed that the 
products of disintegration are electrons of not 
less than 200 kv energy, they could emerge from 
the target and reach the counter after starting 
from a depth as great as that corresponding to 
0.043 gram/cm?, in which thickness the 800 kv 
bombarding electrons are reduced in energy to 
710 kv. We may conclude, therefore, that if 
electrons of approximately 750 kv _ produce 
radioactive nuclei of half-lives measured by a 
few minutes and emitting 8-particles of 200 kv 
energy or over, the efficiency is not more than 
1 disintegration for 10" electrons and that the 
collision cross sections range in value from not 
greater than 10-* cm? for hydrogen to 10-* cm? 
for uranium. 

Pokrowski> has stated that x-rays of 90 to 
140 kv produce radiations from lead ; the activity 
reported was irregular and although it decreased 
rapidly in the first few hours following exposure, 
it apparently accumulated in intensity over 
many days of intermittent activation. Therefore 
in the present work not only was lead examined 
in the routine manner of the other elements, but 
a piece of lead from the shielding around the 
vacuum tank was investigated, after it had been 
exposed to the x-rays produced during several 
weeks operation. No counts above background 
were observed even when the lead was placed so 
close to the counter that particles of 2-mm 
range would certainly have been detected. 

It is a pleasure to acknowledge the friendship 
and enthusiastic support of Professor E. O. 
Lawrence. One of us (A.H.S.) is grateful to the 
Royal Commissioners of the Exhibition of 1851. 
The experiments have been aided by grants from 
the Research Corporation and the Chemical 
Foundation. 


5 Pokrowski, Phys. Rev. 38, 925 (1931). 
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Two Radioactive Substances from Magnesium after Deuteron Bombardment* 


M. C. HENDERSON, Radiation Laboratory, Department of Physics, University of California 
(Received September 16, 1935) 


Under bombardment by 3.3 MEV deuterons, magnesium 
becomes radioactive with the emission of negative electrons 
and gamma-rays, and with two decay periods: 10}+} min. 
and 15.8+0.5 hours. These periods are characteristic of 
Mg”? (Fermi) and Na* (Lawrence). The reactions produc- 
ing these bodies are: 

122Mg**+,H? =;2Mg??+,H!; 12Mg?? =,;Al"7+8+y7. 

i2Mg**+,H? =,,Na*+ 2Het; uNa* = ;2Mg*+8+7. 
The protons have not been looked for. The alpha-particles 
have a range of about 7 cm. The Mg?’ beta-rays have 2.0 


MEV maximum energy, and the accompanying hard 
gamma-ray has about 1.3 MEV energy. The voltage 
excitation functions of these two activities have been 
measured, and have been found to follow different laws. 
The Mg*’ following the Oppenheimer and Phillips function 
and the Na* the older Gamow function. At 3 MEV one 
atom becomes Na*™ for each 9.6 atoms going to Mg?’. 
The thick target yield of Mg?’ atoms is about 5 per 10’ 
deuterons. The nuclear cross section (Mg**) at 3 MEV is 
about 2.4 10~**, 





1. INTRODUCTION 


MMEDIATELY after the discovery by I. 
Curie and F. Joliot' that artificial radio- 
activity could be produced by alpha-particle 
bombardment, several laboratories followed up 
their suggestion that it might be produced by 
deuteron and proton bombardment as well. The 
first communication from this laboratory on this 
subject reported? the results of preliminary work 
with elements in the first two rows of the periodic 
table. It was found at that time that several of 
them showed strong radioactive properties after 
bombardment with 3-MEV deuterons. The 
ionizing radiations emitted were ‘‘particles of 
electronic mass,’’ but were not identified as to 
sign. For magnesium the half-life obtained was 
about nine minutes. 

Since these first observations the apparatus 
has been modified* so that it is now possible to 
bombard the target with the desired ion, either 
in vacuum or in air, and to remove it immediately 
for measurement without disturbing the vacuum 
in the ion production chamber. Last spring work 
with magnesium was resumed, with the object of 
measuring the half-life more accurately and in- 
vestigating the activity more in detail. 


* Reported at the Los Angeles Meeting, June, 1935. 
Abstract 5. 

11. Curie and F. Joliot, Nature 133, 201 (1934). 

2M. C. Henderson, M.S. Livingston and E. O. Lawrence 
Phys. Rev. 45, 428 (1934). 

3 E. O. Lawrence and M. S. Livingston, Phys. Rev. 45, 
608 (1934). See also diagram in Phys. Rev. 47, 453 (1935). 


2. APPARATUS 


The source of deuterons was the large cyclotron 
in the Radiation Laboratory. At the time the 
magnesium work was carried on, it was adjusted 
to furnish deuterons of energy of about 3.3 
MEV. The energy of the ions produced by a 
cyclotron can be approximately calculated from 
the magnetic field or from the wave-length 
employed, but for experimental purposes it is 
always measured by determining their range in 
air, after they have passed out of the apparatus 
through a thin aluminum window. During most 
of the runs this range was 11.5 cm, taken as 
corresponding to deuterons of about 3.36 MEV. 

The target was in all cases a piece of mag- 
nesium ribbon of 99.97 percent purity. Before 
bombardment the surface was cleaned with 
emery paper and washed with petroleum ether 
to remove oil film. The first sample was bom- 
barded in air, but all subsequent ones were 
bombarded in vacuum to avoid oxidation and 
contamination of the surface by radioactive 
oxygen.# 

‘ All the measurements of activity were carried 
out on the Lauritsen-type quartz fiber electro- 
scope of simple design that has been used for 
nearly all the electroscope measurements in this 
laboratory. A thin aluminum window, against 
which the radioactive source was placed, ad- 
mitted the beta and gamma-rays to the ioniza- 
tion chamber, and the absorption measurements 


4E. M. McMillan and M. S. Livingston, Phys. Rev. 47, 
452 (1935). 
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were carried out by placing aluminum or lead 
sheets between the source and window. Un- 
fortunately the intensity of the radiations was 
not great enough to allow the source to be more 
than a few centimeters from the window, and 
consequently the geometry of the arrangement 
left something to be desired. 


3. THe Decay PERIODS 


A thick target of magnesium was bombarded 
in air for about 10 minutes by an ion beam of 
3 MEV energy and about 0.9 microampere 
intensity. At the end of the exposure the target 
was at once placed at the window of the electro- 
scope. The measured activity, as a function of 
time, is plotted on a logarithmic scale in Fig. 1, 
upper curve. The half-life determined from the 
linear portion, after due allowance for the 
natural leak and background activity of longer 
life, is 10.25+0.25 min. The points taken in the 
first eight minutes show a more rapid decay as a 
result of the presence of radio-oxygen, which 
was deposited from the air around the target 
during the bombardment. 

After about an hour and a half the successive 
readings showed definitely the presence of some 
substance with a longer life: in fact the activity 
two hours after exposure was still 20 times the 
natural leak. The target was therefore replaced 
in the cyclotron and the bombardment resumed 
with a view to intensifying the long-lived 
activity. After waiting two hours for the ten- 
minute substance to decay completely, a reading 
of the residual activity was taken and the decay 
was then followed for the next 24 hours. The half- 
life proved to be 15.8+0.5 hours. These two 
periods were checked several times in the course 
of later runs that were made for other purposes. 
The errors are conservatively estimated. 

To determine the sign of the ionizing particles 
they were observed visually in a Wilson cloud 
chamber that was placed in a magnetic field of a 
few hundred gauss.‘ All the particles, both those 
from the ten-minute and those from the sixteen- 
hour activity, proved to be negative electrons. 
No positrons could be detected. This fact rules 
out a carbon contamination as a possible source 


5] am indebted to Dr. Kurie for the use of the cloud 
chamber. 


HENDERSON 
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Fic. 1. Upper curve: decay of the activity of a magnesium 
target after deuteron bombardment. Lower curve: ab- 
sorption of the beta-rays from Mg” in Al. 


of the ten-minute activity, since after deuteron 
bombardment carbon becomes radioactive with 
the emission of positrons. 


4. THE AcTIVE SUBSTANCES 


It is thus apparent that magnesium combines 
with deuterons to produce at least two radio- 
active substances. Their identity is at once 
apparent from the work of others, since the 
periods are characteristic of two substances 
already discovered: radio-sodium, Na”, de- 
scribed by Lawrence,* having a half-life of 
15.5 hours and produced by the action of 
deuterons on sodium; and radio-magnesium, 
Mg?’, discovered by Fermi,’ half-life ten minutes 
and produced—among other ways—by the action 
of slow neutrons on magnesium. 

Known types of reaction, when applied to the 
Mg’* isotope, lead at once to Mg?’ and Na*. 
Accordingly, we have: 


12M g**+ |H? = 12Mg”’+,H!; 
wMg?? = ;Al"+8+7; 

12Mg?*+ ,H?= ,,Na**+ He; 
11Na*4= »»Mg*+8+y. 


6 E. O. Lawrence, Phys. Rev. 47, 17 (1935). 
7 Fermi, Amaldi and others, Proc. Roy. Soc. A146, 483 
(1934); Proc. Roy. Soc. A149, 522 (1935). 
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RADIOACTIVE SUBSTANCES FROM MAGNESIUM 857 


Many cases are now known in which the re- 
action between two nuclei leads to several pairs 
of end products. The two reactions just given 
are, however, the first in which two radioactive 
substances have been produced from a single 
substance by the use of a charged particle as a 
projectile. It seems reasonable to believe that 
this type of branching should be quite general 
and that not only Mg*® but also O'8, Ne”, Si*® 
and other atoms may also yield more than one 
radioactivity under deuteron bombardment. 

The protons emitted in the first reaction have 
not been looked for as yet. The alpha-particles, 
on the other hand, have already been observed, 
although their number was too small for accurate 
measurement.® 


5. THE RADIOACTIVE RADIATIONS 


The beta and gamma-rays from Na™ have 
been studied in detail by Lawrence (reference 6 
and also unpublished work in this laboratory) 
and there seemed to be no point in repeating his 
experiments with an intrinsically weaker source. 
On the other hand, the beta-rays from Mg?’ 
were measured by Fermi using a source obtained 
by neutron bombardment, and therefore neces- 
sarily much weaker than one obtained by the 
direct deuteron-magnesium reaction in the 
cyclotron. 

An absorption curve of the radiation from a 
Mg?’ source is shown in Fig. 1, lower curve. All 
points are corrected for decay and for natural 
leak. It will be seen that there is a component 
of hard radiation present along with the easily 
absorbed beta-rays. When this component is 
subtracted from the measured total intensity the 
steep linear curve results. Extrapolating this 
linear curve to the level of the natural leak gives 
for an “end point’’—the maximum beta-ray 
range—the point marked with a cross. This end 
point falls at 0.95 g/cm? of Al. Using Feather’s 
empirical rule—E =(R+0.091)/0.511, where E 
is in MEV and R in g/cm?—the maximum beta- 
ray energy is 2.05 MEV. Fermi found a half- 
value absorption thickness of 0.07 g/cm? in Al, 
in fair agreement with the data in Fig. 2. 

The point corresponding to 2.05 MEV and 


8G. N. Lewis, M. S. Livingston and E. O. Lawrence, 
Phys. Rev. 44, 55 (1933). 
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DEUTERON RANGE IN CM 


Fic. 2. The voltage transmutation functions of Mg?’ 
and Na*. 


10.25 minutes, when placed on Sargent’s curves 
of beta-ray periods and energies, lies midway 
between two branches. If the maximum energy 
were 3.8 or 1.2 MEV it would lie on one or the 
other of them, but the latter values seem to lie 
well outside the limits of experimental error. 
The absorption coefficient of the gamma-rays 
was determined in a separate experiment. A 
value of u =0.74 (cm~") of lead was obtained. The 
geometrical conditions were poor, but this value 
probably corresponds to an energy of about 
1.3 MEV. From the intensity of the ionization 
produced by the gamma-rays, nearly 2 percent 
of that produced by the beta-rays, it appears 
that one gamma-ray is emitted for each beta-ray. 


6. THE VOLTAGE TRANSMUTATION FUNCTION 


The variation in the yield of transmutations as 
a function of the energy of the incident particle 
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has been called the excitation function, or more 
properly the voltage transmutation function. 
The first theoretical treatment of this variation 
and a method.for calculating it in any given case 
was first given by Gamow in Atomic Nuclei and 
Radioactivity. The formulas are obtained by 
considering, among other things, the variation 
with energy in the probability that the projectile 
as a whole penetrate the nucleus. For most of 
the reactions thus far studied, in which protons 
and alpha-particles are the projectiles, the 
formulas thus obtained fit the data well. How- 
ever, it has been found in this laboratory® that 
the yield of transmutations caused by deuterons 
does not fit the Gamow curve at the higher 
values of deuteron energy, at least in the re- 
actions in which a proton is emitted. 

To meet this difficulty, Professor Oppenheimer 
and M. Phillips’? have recently developed a 
theory for reactions of this type, one of the under- 
lying assumptions of the theory being that the 
deuteron as a whole does not penetrate the target 
nucleus. Only the neutron in the deuteron enters 
the nucleus, the proton being “‘rejected.”’ Accord- 
ing to their calculations the transmutations 
should begin at a lower voltage than on the 
Gamow theory and should not increase as rapidly 
in number with increasing deuteron energy." 

Transmutations in which a radioactive body 
is produced are particularly convenient for 
studies of the transmutation function. Conse- 
quently, as soon as it was established that two 
radioactive bodies were obtained from mag- 
nesium, the attempt was made to measure the 
transmutation functions of both reactions. The 
measurement is particularly important because 
one of the bodies, Na*, is produced by the 
penetration of the deuteron as a whole and the 
subsequent expulsion of an alpha-particle, while 
the other, Mg?’, is produced by the ‘‘rejection’”’ 
of a proton and the absorption of a neutron. 
Thus the two activities should be different in 
their dependence upon voltage, the one following 
the older Gamow formula and the other the 
Oppenheimer theory. 


®E.O. Lawrence, E. M. McMillan and R. L. Thornton— 
TON in Science 81, 421 (1935); Phys. Rev. 48, 493 
(1 ' 

” Oppenheimer and Phillips, Phys. Rev. 48, 500 (1935). 

" T am indebted to Professor Oppenheimer for calculat- 
ing the yield curve according to his theory. 
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Several runs were made to determine the 
transmutation functions by using the now stand- 
ard method of exposing a pile of thin magnesium 
foils to the action of the deuteron beam and 
then separating the foils and measuring their 
activity separately. This method has several 
advantages: It gives the ‘differential’ trans- 
mutation function directly, which is easier to 
compare with theory than the integral function ; 
and it also ensures that all foils receive exactly 
the same exposure, since the deuteron beam 
passes through them all and since the number 
of deuterons that react in any given foil is a 
minute fraction of the total number passing 
through. The top foils in a pile thus act as 
absorbers to reduce the energy of the deuterons 
that reach the lower foils. 

Magnesium is sufficiently malleable to be 
reduced with ease to foils a few ten-thousandths 
of an inch thick. The thickness was measured 
directly with a micrometer, and also more 
accurately, and for these experiments, more 
significantly by measuring their stopping power 
for polonium alpha-particles. For magnesium 
0.001 inch is the equivalent of 2.5 cm of normal 
air (760 mm and 20°C). The thickness found 
convenient was about 1.4 cm of air, 0.00056 
inch, or 0.0014 cm. 

Before making an exposure each foil was 
carefully measured for thickness and was then 
stacked in a pile with a ‘“‘mask”’ of thin copper 
sheet. This mask was placed on the pile to pro- 
tect the edges of the foils from the deuterons and 
to prevent irregularities or cracks in the edges 
from allowing deuterons to leak through to lower 
members of the pile. After exposure in vacuum 
to about 0.8 microampere for 13 hours, the pile 
was removed, separated and the activity of 
each foil measured; the weakest first. After a 
lapse of two hours, sufficient to permit the Mg?? 
to decay completely, the foils were again meas- 
ured for the Na* activity. Since the Na™ activity 
was never more than about 1 percent of the Mg?’ 
activity at the time the latter was measured, no 
correction has been applied for it. The activity 
of each foil, of each kind, was then corrected for 
decay since the end of the exposure and also 
corrected for any difference in thickness there 
might be between foils, since a thicker foil is 
proportionately more active. 
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The results of the last three runs are shown 
in Fig. 2. To bring out the difference in transmu- 
tation function, all the activities have been made 
comparable by making the points at highest 
energy coincide. There are thus six experimental 
points at 11.5 cm and 18 div/sec. The ordinates 
for Na* have been multiplied by approximately 
200 to make this fit. The two solid curves are 
drawn in accordance with the two theories of the 
transmutation function, and have been adjusted 
to pass approximately through the mean of the 
points around 9 cm. The reason for using these 
points instead of the 11.5-cm point is that the 
top foil in a pile loses active atoms by recoil to 
the foil next beneath it, while the lower foils 
gain by this mechanism approximately as much 
as they lose. 

It is readily seen that the transmutation func- 
tion of the reaction producing Mg?’ is quite 
different from the corresponding Na*™ function, 
and the goodness of fit furnishes strong support 
to both theories. 

The Gamow formula for the variation in 
number of transmutation, N, with voltage, V, 
is the following: N=kV—e-S, where k is a 
constant of proportionality and S is defined by: 


4re* Mp i 
s-()[au( 2) 
h 2Ep 


+ ZaLa( 





M,: Mya )| 
2(Mat+Mya)(Q+£Ep*) 

Z and M are the atomic numbers and masses, 
in grams, of the atoms indicated by the sub- 
scripts; Q, Ep and Ep* are the energies, in ergs, 
of the transmutation, of the deuteron, and of 
the deuteron with respect to the center of mass 
of the system: i.e., 


Ep*=Ep: Mug (Mp+Mm;,). 


The first term inside the bracket takes account 
of the penetration of the deuteron through the 
nuclear barrier, the second term the penetration 
of the alpha-particle outwards through the 
sodium nuclear barrier. In actual numbers, 


S=17.0/V'+42.5/((14/13)0+ V)!, Vin MEV. 


In the calculation it is necessary to assume 
some value of Q, the energy of transmutation. 


An accurate knowledge of the range of the alpha- 
particles produced would of course lead at once 
to the value of 0. However, at the time the alpha- 
particles were found® the energy of the deuterons 
was only 1.2 MEV and the currents of the order 
of 10-* ampere. A calculation from the Gamow 
formula and the yield found at higher energies 
shows at once that the number of alpha-particles 
to be expected under the old conditions was 
just on the limit of measurement, 2 or 3 a minute. 
Consequently no great accuracy is to be ascribed 
to the value found: 6 cm. A reconsideration of 
the data leads to 7 cm as the most probable 
value, corresponding to a Q of 8 MEV. More 
data are needed on this point. 

In Fig. 2 two curves have been drawn for 
different Q’s: the solid curve for Q=10 MEV, 
the dotted for 5 MEV. From the fit Q=10+2 
MEV. Imperfections in the magnesium foils will 
lead to too high a value of Q. 

One may also estimate Q from the masses and 
energies concerned. By writing the complete 
equations of the Na* reaction and adding them: 


12 Mg?*+ )H? = 1.Na**+.He*+(Q. 
1.Na*4= i2Mg*4*+8+7. 


(12Mg?*— 12Mg™) +,H?—»He'=Q+8+y. 


The difference in mass between the two mag- 
nesium isotopes has not been accurately meas- 
ured. In the parallel case of neon 20 and 22 
the difference in mass is 1.998+0.001."% One 
would offhand expect the Mg*® isotope to be 
about 2.000 heavier than Mg*, since there are 
two more neutrons in the nucleus than are needed 
to form complete alpha-particles; but on the 
other hand, both the number of protons and 
the number of neutrons in both nuclei are even, 
so the general trend of the mass defect curve 
would lead to about 1.998 for the difference. 
Taking 2.000+0.002 for this difference and 
Oliphant’s values for H? and He‘, one gets 2.000 
+2.014—4.003 =Q0+6+y- =0.011+0.002. From 
Lawrence’s work® and from unpublished work in 
this laboratory, 8+y is about 4 MEV, and hence 
Q is about 7+2 MEV, in rough agreement with 
the value given by the Gamow curve. 

In contrast with the Na*™ reaction, no sensible 


2 Bainbridge, Phys. Rev. 43, 424 (1933). 
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adjustment of the Gamow formula gives a fit with 
the Mg?’ data. Even taking Q as infinite, tanta- 
mount to neglecting the second term in S, still 
gives a curve rising too steeply. The Gamow 
formula after making appropriate changes in 
the Z’s and M’s, becomes S=17.0/V!+12.3/ 
((14/13)Q+ V)!. The solid curve drawn through 
the Mg?? points has been calculated from the 
theory of Oppenheimer and Phillips! according 
to the equation (Eq. (12)). N=kV—e-*’, where 
S’ is defined as 


4Ze / My} 
S= (—) F(W/1). 
h I 








In this equation, F(W/J) is the minimum value 
of F(W/I), the values of which are taken from 
their plot of this function. J is the binding 
energy of the deuteron, 2. its mass, and W its 
kinetic energy. The only adjustable constant is 
the binding energy of the deuteron. This has 
been taken to be 2.0 MEV. 

Although Q for the reaction in which Mg?’ is 
produced does not enter the formula for the 
transmutation function, it is of interest to esti- 
mate it. As before, writing the equations of 
production and decay, one has: 


12Mg?*+ ,H?= 12Mg??+,H'+Q, 
12Mg?? = ,3Al"+8+y. 
and adding them: 
(12M g?®— ,3Al?’) + ,H?—,H'=Q+8+y7, 


Al*?? should be at least 1.000 heavier than Mg*® 
from the mass defect curve. Thus one has 
— 1.000+ 2.014 — 1.008 = 0.006. 8+7 is about 3.3 
MEV. Q is therefore about 3 MEV. This con- 
sideration supplies an upper limit to the mass of 
Al*’, since if Q@=0, Al??? — Mg?*= 1.003. 


7. THE BRANCHING RATIOS AND THE ABSOLUTE 
YIELDS 


When a target is exposed to a constant beam 
of deuterons, any radioactive atoms that may be 
formed increase in number with the expression 
(1—e-7/'av), Consequently, after an exposure of 
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two hours the number of Mg?’ atoms in a mag- 
nesium target has practically reached a constant 
maximum value, while the number of Na™ 
atoms is only 8.5 percent of the maximum 
number. This maximum number represents, for 
each activity, a state of affairs in which as many 
atoms are being formed per second as are de- 
caying. Since the number decaying per second is 
what is measured as “‘activity,’’ the ratio of the 
final maximum activities is the ratio of the prob- 
ability of formation of the two kinds of atom. 

In one experiment the activities in the top 
foil, 1.4-cm air equivalent, at the end of a two 
hour exposure to 0.8 microampere were 19.5 
div/sec. due to Mg?’ and 0.18 due to Na*™. 
The saturation values of these activities, allowing 
for the finite time of exposure, are therefore 19.5 
and 2.1. Hence for each atom of Na*™, 19.5/2.1 
=9.3 atoms of Mg?’ are produced. While this 
ratio is different at lower voltages, its variation 
with voltage may be readily calculated from 
Fig. 3.% 

The absolute number of radioactive atoms 
produced is much more difficult to estimate 
accurately, depending as it does on the constancy 
of the deuteron beam, the transmission of the 
window, the calibration of the electroscope, and 
other factors. However, assuming a transmission 
of }—certainly too high—and the electroscope 
calibration of 2X10‘ beta-particles/sec. for 1 
div/sec. ; one finds that there are 1.6 Mg?’ atoms 
produced for each 107 deuterons. The integrated 
yield, for a thick target, is about 3 times this, 
or 5 per 10’. 

Transmutation probabilities are often ex- 


13 In this and the other similar calculations no allowance 
has been made for the different energies of the Mg?’ and 
Na*‘ beta-rays. 
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pressed in terms of the cross section of the 
nucleus, assuming that every deuteron which 
strikes within that cross section produces a 
transmutation and that all others do not. For 
Mg”**, 11 percent of the total number of atoms of 
Mg present, this cross section, from the data 
given, is 2.4X10-*° cm? at 3 MEV. Both this 
cross section and the yield as given above are 
certainly too low. 


I am very greatly indebted to Professor 
Lawrence for the opportunity of using the 
cyclotron and for the facilities of the laboratory. 
To Professor Oppenheimer and Dr. McMillan I 
owe many thanks for helpful discussions and 
cooperation with some of the calculations. 
Thanks for financial assistance to the laboratory 
are also due the Research Corporation and the 
Chemical Foundation. 
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The infrared spectrum of GeH, has been measured to 
determine the fundamental vibration frequencies and to 
investigate the rotational structure of the bands. Intense 
absorption regions were found, enumerated in the order of 
their intensities, at 4.74u (2110 cm~), 10.74 (934 cm™), 
3.4u (2922.7 cm™ and 3031 cm™) and 2.34 (4300 cm) 
which have respectively been identified as v3, v4, vi+v4, 


HE recent measurements on the infrared 

spectrum of methane! and silane? have 
shown that the observed bands conform very 
well in a general way with the predictions of the 
theory for molecules of the regular tetrahedral 
type, but that there exists a definite complexity 
of the rotational structure in many of the bands 
which is not readily accounted for. This lack of 
agreement between observation and theory and 
the obvious desirability of finding the values of 
the fundamental vibrations has seemed sufficient 
to warrant investigating the infrared spectrum 
also of germane (GeH,) a molecule which also 
must certainly belong to the tetrahedral XY, 
type. 

A sample of chemically pure GeH, was ob- 
tained as a gift from Professor Warren C. 
Johnston of the Department of Chemistry, Uni- 
versity of Chicago. Like silane (SiH,) this gas 
should be manipulated in the absence of air or 


1A. H. Nielsen and H. H. Nielsen, Phys. Rev. 48, 864 


(1935). 
2 W. B. Steward and H. H. Nielsen, Phys. Rev. 47, 828 
(1935). 


v3+v4, and 2v;. A much weaker absorption peak near 5.74 
is identified as »2+»4. From the spacings between the 
principal rotation lines in the fundamentals »; and » a 
value 7.0X10~-* g-cm? for the moment of inertia of the 
molecule and values 1.37X10-§ and 2.06X10°* cm, 
respectively, for the Ge-H and the H-H distances are 
indicated. 


oxygen to guard against explosions. Two absorp- 
tion cells of glass fitted with windows of polished 
rocksalt were used, one 30 cm long filled with gas 
to a pressure of 50 cm of Hg for measurements on 
weaker bands and another 10 cm long filled with 
gas to 10 cm Hg pressure for measurements on 
the more intense regions of absorption. 

A set of preliminary measurements extending 
from 1.0u to 13.0u was made with a prism 
spectrometer of the Wadsworth type, to de- 
termine the approximate positions of the prin- 
cipal regions of absorption and their relative 
intensities. Four fairly intense regions were 
found, enumerated in the order of their intensities 
at 4.74y, 10.7y, 3.44, and 2.34. A much weaker 
peak was found near 5.74. Of these the three 
regions at 3.4u, 4.744 and 10.74 have been 
studied under high dispersion using a prism 
grating spectrometer described in relation to 
earlier experiments, and two different echellette 
gratings ruled by Wood, one with 4800 lines to 
the inch for the 4.0u region and another with 800 
lines to the inch for the 10y region. 
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In these regions, spectrometer slits equivalent case. Moreover, all three of these regions are 
to frequency intervals in the spectrum of 0.8 quite free from lines of atmospheric absorption 
cm“, 0.8 cm™ and 0.6 cm™, respectively, could so that the data have been recorded throughout 
be used and galvanometer deflections were taken in terms of galvanometer deflections. The absorp- 
at divisions on the circle equal to the spectral tion patterns for these regions obtained by 
range subtended by the spectrometer slit in each plotting galvanometer deflections as ordinates 
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TABLE I. 3.3u region of Gell,. 


TABLE II. 4.74 region of Gell,. 











LINE FREQ. LINE FREQ. | LINE FREQ. 
No. (cm~) No. (cm™!) No. (cm™'!) 
2841.6 8 2956.7} 17 3065.8 

1 2850.3 9 2962.8 3071.8 
2 2864.8| 10 2970.1 3078.1 
2868.0} 11 2977.2| 18 3084.3 

3 2875.7 2983.7 3090.0 
2885.0} 12 2986.0 3095.8 

+ 2890.0 2990.2} 19 3102.8 
2894.9} 13 2996.1} 20 3109.2 

5 2901.1 3002.0 3113.7 
2905.0 3006.8 3117.6 

2908.3 | 14 3013.4} 21 3121.5 
2915.1 3023.7 3126.0 
6 2922.7| 15 3031.0} 22 3131.9 
2936.0} 16 3054.9 3138.3 
7 2945.5 3059.9 3141.0 











against a scale of frequencies and wave-lengths 
are given in Figs. 1, 2 and 3 and accompanying 
each figure is a tabulation of frequencies of th: 
principal lines (Tables I-III). 

The region near 3.34 apparently corresponds 
to the 2.34 region in the methane spectrum and 
the 3.2u region in the silane spectrum, and like 
these it consists of two badly overlapping bands 
with centers at 2922.7 cm and 3031.0 cm. 
By analogy to the spectra of methane and silane 
these two bands have been identified as »;+ 9, 
and v3+v,4, respectively. Because of the over- 
lapping, it is impossible to find any regularity 
of line spacing except perhaps on the high fre- 
quency side of the band at 3031 cm™ where an 
average over several nearly equally spaced lines 
leads to a value of 6.2 cm for the separation 
between lines. 

The 4.74 band with its center at 2110 cm™ 
undoubtedly corresponds to the 3.34 band in the 
methane spectrum and the 4.5 band in the silane 
spectrum and is accordingly designated as 13. 
It is the most intense of the bands measured in 
the germane spectrum. Like the corresponding 
bands in the spectra of methane and silane, it 
consists of a strong central line with a series of 
sharp rotation lines on either side forming P and 
R branches. As may be seen from Fig. 3, there 
occurs also here a background which we do not 
believe can be ascribed to observational errors 
since it has been repeated almost exactly a 
number of times. The background is here not as 
predominant as in the corresponding band in 


‘the silane spectrum and neither is the Q branch 


so extremely broad, so that as a whole this 


LINE Fre. | LINE FreQ. | LINE FREQ. 
No (cm~!) | No. (cm™'!) No. (cm!) 
1 1998.9; 15 2082.7 | 29 2166.1 
2 2005.1 16 2088.4 | 30 2171.3 
3 2011.7} 17 2094.3} 31 2176.5 
4 2017.5} 18 2099.8 | 32 2181.8 
5 2023.3} 19 2110.1} 33 2187.1 
6 2029.4 | 20 2116.8) 34 2192.0 
7 2035.3} 21 2122.4} 35 2197.1 
8 2041.4} 22 2128.0| 36 2202.4 
9 2047.5 | 23 2133.6| 37 2207.1 
10 2053.4 | 24 2138.9} 38 2212.2 
11 2059.3 | 25 2144.6| 39 2217.0 
12 2065.1 | 26 2150.0} 40 2222.1 
13 2071.2} 27 2155.1) 41 2226.9 
14 2077.1} 28 2160.9 














TABLE III. 10.74 region of Gell,. 




















LINE FREQ. | LINE FREQ. | LINE FREQ. 
No (cm!) | No. (cm=!) | No. (cm™!) 
l 870.0 913.5| 16 981.1 
871.3 917.1 983.5 

2 873.3| 10 921.4) 17 986.9 
875.2 926.6 990.2 

3 877.9| 11 934.0; 18 992.8 
4 881.2} 12 944.5 998.7 
883.8 949.9; 19 1004.7 

5 885.8 951.5 1010.8 
6 889.2} 13 954.6; 20 1017.0 
7 898.4 958.5 1018.8 
901.1 14 963.2 1022.4 

8 904.3 967.2) 21 1028.7 
906.8} 15 971.7 1034.2 

9 912.1 974.9| 22 1037.9 








region resembles the corresponding one in the 
methane spectrum a great deal more than that 
in the silane spectrum. Some convergence of the 
principal lines is observed towards higher fre- 
quencies and the average separation between 
lines taken over the entire band is found to be 
5.56 cm™. 

The long wave-length region has its center at 
934 cm and is taken to correspond to the 7.74 
band in the methane spectrum and to the 
9.54 band in the silane spectrum. It is conse- 
quently identified as v4. It should be pointed out 
that while the high frequency fundamentals in 
the methane and silane spectra appear to be 
less intense than the low frequency fundamentals, 
the opposite is here true, but aside from this the 
10.54 region in the germane spectrum resembles 
the 7.74 methane band a great deal, in fact much 
more so than the corresponding silane band 
which occurs as two overlapping bands. While 
here as in methane some regularity of rotationa 
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structure is discernible, especially on the high TABLE IV. Summary of po svat of the absorption bands 

frequency side where a separation of about 6.0 —— <ideniianemassiniaaaa Nannanaiaiaiainineaiaie 

cm~! seems to establish itself, the structure is) Onserven » COMPUTED v IDENTIFICATION 

quite definitely not the simple structure to be = — : 80cm em 

expected for so symmetrical a molecule. — : sadies Q) 
It is felt that these measurements as well as an 1988.7 sg 

the earlier ones on methane and silane emphasize aoe ~ " 

the importance of remeasuring the bands of all 3031.0 3044.1 ac 

4250 4220.2 2v2. 


the X-tetrahydrides (XH,) at much reduced 
temperatures, to establish whether the complex 
structure observed in many cases can be ascribed 
to overlapping bands arising from transitions 
between higher vibration levels or whether a 
more involved explanation is to be sought. We 
hope soon to be able to report on such measure- 
ments. 

When the spacings in the two optically active 
fundamentals vibrations v3 and vy have been 
determined, it is possible by the method of 
Dennison and Johnston* to compute the value 
for the moment of inertia of the molecule. This 
we have done also for the germane molecule 
taking the spacings between lines in the two 


’D. M. Dennison and M. Johnston, Phys. Rev. 47, 93 
(1935). 








fundamentals v3 and », to be 5.6 cm™ and 6.0 
cm~'!, respectively. These spacings lead to a 
value for the moment of inertia of the germane 
molecule of 7.0X10-* g-cm? and to the values 
1.37 X10-8 and 2.06 10-8 cm for the Ge-H and 
the H-H distances. It must be emphasized that 
the structure of v4 is so irregular and its spacing 
constant so doubtful that the above values are 
highly uncertain. 

For convenience we have restated in Table IV 
our identification of the bands observed in the 
spectrum of germane. 

We wish to express our gratitude to Professor 
Warren C. Johnston for the sample of germane 
gas which he so kindly presented to us. 
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The Infrared Absorption Bands of Methane 


ALVIN H. NIELSEN* AND HARALD H. NIELSEN, Mendenhall Laboratory of Physics, Ohio State University 
(Received September 30, 1935) 


The two fundamental vibration-rotation bands v3 and v4 
of methane have been remeasured with greater resolution 
than that available in the earlier investigation of Cooley. 
Whereas the band at 3.34 showed no essentially new 
details, the fundamental band at 7.74 was found to possess 
a somewhat more complex rotational structure than was 
observed by Cooley. Many of the lines given in Cooley's 
curve are here resolved into groups of closely spaced com- 


ECENT measurements on the absorption 
spectrum of silane! in the infrared have 
shown that the rotation lines are not single 
sharp lines, but are accompanied by satellite 


* Honorary Fellow in Physics, Ohio State University. 
1W. B. Steward and H. H. Nielsen, J. Chem. Phys. 2, 
712 (1934); Phys. Rev. 47, 828 (1935). 


ponents. The separation between the components within 
a group appears to increase with the ordinal number of the 
line counted outward from the center of the band. In 
addition to these measurements, the two combination 
bands », +4 at 4218 cm™ and »3+»,4 at 4315 cm™ have been 
at least partially resolved and appear to have line spacings 
of about 5.3 cm™ and 13.5 cm™, respectively. 


components. Because of the apparent similarity 
between the methane and silane molecules it has 
seemed worth while to remeasure the known 
absorption bands in methane under the im- 
proved conditions made available since the in- 
vestigation by Cooley? was made. 

2 J. P. Cooley, Astrophys. J. 62, 73 (1925). 
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Fic. 1. Absorption pattern of methane in the 3.34 region. 


For the measurements reported here, a prism- 
grating spectrometer with a collimating mirror 
of two-meter focal length was used. There were 
available three echelette gratings ruled by R. W. 
Wood, one with 3600 lines per inch for the 3.5y 
region, one with 4800 lines per inch used in 
second order for the 2.3u region, and one with 
2000 lines per inch for the 7.0u region. The 
methane gas was prepared by dropping distilled 
water on chemically pure aluminum carbide. 
To guard further against impurities the gas was 
fractionally distilled into the evacuated absorp- 
tion cells to atmospheric pressure. For the 3.3u 
and 7.74 bands a cell of 2 cm length was used, 
and for the 2.34 bands the cell was 6 cm in 
length. These cells were fitted with windows of 
polished rocksalt. 


THE 3.34 REGION 
The 3.34 fundamental band of methane, 13, 
is shown in Fig. 1. The curve is plotted in galva- 


TABLE I. The frequencies and wave-lengths of the lines in the 
3.3 band of methane. 











WAVE- WAVE- 
LINE FREQ. LENGTH | LINE FREQ. LENGTH 
No. (cm~') (mu) No. (cm~!) (u) 
—13 2886.1 3.4649; +2 3039.0 3.2906 


—12 2896.7 3.4522} +3 3048.6 3.2802 
—11 2907.0 3.4400 | +4 3058.1 3.2700 
—10 2917.7 3.4274} +5 3067.5 3.2600 
—9 2927.8 3.4155| +6 3077.0 3.2499 
—8 2938.1 3.4036 | +7 3086.3 3.2401 
—7 2948.3 3.3918; +8 3095.4 3.2306 
—6 2958.6 3.3800 | +9 3104.7 3.2209 
—5 2968.8 3.3684 | +10 3113.6 3.2117 
—4 2979.1 3.3567 | +11 3122.6 3.2025 
—3 2989.2 3.3454 | +12 3131.6 3.1933 
—2 2999.7 3.3367 | +13 3140.3 3.1844 
0 3019.6 3.3117 | +14 3149.2 3.1754 
+1 3029.1 3.3013 | +15 3157.8 3.1668 











nometer deflections against a scale of frequencies 
in cm~ and a scale of wave-lengths in yu. The 
band, while here completely resolved, reveals 
nothing essentially different from the earlier 
curve of Cooley except for the very weak back- 
ground observable between the lines in Fig. 1. 
Concerning this background it must be said 
that while it can be reproduced in detail and can 
in no way be ascribed to observational errors, 
it is difficult to identify it definitely as a part of 
the methane spectrum since there is in this 
region a small amount of atmospheric water 
vapor absorption. Measurements have, however, 
been made on this band with the cell in the beam 
and with the cell out of the beam, and when 
the data were plotted as percent absorption this 
same background persisted. It was also found 
impossible to make any correlation whatever 
between the background and the fluctuation of 
the deflections due to atmospheric absorption. 
The slit widths used for measurement in this 
region included about 0.8 cm and galvanometer 
deflections were read at intervals of about 0.8 
cm. Table I gives the new frequency and 
wave-length positions of the lines in the band, 
and numbers corresponding to those in Cooley's 
curve have been assigned. 


THE 7.74 REGION 


In Fig. 2 we reproduce »4, the 7.74 funda- 
mental of methane as given by our measure- 
ments. This curve is plotted in percent absorption 
against a scale of frequencies in cm~! and wave- 
lengths in yw, readings having been taken with 
the cell in the beam and with the cell out of the 
beam to compensate for any atmospheric lines 
in this region. Galvanometer deflections were 
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Fic. 2. Absorption pattern of methane in the 7.74 region. 


read at intervals of about 0.4 cm™ and the slits 
included about 0.45 cm~!. Several of the lines, 
however, were studied separately with some- 
what narrower slits and with the readings taken 
at correspondingly smaller intervals. Cooley’s 
curve quite definitely suggests a more complex 
structure, and this is just what is shown by our 
measurements in Fig. 2. The nature of the effect 
is best seen in the high frequency side of the 
band, the low frequency side being further 
distorted by the very broad Q branch which 
here is apparently partially resolved. On the 
high frequency side no essential difference be- 
tween Cooley’s curve and Fig. 2 is observed 
until the seventh line which we have resolved 
into two components lying 0.8 cm~ apart. The 
components of the following lines are successively 
more widely separated. It is difficult to compare 
Cooley’s curve with Fig. 2 on the low frequency 
side since Cooley’s curve extends only to the 
fifth line. It will be seen, however, that the 
positions in our curve corresponding to lines —2, 
—3, —4 and —5 contain 3, 4, 2 and 2 com- 
ponents, respectively, some of which are probably 
due to the resolved structure of the very broad 
Q branch. Table II gives the frequency and wave- 
length positions of the lines in the band numbered 
as in Cooley’s work except that the components 
of the lines have been given the assignments 
a, b, c, ete. 


THE 2.34 REGION 


The combination bands »;+ 74 and v3+ vy, at 
2.3u, at least partially resolved, are shown in 
Fig. 3. It will be remembered that none of the 


fine structure of these bands was resolved in 
the earlier work. The curve in Fig. 3 is plotted 
in galvanometer deflections against a scale of 
frequencies in cm~ and wave-lengths in yu, the 
deflections having been read at intervals of 
0.9 cm-! by using slit widths which included 
about 0.8 cm~. With the Q branches separated 
by only 97 cm™ the R branch of v;+», and the P 
branch of v3+¥v, overlap to such an extent that 
no estimate of the line spacings of these bands 


TABLE II. The frequencies and wave-lengths of the lines in the 
7.74 band of methane. 











WAVE- WAVE- 
LINE FREQ. LENGTH LINE FREQ. LENGTH 
No. (em~!) (u) No. (cm™~!) (u) 
—8c 1253.0 7.9808 f 1302.4 7.6781 
b 1254.4 7.9719 e 1304.1 7.6681 
a 1256.1 7.9612 d 1305.1 7.6622 
-—7d 1259.1 7.9422 c 1306.1 7.6564 
c 1260.3 7.9346 b 1306.8 7.6523 
b 1261.5 7.9271 a 1307.6 7.6476 
a 1262.6 7.9202 
—6e 1265.2 7.9039 1 1312.2 7.6207 
d 1267.4 7.8902 2 1317.5 7.5901 
c 1268.5 7.8833 3 1322.9 7.5591 
b 1270.3 7.8722 4 1328.1 7.5295 
a 1271.1 7.8672 5 1333.3 7.5001 
—5b 1275.1 7.8431 6 1338.2 7.4727 
a 1277.3 7.8290 7a 1342.7 7.4476 
—4d 1278.4 7.8223 b 1343.5 7.4432 
c 1279.5 7.8156 8a 1347.4 7.4217 
b 1281.7 7.8021 b 1348.7 7.4145 
a 1283.1 7.7936 9a 1352.0 7.3964 
—3d 1286.6 7.7724 b 1353.7 7.3871 
Cc 1287.7 7.7658 10a 1356.4 7.3724 
b 1288.7 7.7598 b 1358.9 7.3752 
a 1289.6 7.7543 lla 1360.8 7.3484 
—2c 1293.0 7.7340 b 1362.3 7.3404 
b 1294.9 7.7226 Cc 1364.9 7.3261 
a 1295.7 7.7184 
—ib 1298.1 7.7036 
a 1300.3 7.6905 
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Fic. 3. Absorption pattern of methane in the 2.34 region. 


can be made from this region. However, certain 
lines stand out more prominently than others in 
the P branch of »,+», and the R branch of 
v3+v, and it seems reasonable to assign a spacing 
of about 5.3 cm and 13.5 cm™ to these two 


TABLE III. The frequencies and wave-lengths of the lines in 
the 2.3 bands of methane. 











WAVE- WAVE- 
LINE FREQ. LENGTH | LINE FREQ. LENGTH 
No. (cm~!) (us) No. (cm~!) (u) 
1 4143.4 2.4135 31 4265.8 2.3442 
2 4150.4 2.4094 |} 32 4269.9 2.3420 
3 4152.8 2.4080 | 33 4272.0 2.3408 
4 4154.4 2.4071 | 34 4275.4 2.3390 
5 4157.1 2.4055 | 35 4280.8 2.3360 
6 4159.9 2.4039 | 36 4287.2 2.3325 
7 4163.0 2.4021} 37 4293.1 2.3293 
8 4166.2 2.4003 | 38 4298.2 2.3266 
9 4171.7 2.3971 | 39 4301.6 2.3247 
10 4174.5 2.3955 | 40 4304.1 2.3234 
11 4178.1 2.3934 | 41 4307.5 2.3215 
12 4181.3 2.3916} 42 4314.4 2.3178 
13 4183.7 2.3902 | 43 4320.4 2.3146 
14 4187.7 2.3879 | 44 4323.8 2.3128 
15 4189.7 2.3868 | 45 4328.5 2.3103 
16 4195.7 2.3834 | 46 4332.4 2.3081 
17 4200.5 2.3807 | 47 4339.6 2.3044 
18 4206.2 2.3774 | 48 4345.9 2.3010 
19 4210.2 2.3752 | 49 4349.1 2.2993 
20 4212.9 2.3737 | 50 4355.9 2.2957 
21 4217.5 2.3710} 51 4360.5 2.2933 
22 4223.2 2.3679 | 52 4365.6 2.2906 
23 4225.7 2.3665 | 53 4378.8 2.2837 
24 4229.4 2.3644 | 54 4382.8 2.2816 
25 4233.9 2.3619 55 4393.0 2.2763 
26 4239.2 2.3589 | 56 4399.2 2.2731 
27 4244.2 2.3562 | 57 4406.4 2.2694 
28 4248.7 2.3537 | 58 4414.5 2.2653 
29 4253.7 2.3509 | 59 4419.5 2.2627 
30 4257.4 2.3489 | 60 4429.9 2.2574 











bands, respectively. Table III gives the fre- 
quencies in cm= and the wave-lengths in y» of 
the lines resolved in this region. For identifica- 
tion arbitrary numbers have been assigned to the 
lines in the two bands. 


DISCUSSION OF EXPERIMENTAL RESULTS 


We do not, in this communication, propose to 
give any adequate explanation for the complex 
structure observed at 7.74. The possibility sug- 
gests itself that it may be due to upper stage 
bands; i.e., to transitions between higher vibra- 
tion states. It would, however, seem unlikely 
that these would be observed as anything but 
second-order effects in this region of the spec- 
trum. Another explanation is the one suggested 
in connection with a similar effect observed in 
silane, that this structure is in some way inti- 
mately associated with the convergence of the 
lines in the bands. This issue can only be defi- 
nitely settled, as was earlier suggested, by re- 
measuring the bands with the gas at greatly 
reduced temperatures. This we hope to do in 
connection with similar measurements on the 
bands of silane and germane. 

May we take this opportunity to express our 
appreciation of a gift from the International Salt 
Company of Scranton, Pennsylvania, of several 
pieces of optical rocksalt used in making windows 
for our absorption cells. A grant-in-aid from the 
National Research Council is also acknowledged 
with gratitude. 
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The Interaction Between Vibration and Rotation for Symmetrical Molecules 
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The explanation first proposed by Teller of the anomalous 
fine structure of the infrared bands of symmetrical mole- 
cules is discussed. The fine structure is closely related to an 
internal angular momentum of magnitude ¢h/2z which is 
due to the vibration of the molecule and arises from the 
degenerate character of the motion. A simple derivation is 
given for the spacing constant [(1—¢)/C—1/A Jh/4x? of 
axially symmetric molecules and for the spacing constant 
(1—¢)h/4x°A of tetrahedral molecules. A detailed calcula- 
tion is made of the internal angular momenta to be asso- 
ciated with the | frequencies v2 and v, of the axial molecule 
YX;. The resulting ¢’s are found as functions of the mo- 
ments of inertia and of the potential constants. The sum 
é2+¢4 is shown to have the value C/2A —1. These results 
are applied to the molecules NH; and ND; and the line 
spacings of v2 and v4 are computed. The axial molecule 
ZYX; is treated and it is proved that the sum of the ¢’s 
characterizing the three | bands ve, vs, and vs is equal to 
C/2A. The moment of inertia C of the methyl halides is 
computed and found to be, 5.61, 5.35, 5.44, and 5.44 « 10° *” 
for methy] fluoride to methyl iodide, respectively. The error 
is estimated to be around 5 percent. A calculation is made 
of the ¢’s to be associated with the overtones of axial 
molecules possessing threefold symmetry. It is found that 


the ¢ appropriate for the overtone 2»; of a | frequency is 
—2¢; while the ¢ for 3»; is ¢; itself. The combination of 
two | frequencies v2+ 4 is next treated and the ¢ is proved 
to be —({2+s)=1-—C/2A. An application of these 
formulae to the observed spacings of the overtone bands of 
NH; and CH;Cl results in a very satisfactory agreement. 
Expressions are obtained for the two ¢’s which determine 
the line spacings of the active fundamentals v; and », of 
the tetrahedral molecule YX4. It is shown that ¢3+¢4=}. 
The positions of the fundamental bands of methane to- 
gether with the line spacings of v3 and v4 yield the moment 
of inertia A =5.47 X10-*° together with the five potential 
constants describing the molecule. A good agreement is 
found between these constants and those obtained by 
Ginsburg and Barker from the spectrum of methyl deu- 
teride. The ¢’s suitable for the overtones of YX, are next 
treated. It is found that the ¢ to be used with 2»;, 1=3, 4 
is —¢; while the ¢ belonging with the combination band 
vst+vs is —}3(f3+¢4) = —}. The observations on the over- 
tones of methane are in very good accord with these 
expressions. Less satisfactory agreement is obtained in the 
case of the silane spectrum where it appears that the higher 
order perturbation terms play a larger role. 





§1. INTRODUCTION 


URING the past several years, great prog- 

ress has been made in interpreting the 
spectra of polyatomic molecules. The molecule 
has been represented by a system of mass points 
whose equilibrium positions constitute the con- 
figuration of the molecule. It is assumed that, in 
zeroth approximation, the mass points are bound 
by elastic forces and that throughout their 
vibratory motion they remain essentially in the 
neighborhood of their equilibrium configuration. 
From this model it has been possible to under- 
stand the spectra of many molecules and to com- 
pute such physical constants as the potential 
energy function and the dimensions of the sys- 
tem. However, one grave difficulty has appeared 
in that the simple theory predicts that for a given 
molecule all the infrared bands which correspond 
“to a particular type of vibration should show the 
same fine structure. This prediction is fulfilled 
for many types of vibration, for example, the || 
bands of ammonia. However, it is flagrantly 
violated for others, notably the bands of methane 


and the | bands of the methy! halides. 

A solution to this problem has recently been 
given by Teller! who has pointed out that those 
vibrations which show anomalous fine structure 
are degenerate and may possess an angular 
momentum arising from the vibratory motion 
itself. The interaction between this angular mo- 
mentum and that due to the rotation of the 
molecule will just account for the observed 
effects. The purpose of the present paper is to 
extend Teller’s theory and to make detailed cal- 
culations for several molecular models. 

The algebra involved in the computations is 
often rather complicated and may hide the sim- 
plicity of the underlying ideas. The vibratory 
motion of the nuclei may be regarded as a motion 
taking place on a slowly rotating turntable cor- 
responding to the rotation of the molecule. This 
rotation of the system introduces a Coriolis force 
into the vibrational motion which in turn makes 
a contribution to the energy. In general, the con- 


1E. Teller, Hand- und Jahrbuch d. chem. Physik (1934), 
Vol. 9, p. 125. 
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tribution is very small (of second order) since the 
motion of each particle relative to the turntable 
is, to a high approximation, simple harmonic 
motion along a straight line. Thus the average 
Coriolis force is zero. However, for a degenerate 
vibration, which corresponds to the existence of 
two or more harmonic motions of equal frequency 
at right angles to each other, this is by no means 
the case. A suitable choice of amplitudes and 
phases renders the motion for each particle circu- 
lar or at least elliptical. Here the average Coriolis 
force does not vanish but has the value +kw 
where & is a proportionality factor and w is the 
angular velocity of the turntable. The plus and 
minus signs relate to whether or not the circular 
motion of vibration is in the same direction as the 
rotation of the turntable. We shall now compute 
in detail the first order change in energy caused 
by the rotation. 

It is convenient to introduce a moving coordi- 
nate system x, y, whose axes are defined by the 
equilibrium positions of the nuclei and which 
coincide with the principal axes of inertia A, B, C. 
If the molecule possesses degenerate frequencies, 
there will be in general an angular momentum 
produced by the vibrational motion itself. The 
components of this internal angular momentum 
p, along the axes x, y, 2 will be designated by 
Pz, Py, pz. The eigenvalues of the internal angular 
momentum are not necessarily integral multiples 
of h/2x since p relates to only part of the system, 
namely, the vibration, and does not take account 
of the rotation of the molecule. Following Teller 
we shall introduce ¢ as the proportionality factor. 
In molecules having an axis of symmetry, such 
as the methyl halides, ammonia, etc., p lies along 
this axis and has the magnitude ¢h/27. For mole- 
cules having tetrahedral symmetry (methane) 
the vibrations in which we are interested possess 
spherical symmetry and consequently, 


pe=M(I+1) eh? /49?. 


The quantity ¢ is essentially dependent upon 
the vibrational properties of the molecule, that 
is, upon the relative dimensions, the masses, and 
the force constants. It has different values for the 
various normal modes of vibration. For this rea- 
son the individual bands of the same molecule 
may show different fine structure spacings. A 
knowledge of the {’s as determined from the fine 


structure enables us to obtain information regard- 
ing the force constants of the molecule. 

It is now necessary to consider the total angular 
momentum of the molecule P which includes the 
contributions due to the rotation and to the vi- 
bration. It has the components P,, P,, P, along 
the moving axes and must be quantized. The for- 
mula for the rotational energy of a rotator in- 
volving an internal flywheel has appeared many 
times in the literature and has the form, 


Trot = (P:—pz)*/2A + (Py — p,y)?/2B 
+(P,—p.)?/2C. 


This expression has recently been carefully re- 
examined by Eckart,? Van Vleck,’ and others, 
and appears to be entirely correct. 

The rotational energy levels of the symmetric 
rotator, A=B, C may be readily obtained. The 
internal angular momentum for such a molecule 
must lie along the symmetry axis C and hence 
p:=p,=0 and p.=+), the double sign indicat- 
ing whether ? is parallel or antiparallel to the C 
axis. We may write the rotational energy 


Hrot = (P2+P3+P 2)/2A 
+P 2(1/2C—1/2A)¥ pP./C+p*/2C. 


The Hamiltonian may be diagonalized by 
using the representation for P such that P, has 
only diagonal elements. P,=Kh/2x where K 
is a positive or negative integer. The total angular 
momentum is of course diagonal. 


P?2+P7+P2=J(J+1)h?/4r°. 
Hence 
Wr = J(J+1)h?/82°A 
+(1/C—1/A)K*h?/82? ¢Kh?/47°C. 


The last term p*?/2C may be omitted since it does 
not depend upon the rotation and hence may be 
absorbed into the vibrational energy. 

In order to find the actual frequencies of the 
fine structure lines we must make use of the selec- 
tion rules. These may be obtained by the rather 
laborious process of examining the matrix ele- 
ments of the electric moment. However, they 


2C. Eckart, Phys. Rev. 47, 552 (1935). 
3 J. H. Van Vleck, Phys. Rev. 47, 487 (1935). 
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may be easily understood from the following ar- 
gument. The electric moment depends upon the 
positions of the nuclei. In the coordinate system 
x, y, z which moves with the molecule, the electric 
moment executes circular motion about the axis 
of symmetry. Suppose that the direction of rota- 
tion of the electric moment is in the same sense 
as the direction of the angular momentum of vi- 
bration p. Now if p and P, are in the same sense 
(in which case we use the upper or minus sign in 
the expression for W,..) the motion of the electric 
moment may be described as a rotation about the 
symmetry axis with a frequency vo+ x and a pre- 
cession about the P axis with a frequency vy. 
Thus, if the vibrational quantum number goes 
from 0-1 we have K-K+1 and J—>J-1, 
J—J+1, or JJ as the possible transitions. In 
the case where p and P, are in opposite senses, 
the frequency of rotation is clearly vp) — vx and the 
selection rules are, vibrational number 0-1, 
K-—K-—1 and J>J-—1, J->J+1, or J>J. We 
must here use the plus sign in Wrot. 

In the 1 bands of axial molecules it usually 
happens that only the so-called zero-branch lines 
for which J—J are observed. The frequency of 
the fine structure lines under either of the above 
sets of rules (p parallel or antiparallel to P.) is 
then 


v=v—(1/C—1/A)h/8r? 
+[(1-—¢)/C—1/A ]Kh/4r’, 
K=0, +1, +2, ->>. 


In deriving this formula it was assumed that 
the direction of rotation of the electric moment 
was parallel to the angular momentum ~. How- 
ever, it may occur that the rotation of the elec- 
tric moment and the internal angular momentum 
are antiparallel. In this case the selection rules are 
just reversed and the spacing of the zero-branch 
lines may be shown to be [(1+¢)/C—1/A] 
X (h/4x*). It is convenient to retain only the first 
formula and to consider that ¢ is positive or nega- 
tive according to whether the rotation of the 
electric moment and the vibrational angular 
momentum have the same or opposite senses. 

Finally it must be remarked that we have here 
considered a transition from the normal state, 
where there is no vibrational angular momentum, 
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and consequently ¢=0, to an excited state ¢=¢. 
These bands are the ones most frequently ob- 
served experimentally. However, transitions may 
occur between two excited levels, each of which 
has a ¢ different from zero. In this case the above 
formula for the fine structure lines would be 
altered. 

The second class of molecules to be treated in 
this paper is the group of molecules which have 
tetrahedral symmetry. For these all three mo- 
ments of inertia are the same. The Hamiltonian 
may be written in the simple form 


Hot = P2/2A — (P+ p)/A+p?/2A. 


The first term becomes J(J+1)h?/82°A in vir- 
tue of the diagonal nature of the total angular 
momentum P?, while the third term may be ab- 
sorbed into the vibrational energy. The second 
term, which contains the scalar product of the 
total angular momentum P and the internal an- 
gular momentum #, represents the interaction 
which we wish to evaluate. 

The fundamental frequencies of the methane 
molecule which are active in the infrared are to 
be associated with a transition from the normal 
state where there is no internal angular momen- 
tum to the first excited state of a three-dimen- 
sional isotropic oscillator. This latter state 
possesses an internal angular momentum p which 
may have the three eigenvalues +¢{h/2z7 and 0. 
Thus the interaction term J may be written 
I= —(II-2)th?/42°A where II is a vector having 
eigenvalues J, J—1, --- —J and z has the values 
1, 0, —1. The expression for the eigenvalues of 
(11-7) is well known in spectroscopy and is 


2(II- r) =J(J+1)+2—R(R+1). 


R is a vector which, when compounded with 7, 
will give the total momentum II. Evidently it 
may take on the values J—1, J, and J+1 thus 
yielding the following three eigenvalues for 
(II-x); J+1, 1, and —J. Now the vector R cor- 
responds to a rotation of the molecular frame- 
work and the three values for (II-7) are to be 
correlated with the three possibilities, p parallel 
to P, p perpendicular to P, and antiparallel to 
P, respectively. These considerations furnish a 
key to the selection rules and show that in mak- 
ing a transition from the normal to the excited 
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state we must use the first, second, or third 
eigenvalue of (II-7) according as J—1—J, 
J—J, or J-J—1. We thus obtain the following 
formulae for the fine structure lines 

Positive branch v= v9—h{/42°A+Jh(1—¢)/47°A 

J-1—-J 
Zero branch 

J—J 
Negative branch v= v9p—h{/42°A —Jh(1—¢£)/47°A. 

J—J-—-1. 

It thus appears that the form of the vibration- 
rotation bands has not been changed by the 
introduction of an internal angular momentum. 
However the actual spacing of the fine structure 
lines is altered by the factor (1—¢). These re- 
sults, which have also been obtained by Teller, 
we believe to be correct only up to the first-order 
approximation here treated. In higher approxi- 
mation where the particular symmetry of the 
molecule must be taken into account it seems 
probable that the fine structure lines will be di- 
vided into multiplets.‘ It is hoped shortly to 
make calculations on this effect of which there 
are already some experimental indications. 

The formulae which have been developed 
show that the line spacing depends essentially 
upon a quantity ¢ without a knowledge of which 
the experimental data cannot be interpreted. 
th/2x represents the angular momentum due to 
the vibratory motion of the nuclei relative to a 
coordinate system defined by the equilibrium 
positions of the nuclei. This internal angular 
momentum may always be computed as a func- 
tion of the force constants, the masses, and the 
relative dimensions of the molecule, but in many 
cases the actual calculation becomes very labori- 
ous and involved. There exists, however, an im- 
portant theorem discovered by Teller! which 
states that the sum of the ¢’s of all the first ex- 
cited states belonging to a given symmetry class 
is independent of the force constants and is a 
function only of the masses and the relative di- 
mensions of the molecule. The quantity 2¢ is thus 
relatively simple to compute since in making the 
calculation we may use any type of force field, 
in particular, such forces that the motion of the 
system may be readily visualized. Immediate use 
may be made of >¢ since the sum of the fine 


v=vy—ht/4r°A 


‘See the concluding portion of the present paper. 
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structure spacings of all the fundamental bands 
of a given category depends only upon Y¢ and 
not upon the individual ¢’s. This point will be 
illustrated later in the present paper. 


$2. AXIALLY SYMMETRIC MOLECULES 


We begin our discussion of molecules having an 
axis of symmetry by considering the system Y Xs, 
of which ammonia is a representative. It is well 
known that the regular pyramid YX; possesses 
four normal frequencies which are separated into 
two symmetry classes, the two || frequencies and 
the two | frequencies. Since these belong to 
different symmetry classes the Hamiltonian falls 
into two parts and there exist no interaction 
terms. It is only the | frequencies which are de- 
generate and which may exhibit an internal an- 
gular momentum. Consequently we shall treat 
them exclusively in what follows. In separating 
the |! and | vibrations it is only necessary to in- 
troduce symmetry coordinates; coordinates which 
are capable of specifying all possible configura- 
tions of the system which belong to the desired 
symmetry class. 

We shall begin by obtaining the kinetic energy 
and we will make use of the theorem that the 
kinetic energy of any system may be regarded as 
the sum of two terms, the energy relative to the 
center of gravity of the system and the transla- 
tional energy of the center of gravity. We shall 
divide the YX; molecule into two parts, the Y 
atom and the X; group. The energy may thus be 
written as the energy of the Y atom, the trans- 
lational energy of the center of gravity of the X; 
group and the energy of the X; group relative to 
its center of gravity. The method of calculation is 
made clearer by Fig. 1. The dotted horizontal line 
is the equilibrium position of the figure axis. The 
center of gravity of the X atoms is shown by the 
letter G and the quantity y; is the sum of the dis- 
placements of the Y atom and G. The inclination 
of the plane of the X atoms to the symmetry line 
is indicated by the angle @. 

The kinetic energy of the whole system is then 


T= eye + 316+ To, 


where yu is the reduced mass of the Y atom and G, 
that is, ~x=3mM/(3m+M); m and M being the 
masses of the X and Y atoms, respectively. The 
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second term represents the energy of tipping of 
the X; triangle and J is the moment of inertia of 
the triangle about an axis through G lying in the 
plane of the triangle. It will be more convenient 
to use the substitution /=C,/2 where C is the 
moment of inertia of the X3 atoms, and hence of 
the whole YX; molecule about the figure axis. 
7) denotes the kinetic energy of the X atoms 
within their plane and relative to G. We shall 
presently evaluate 7); the only point to consider 
for the instant is that there can be no interaction 
terms between 7» and the tipping term }/@ since 
the displacements of the particles for the two 
motions are always perpendicular. 

The first two terms of 7 may be combined be- 
cause of the condition that the angular momen- 
tum of the system perpendicular to the symmetry 
axis must remain zero throughout the motion. 
(This is equivalent to the demand, explained 
earlier, that we examine the motion in a moving 
coordinate system x, y, 2, which is fixed by the 
equilibrium configuration of the system.) Thus, 
uhy,= 16 where h is the height of the pyramid. 
But uwh?+7=A, the moment of inertia of the 
molecule perpendicular to the symmetry axis. 
By substitution we obtain for the first two terms 
of the kinetic energy »Ay,;’/C. It has appeared 
simplest to develop the kinetic energy as a func- 
tion of y; and one might continue to treat the 
motion using this coordinate. Such a treatment 
would suffer a disadvantage, however, in that y; 
does not specify a geometric property of the sys- 
tem, since the actual geometric configuration is 
determined both by y; and by the tipping of the 
X; triangle. That is to say, it involves the mo- 
mentum relation and hence implies that the posi- 
tions of the particles are a function both of y; and 
of the masses. This difficulty, which becomes par- 
ticularly important when we examine isotope 
effects, may be obviated by choosing a geometric 
coordinate, one which characterizes a geometric 
property of the system and is independent of the 
masses. For example we may use the perpendicu- 
lar distance between the Y atom and the normal 
to the X; plane through G. Let this be called y. 
Clearly 


y=ythd or y=(14+2ph?/C)y,=(24/C)y1 


by reason of the momentum relation. 
Up to the present we have considered a dis- 
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placement y perpendicular to the symmetry axis 
and lying in the plane of the paper (Fig. 1). Let 
the corresponding displacement perpendicular to 
the plane of the paper be x. Making the indicated 
substitutions we obtain 


T= 3(uC/2A)(#+y*) +To. 


The kinetic energy 7» of the X atoms within their 
plane and relative to their center of gravity G has 
appeared many times in the literature’ since it is 
the kinetic energy for a triatomic molecule. We 
begin by expressing it as a function of qi, go, 
and g3, the changes in distance between pairs of 
atoms. 


To= 3m { (5/9) (G1? +42? +43") 
— (2/9) (Gide +4193 + 429s) }. 


A linear transformation may now be made 
which diagonalizes 7) and at the same time intro- 
duces symmetric coordinates. Let 


£=(qi—q2)/3', n=3(gitge—29s), 


$= 3(qitg2+9s), 


in which case 7p>=3m(#+7+). The displace- 
ments corresponding to é, 7, and ¢ are shown in 
Fig. 2 and it is at once evident that £ and » form 
a degenerate pair capable of describing a | vi- 
bration. ¢ belongs to the other symmetry class 
and hence for our present purpose it may be dis- 
regarded. 

The motions corresponding to ~ and 7 are 
subject to the condition that the linear and angu- 
lar momenta of the system must remain zero. 
This serves to determine both the direction and 
magnitude of the actual space displacements of 
the particles which may be characterized as fol- 


5 N. Bjerrum, Verh. d. D. Phys. Ges. 16, 737 (1914). 
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lows. (1) The displacement of any individual 

atom in the £ motion is perpendicular to its dis- 

placement in the 7 motion. (2) The displacements 

of all the particles are the same and are given, 

respectively, by d;=£/3! and d;=n/3' where 7 is 

the number of the atom and runs from 1 to 3. 
The total kinetic energy thus becomes 


T=3mla#+H)+E+#), 


where a=uC/2mA. 

The potential energy may be written as a 
general quadratic form involving the arbitrary 
coefficients a, b, and c. 


V= 3 la(eP+y’) +0(P +9?) +2c(xé+yn)}. 


Although there are ten coefficients entering a 
general quadratic form containing four variables, 
the above expression is the most general one 
which may be used for the potential energy. This 
is due to the axial symmetry of the molecule, 
which insures first, that x and y and also £ and n 
must enter on the same footing, and second, 
that the coefficients of such terms as xy, xn, etc., 
must vanish. 

The Hamiltonian clearly may be divided into 
the sum of two identical Hamiltonians, the first a 
function of x and é and the second of y and 7». 
Either of these may be treated by the method 
of normal vibrations, which yields the following 
relations between the X’s, where the normal 
frequency v;=),3/27: 


m(Ao+A4) =a/at+b, mod == (ab—c*)/a. 


This result is in agreement with the expres- 
sions obtained by J. E. Rosenthal® for the fre- 
quencies of the YX; molecule although the nota- 
tion is somewhat different. We reproduce it here 
merely because it is one of the steps in the calcu- 
lation of the internal angular momentum factors 
fo and &%. 
~ 6 J. E. Rosenthal, Phys. Rev. 47, 235 (1935). The rela- 


tions between our constants and those of Dr. Rosenthal 
are A =a, 3D=5), and 3E*=c?. 
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It is now necessary to transform to the normal 
coordinates x2, Xs, Ye, and y, whereby 


IT = H1,(x2, x4) +171(y2, ya), 
IT \(xe, X4) == B( PrP t+ pr +roxe? + A4x,’). 


The normal coordinates x2, x4 are related to x 
and & by means of a linear transformation whose 
coefficients are proportional to the first minors of 
the determinant which gave the \’s. We find 


x=TA.x;, §=ZTB xi, i=2, 4, 
where A;=(md;—b)/[ma(md;—b)?+mce?}}, 
B;=c/[ma(md;— b)?+me}}. 


Identical expressions connect ye and y, with y 
and ». 

The angular momentum of the entire system 
might be determined by finding the angular mo- 
mentum of each particle and summing. However 
it will be more convenient to divide the molecule 
into two parts, the Y atom and the X; group. 
Now, as is well known, the angular momentum of 
any group of particles (say X;) is equal to the 
angular momentum of the particles relative to 
their center of gravity G, plus the contribution 
due to the center of gravity. Thus Y and G 
together, which constitute a two particle group, 
will possess an angular momentum u(x — yids) 
=yu(C/2A)?(xy—yaz). The angular momentum of 
the Xs group relative to G may be found by 
considering the momentum of any one particle 
and multiplying by 3. We obtain —m(§)—n€). 
This expression is preceded by a minus sign since 
it is oppositely directed from the momentum of 
Y and G. The reader may easily verify this state- 
ment by drawing several phases of the motion 
illustrated in Fig. 2. The momentum due to Y 
and G has the same direction as the rotation of 
the electric moment, since Y and G are the 
carriers of the electric charges. 

We may now transform to the normal coordin- 
ates and (xy — yd) becomes A 2?(x2t/2— Vote) 
+A? (xis — Ysts) + AA a(X2is — Yate t+ Xue — Vods). 
The eigenvalues of these terms may be found 
very easily when we recall that the wave func- 
tions are the usual Hermitian orthogonal func- 
tions in the variables 2r(v;/h)!x,; and 2r(v;/h)*yi. 
We obtain A2’h/22 or Agh/2x for the first ex- 
cited states of ve and », respectively. Thus the 
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contribution to ¢; from the motion of Y and G is 
u(C/2A)*A 2. From similar reasoning one obtains 
the remaining part of ¢;. which is due to the 
momentum of the X; particles relative to G. 

Adding these two terms and substituting for 
A, and B; we have finally 


(:=(C/2A —6,;)/(1+46,) 
6;=2mAc?/uC(mdr,;—b)?. 


i=2, 4, 
where 


The sum (2+; may now be determined by 
direct substitution. The resulting complicated 
expression may be reduced by means of the rela- 
tions between 2 and dy to the simple form 
Y¢=C/2A—1. The fact that LF turns out to be 
independent of the force constants a, }, and ¢, is 
both an example of the general theorem relating 
to S¢ and a verification that our calculations 
have been free from algebraic errors. 

These formulae may be used in the analysis of 
the spectra of NH; and ND3;. The experimental 
data are as follows. 

The long wave perpendicular band », of NH; 
has been mapped by E. F. Barker’ and has its 
center at 1631 cm~. It shows an intense central 
region of absorption which is between 30 and 40 
cm wide. This no doubt corresponds to a 
clustering of zero branch lines whose spacing 
Av, may be estimated at around 3 cm. The spec- 
trum of ND; has been studied by M. V. Migeotte 
and E. F. Barker® who have found both of the 
fundamental perpendicular bands, ve=2556 and 
vy=1188. The short wave bands consist of a set 
of equally spaced lines (the zero-branch lines) 
with a spacing Ave=5.2 cm='. The long wave 
band », exhibits a strong central absorption cor- 
responding to a spacing of from 1 to 2 cm~. We 
wish to express our gratitude to Professor Barker 
and Dr. Migeotte for allowing us to reproduce 
here a portion of their unpublished data. 

We shall adopt the following procedure. By 
means of the observed band centers v5 of NH; and 
vo, vs Of NDs, the potential constants a, b, and c 
may be determined. From a knowledge of these 
constants, the ¢’s and finally the fine structure 
spacings may be computed and compared with 
the observed spacings. In making all these calcu- 
lations we require the ratio of the moments of 


‘ To be published in the near future. 
* M. V. Migeotte and E. F. Barker, Phys. Rev., in press. 
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TABLE I. Values of the ¢'s and line spacings. 
A $ g 











fe Ave Comp. Ave OBs. vs AvsComp. Avs OBs. 
NH; -—0.004 7.15 - —0.171 5.76 ~3 to4 
2.66 ~lto2 


ND; + .104 4.81 5.2 — .257 


inertia C/2A. To find this we shall use the line 
spacings of the pure rotation spectra of the am- 
monias. Wright and Randall® obtained for NH; 
that h/47°A =19.88 whereas Barnes, Benedict 
and Lewis!’ found for ND; that #/47°A =10.22 
cm~!. These data yield C’/2A =0.825 and 0.847 
for NH; and NDs, respectively. Unfortunately 
this ratio is extremely sensitive to small changes 
in the line spacings and appears to be determined 
with an accuracy of not better than 5 to 10 per- 
cent. The potential constants are then found to 
be a=9.42, b=4.12, c=+3.97 X10° dynes /cm. 
Here again the errors are considerable, amounting 
to at least 5 percent. This is due to two causes: 
first to the inaccuracy in C/2A; second, to the 
fact that we have used the observed positions of 
the fundamental bands rather than the normal 
frequencies which are demanded by the theory. 
The normal frequencies themselves may be found 
only when all the anharmonic constants are 
known. 

The results of the calculations on the ¢’s and 
the line spacings are shown in Table I. 

The agreement between the computed and 
observed Ay is not at all satisfactory although it 
appears to be best in the one case, Ave of NDs, 
where the experimental data are most precise. 
We feel that these results do not constitute any 
real test of the theory. They show rather that the 
Av are sensitive to the potential constants and to 
the ratio C/2A, and consequently it may be 
possible through them to obtain more exact values 
for the molecular constants. This must wait how- 
ever in the case of the ammonias until further 
measurements are available on the line spacings. 

The next class of molecules to be considered is 
the axial system ZYX; to which the methyl 
halides belong. There exist three perpendicular 
bands which are determined by six potential 
constants. The task of computing the individual 
¢’s in terms of the six constants becomes rather 


9N. Wright and H. M. Randall, Phys. Rev. 44, 391 


(1933). 
1 R. B. Barnes, W. S. Benedict and C. M. Lewis, Phys. 


Rev. 45, 347 (1934). 
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arduous and for the present we shall merely ob- 
tain an expression for 2¢ over the three perpen- 
dicular fundamental bands. Now 2X¢ is inde- 
pendent of the potential constants and therefore 
may be calculated by using any limiting force 
field. The result is then correct for all force fields. 

We have chosen the limiting field to be such 
that the YX; group is held together by very 
strong forces while the Z particle is very weakly 
bound to the YX; group. The three perpendicular 
frequencies may now be characterized by v2 and 
vs which correspond to a motion within the YX; 
group itself, the Z atom standing still, and yg, 
which is to be correlated with a vibration of the 
Z atom against a rigid pyramid YX;. In this 
latter mode of vibration, the Z atom moves in a 
path perpendicular to the symmetry line against 
the center of gravity G of the YX; pyramid. This 
motion is shown in Fig. 3. In order to preserve the 
constancy of the angular momentum the pyramid 
must tip through the angle @. 

In this extreme type of force field the frequen- 
cies vp and »; belong to a motion of the YX; group 
alone and hence we have {2+{,= C/2A’—1 where 
A’ is the moment of inertia of the YX; pyramid 
about an axis perpendicular to the symmetry 
axis and drawn through G. C is the moment of 
inertia of the YX; group, or, in fact, of the entire 
molecule ZY X3 along the symmetry axis. It re- 
mains now to compute {«. The total kinetic en- 
ergy for the motion v¢ is the sum of the energies of 
Z, of G and of the tipping of the pyramid. 


T= 3(miy?+A’'6), 


where y is the sum of the displacements of Z and 
of Gand yy; is the reduced mass p; = M,(Me2+3m)/ 
M,+M2+3m; Mi, M2 and m are the masses of a 
Z, Y and X atom, respectively. 

The coordinates y and @ are not independent 
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but are connected by the momentum relation 
u,by=A'O 


where 0 is the equilibrium distance between Z 


and G. Thus 
T =}(uiturb?/A’)y. 


This expression may be somewhat simplified by 
recalling that A’+y,5?=A, the moment of inertia 
of the whole molecule about an axis perpendicu- 
lar to the symmetry axis. We shall also introduce 
the coordinate x, corresponding to a displacement 
of Z relative to G, out from the plane of the paper 
in Fig. 3. 


T=}3(mA/A’)(? +9). 


The kinetic energy may however be written 
as a function of the normal coordinates x¢ and ¢. 


T=3(t68 +96) 
and therefore we have the relations 
x=(A'/p1A) x, and y=(A’, 111A) v6. 


The total angular momentum of the system 
about the symmetry axis is most easily obtained 
by decomposing it into two parts. First, consider 
the contribution due to the motion of Z and G. 
This is clearly y:(xy— ya). The second contribu- 
tion is that coming from the motion of the YX; 
pyramid relative to its center of gravity G. This 
becomes p2(x1j1—i41) where pe is the reduced 
mass 3mM2/(3m+ M2) and x, y, are the dis- 
placements of the Y atom relative to the center 
of gravity of the X; group, perpendicular to and 
in the plane of the paper (Fig. 3), respectively." 
The coordinate y,;=h@ where, as earlier, h is the 
height of the YX; pyramid. Hence, through the 
momentum relation y,; is connected with y. 
Thus y:=(u:bh/A’)y and the total angular mo- 
mentum of the ZYX; molecule becomes 


(mitme(uibh/A’)*) (xy — yz) 
or in terms of the normal coordinates 


1 Tn considering this second contribution we have treated 
the YX pyramid as two units, the Y atom and the X; 
group. It might be thought that there would be a third 
contribution to the angular momentum coming from the 
motion of the X; group relative to its center of gravity. 
This term is however zero since, relative to their center of 
gravity, the displacements of each X atom are always 
parallel to the symmetry axis. 
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[uitme(uidbh/A’)* \(A’/u1A)(xeije—Yeorts)- 


The eigenvalues of (xeije— vets) may be found 
by the method indicated in the calculations on 
the YX; molecules. For the first excited state 
they are +h/27. Therefore 


fo=([uitue(uibh/A’)? }(A’/piA). 
A simple calculation reduces this to the form 
f6=C/2A —C/2A’+1. 


The sum over all the ¢’s therefore becomes 
~¢=C/2A. We wish to emphasize again that, al- 
though this result has been obtained from a 
computation using a limiting and physically un- 
real force field, it will hold true for any actual 
force field. The conditions for its strict validity 
are, harmonicity of the motion and no resonance 
relations between the frequencies. It is clear that 
this method of the limiting fields may be used to 
determine the Y¢ for any axial molecule. 

The fundamental perpendicular bands of the 
methyl halide molecules have all been observed 
by Bennett and Meyer” and their fine structure 
measured. Since in the present theory ¢ and 
not the individual ¢’s themselves have been com- 
puted, it will be possible to make use of the sum 
of the line spacings only. Adding the line spacings 
of the three perpendicular bands we obtain 
evidently" 


YAv=[(3—E¢)/C—3/A lh/4x? 
=(3/C—7/2A)h/4n°. 


The observations of Bennett and Meyer yield 
LAv = 24.65, 27.15, 28.32 and 28.5 cm™ for CH;F, 
CH;Cl, CH3;Br and CHs3lI, respectively. The mo- 
ment of inertia A has been measured directly from 
the line spacing of the parallel type bands in the 
case of methyl fluoride” and methyl chloride“ and 
found to be 39.5 and 50 X 10°. For the remaining 
methyl halides we have only the estimates based 


2 W. H. Bennett and C. F. Meyer, Phys. Rev. 32, 888 
(1928). 

13 In the preliminary note, Dennison and Johnston, Phys. 
Rev. 47, 93 (1935), an error appeared in this formula which 
slightly affected the resulting values of the moments of 
inertia. Since the publication of the note, a more precise 
measurement has become available for the moment of 
inertia A of methyl chloride which is incorporated in the 
present paper. 

“ME. F. Barker and E. K. Plyler, J. Chem. Phys. 3, 
367 (1935). 

15S. L. Gerhard and D. M. Dennison, Phys. Rev. 43 
197 (1933). 
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upon the envelopes of the parallel bands, namely, 
89 and 99 X10-*°, respectively. These values are 
admittedly not very accurate, but small errors in 
them will have almost no influence upon the 
computation of C since the term 3/C is so much 
larger than 7/2A. 

The indicated substitutions are easily made 
and give the following numbers for the moment of 
inertia C of methyl fluoride to methyl iodide: 
5.61, 5.35, 5.44, and 5.44x10-*°. The corre- 
sponding H —H distances are 1.84, 1.80, 1.81 and 
1.81 <10-° cm. In the next section we shall find 
the moment of inertia of methane A = 5.47 X10~*° 
and H—H=1.81X10-* through an analogous 
method. Direct measurements on the moment of 
inertia of methane through the line spacings of 
methyl deuteride’® are now available and give 
A =5.298 and H—H =1.785. We feel that these 
results are to be interpreted by saying that, within 
the present limits of error, the dimensions of 
the H; triangle and presumably of the CH; group 
remain unchanged throughout the series of me- 
thane and the methyl halides. There are two 
principal sources of error: on the experimental 
side the line spacings of many of the methyl 
halide bands are not very accurately known; on 
the theoretical side we have taken no account of 
the departures from harmonic motion. The total 
error in C could be as high as 10 percent although 
probably it is nearer 5 percent. We do not mean 
to imply that the CH; group undergoes no change 
as we proceed from CH, to CHs3lI, but only that 
our calculations show that whatever change does 
occur, must be smaller, possibly much smaller, 
than 5 to 10 percent. 


§3. THE OVERTONES OF AXIALLY SYMMETRIC 
MOLECULES 


Transitions from the normal state to excited 
states higher than the first result in overtone 
bands. The fine structure spacing of the perpen- 
dicular bands will evidently be given by the 
same formula as that for the fundamentals 
Av=[(1-—¢)/C—1/A Jh/4n°*. It is our purpose to 
compute the ¢’s appropriate for the lower over- 
tones of the molecules YX; and ZYX3. We shall 


begin by treating the overtones 2y;. As usual the 


16 N. Ginsburg and E. F. Barker, J. Chem. Phys. 3, 668 
(1935). 
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wave function in first approximation is separable 
in the normal coordinates. Consider that part 
which depends upon x; and y;. It is now simpler 
to transform to polar coordinates p;, 6; where p; 
is proportional to a displacement say of the Y 
atom perpendicular to the symmetry axis, and 6; 
gives the angle of turning about the symmetry 
axis. Thus p*?7=x;?+y and tan 0;=y;/x;. These 
coordinates and the resulting wave functions are 
identical with those employed in describing the 
CO: molecule.” The wave function of the normal 
state (at least that part depending upon p; and 
6;) is Ry,)®° while the excited state in question is 
characterized by the three functions R,,)%e*?# 
and R,,)*°. The coefficients of 70 show that the 
eigenvalues of (x,yj;—y.i%,) for the excited state 
are +2 and 0. 

Now the electric moment perpendicular to the 
symmetry axis is proportional to x; or to y;, that 
is, to pe*. It is thus clear that the matrix ele- 
ments of the electric moment all vanish since the 
integrand is periodic in 6. This is to be expected, 
for up until now we have considered only the 
lowest order approximation, in which the motion 
is strictly harmonic and in which no overtone 
bands can appear. We must therefore consider 
the modifications to the wave functions caused by 
deviations from harmonicity. The perturbing 
function AV’ may be expanded as a periodic 
series in 30, the factor 3 signifying the threefold 
symmetry of the X; group. Thus 


r J ad = > F; p) medimd | 
m 


For the present discussion where we are inter- 
ested in the possibility but not the actual prob- 
ability of a transition, it is sufficient to consider 
the effect of the perturbation on the wave func- 
tion of the normal state. By the usual formulae 


of perturbation theory we obtain R°%+.A > S™ 


m=—oo 


Xe”? where S” is a function of p. It is now 
evident that no transitions resulting in perpen- 
dicular bands occur to the state R®°. The matrix 
element of the electric moment does however re- 
ceive a contribution from the two wave functions 
Stet? and R*et®® and also from the terms 
S'e—*8? and R”e~*. These considerations give us 
an answer to the question of the appropriate ¢. 


'7 A. Adel and D. M. Dennison, Phys. Rev. 43, 716 
(1933). 


It must have the magnitude 2¢; (from the magni- 
tude of the eigenvalue of (xiyi;—yi%,) and it 
must have the negative sign. This last conclusion, 
namely, that the direction of rotation of the 
electric moment is reversed may be seen from 
the fact that in the fundamental bands we con- 
nect R® and Re’, and the coefficient of 70 
goes 0— 1. In the present case the coefficient of 
10 goes 3—>2 thus giving a reversal of direction. 

The next overtone 3v; corresponds to a transi- 
tion from the normal state to an excited state 
described by the four wave functions R*e*# 
and R**e*#*, By an application of the reasoning 
just employed, it appears that a transition re- 
sulting in a perpendicular type band may occur 
only to R*e**. The ¢ characterizing the fine 
structure spacing is therefore just ¢; itself. In 
this manner we could compute the ¢ for any 
harmonic nv; but the result would be rather 
meaningless for large values of m since here the 
influence of the anharmonic terms will become 
appreciable. The reader must not confuse the 
fact that although we have made use of the 
existence of the anharmonicity in order to find 
the selection rules, we have not taken account of 
its effect on the fine structure. 

One further type of overtone may be discussed 
profitably, namely, the combination band v;+». 
Three possibilities exist. First »; and », may both 
be parallel bands, in which case v;+, is also 
parallel and hence not included in the present 
treatment. Secondly, one of the bands, say vj, 
is perpendicular while the other »% is parallel. 
The presence of energy in a parallel mode of 
vibration does not at all change the internal 
angular momentum and hence the ¢ for the over- 
tone v;+ is just ¢; The third case where »; 
and » are both perpendicular is somewhat more 
complex. The excited state is fourfold degener- 
ate and is specified by the wave functions 
RUR,Me* eit) and RR Me* 1-6), The in- 
ternal angular momentum is additive and de- 
pends upon the coefficients of 6; and 4; hence 
the angular momentum belonging with these 
wave functions is +(¢(;+¢,)h/27 and +(fi—fi)h 
/2m, respectively. The electric moment with 
components yw, and yp, is given by pwetip, 
=a;p.e"*'+a,p.e** where a; and a, are propor- 
tionality factors. We must now consider the in- 
fluence of the anharmonicity on the wave func- 
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tions. The perturbing energy AV’ is composed of 
a sum of terms of the type 


AF (pipx etl eit (3m2—my) Ox) 


where m, and m, are integers. The threefold sym- 
metry of the molecule shows itself in the fact that 
the sum of the coefficients of 6; and 6 must be a 
multiple of three. 

As before it will be sufficient to obtain the 
form of the perturbed wave function of the nor- 
mal state. This is evidently 

ROR+XA = 
m 


1, m2 


Sm, m2pi[m iit (3m2—m})6;) 


A simple trial now makes it clear that no transi- 
tions of the perpendicular type may occur from 
the normal state to the wave functions 
R;R,e**-*®) since the integral is always 
periodic either in 6; or 6. Transitions may occur, 
however, to the states R;R,e**?+%), The par- 
ticular terms of the perturbed normal state wave 
functions which furnish a nonperiodic integrand 
have the following dependence upon the angles, 
e* (0.420) or e*i8i+%) From these considerations 
we see that the ¢ appropriate for the band v;+», 
has the magnitude ¢;+¢, and that it has the 
negative sign. The reasoning which fixes the sign 
is identical with that employed in the case of the 
overtone 27;. 

An application of this result may be made to 
the band of NHs3 which lies at about 2u. The 
band,'® which is of the perpendicular variety, is 
composed of a set of equispaced lines with a 
spacing constant Av=9.98 cm. The position of 
the band center would be consistent with its 
identification as either »;+v4 or ve+v, where v2 
and y, are the perpendicular fundamentals and 
is the high frequency parallel band. In the first 
case 11+, the ¢ would be equal to ¢; and the 
fine structure should resemble that of the funda- 
mental vy. However, this is not what is observed 
since vy shows a strong central region of absorp- 
tion corresponding to Ay~3 or 4. Under the 
second identification, ve+v»,, the appropriate ¢ 
would be —({2+{s) =1—C/2A. Substituting in 
the formula for the line spacing we obtain 
Av= —h/8r°A. That is to say, the spacing should 
be just half the spacing of a parallel type band 
or 9.94 cm=', in very good agreement with the 
observations. 


18 Stinchcomb and Barker, Phys. Rev. 33, 305 (1929). 
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The only other existing data on the fine struc- 
ture of overtones of axial molecules appear to be 
those by Nielsen and Barker’? on methyl 
chloride. The methyl! chloride molecule possesses 
three fundamental perpendicular bands v, v4, 
and v.. The fine structure of each has been meas- 
ured by Bennett and Meyer™ and found to be 
8.2, 12.0 and 6.95 cm™, respectively. From these 
values, together with the moments of inertia A 
and C already cited, we obtain ¢.=0.100, 
a= — 0.268 and ¢,=0.220. 

Nielsen and Barker observed a perpendicular 
band at about 6100 cm™ which they identify with 
2v.. The appropriate ¢ for this overtone we com- 
puted to be —2¢=—0.200, which yields 
Av=11.3 cm~'. The experimental value of the 
spacing was 12.0 cm“. 

The next perpendicular bands which they re- 
solved lie at about 4360 and 4460 cm~! and were 
labeled ve+yv3 and ve+2y;. These correspond to 
the sum of a perpendicular v2 with a parallel fre- 
quency v3 or 2v; and accordingly should have the 
¢=f2and Av=8.2 cm“. Actually they were found 
to have the spacings 8.4 and 8.6 cm™, respec-. 
tively. 

The final perpendicular region observed by 
Nielsen and Barker lies at 2462 cm and shows a 
spacing of 9.7 cm~'!. They identified this band as 
ve+2v; in which case {=f and Av=6.95 cm™ in 
bad agreement with the measurements. A second 
identification which would yield essentially the 
same position for the band center is »s+¥v.5. The 
appropriate ¢ for such a combination is — (¢i+£¢) 
= +0.048 giving Av=8.74 cm™. 

The agreements which have been obtained be- 
tween the predicted and observed line spacings of 
the perpendicular overtone bands of ammonia 
and methyl chloride we regard as very satis- 
factory. The deviations which are in no case 
larger than 10 percent are no doubt caused by the 
effect of the anharmonic terms whose influence is 
of course greater on the overtone bands than on 
the fundamentals. 


§4. SPHERICALLY SYMMETRIC MOLECULES 

The symmetrical molecule YX, possesses four 
fundamental frequencies 7, ve, v3, and vs which 
are single, doubly, triply and triply degenerate, 


1 A. H. Nielsen and E. F. Barker, Phys. Rev. 46, 970 
(1934). 
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respectively. The first two, v; and v2 involve no 
motion of the Y atom and are consequently in- 
active in the infrared. On the other hand », and », 
are both active and correspond in first approxima- 
tion to isotropic harmonic oscillators in space. 
Each of these latter may possess an internal angu- 
lar momentum {,4/27 and we shall first sketch 
the method of calculating £3; and {, as functions 
of the potential constants. 

The tetrahedral system YX, has nine degrees 
of vibrational freedom and hence must be de- 
scribed by nine coordinates. It seems very con- 
venient and straightforward to specify them in 
the following manner. The positions of the four X 
particles may be fixed relative to their center of 
gravity G by means of the mutual displacements 
91, 92, ** * Gg Of the particles along the edges of the 
tetrahedron defined by the equilibrium con- 
figuration X,4. Let gi and qe, gs and gs, gs and ge 
relate to opposite edges of the tetrahedron. 
Finally, the displacement of Y relative to G may 
be given by three Cartesian coordinates x, y, 2. 
These rectangular axes may be chosen in any 
manner, but it is convenient to let x, y, and z 
be lines perpendicular to the opposite edges of 
the equilibrium tetrahedron X, 1 and 2, 3 and 4, 
and 5 and 6, respectively. This numbering of the 
edges corresponds with the subscripts of the q’s. 

The potential energy of the system V can now 
be written as a general quadratic form in the nine 
coordinates. The number of independent con- 
stants is greatly reduced, however, by reason of 
the tetrahedral symmetry and has been shown to 
be only five.2° The number and the nature of 
these constants may be explained by the follow- 
ing description. The axes x, y and z are all equiv- 
alent and hence x’, y’, and 2? have the same co- 
efficient, say a. Similarly all the g? have the same 
coefficient, b. There are two types of cross prod- 
ucts between the q’s, first, cross products between 
adjacent edges (for example qigs3) and second, 
cross products between opposite edges, as 4192. 
These two classes require the constants c and d, 
respectively. Finally, the tetrahedral symmetry is 
unchanged by reversing the direction of one of 
the Cartesian axes, say the x axis, and simultane- 
ously interchanging the two opposite edges in this 
case g; and gz. Thus there may be a cross product 


2 J. E. Rosenthal, Phys. Rev. 45, 538 (1934). 


term of the sort x(g;—g2) but no terms such as xy. 
Hence 


V=pla(x?+y?+2*) +b(g? +9? +9" +9? +95 +4") 
+2c(9igst+Qigst 9195+ Gige+ 929s +9295 
+9295 +9296 +9395 +9396 +9495 +9495) 
+2d(9192+9391 +9596) +2e[ x(q — 92) 

+9(¢3— 4) +2(95—9e) }}. 


The kinetic energy of the YX, molecule was 
found by Dennison! using coordinates identical 
with those defined here, except for a nonessential 
charge in the orientation of the xyz axes. 


T=} {u(?+y?+2*) + (Sm/8) (G1? +422 +43" 
+g +95? +46?) —2(m/8) (Gidstqigst ids 
+9196 + 929s + 9291 +9295 + 9296 + 9395 +9396 
+GsGs+Gago) +2(m/8) (Gige+Gsdst9s9e) }, 


where p=4mM/(4m+M) and m and M are the 
masses of the X and Y particles, respectively. 
The determinant giving the normal frequencies 
v;=,)/22 may be readily solved and yields, 


md,;=4b+4d+16c, mdrA.=b+d—2c, 
m(A3+A4) = ma/u+2b—2d, 
m?(A3A4) = 2m(ab — ad — 2e”) /u. 


The normal frequencies of the YX, molecule 
have already been computed by J. E. Rosenthal?° 
as a function of five constants, A, B, C, Dand E 
which are connected with our constants by the 
linear relations, A =a, B=(b—d)/4, C=(b+d)/4 
—c/2 and D?=e/2. E=(b+d)/4+¢. 

The internal angular momentum associated 
with the frequencies v3 and vs; may be calculated 
by a method essentially similar to that employed 
in the treatment of the YX; molecule. The first 
minors of the normal determinant will give the 
coordinates x, y, 2, g1 «+: ge as linear functions of 
the normal coordinates. These functions in the 
case of the motions v3 and », are indeterminate to 
the extent of the threefold degeneracy. This de- 
generacy may be removed if we assume that the 
angular momentum lies along some definite axis, 


tT). M. Dennison, Astrophys. J. 62, 84 (1925). 
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say the z axis. This means that the amplitude of 
z is zero while the amplitude of x and y are equal 
but ninety degrees out of phase. The amplitudes 
of all the g's, and hence the actual space displace- 
ments of all particles, are then completely de- 
fined. The angular momentum may be divided 
into the two parts: that caused by the motion of 
the Y particle and G, namely, u(xy— yz) and that 
caused by the motion of the X particles relative 
to G. These two contributions are then added 
and evaluated in terms of the normal coordinates 
exactly as was done in the case of the YX; 
molecule. The calculation is somewhat compli- 
cated but eventually leads to the simple result” 


6:=(6;—3)/(6;+1), 


where 6;=4ye?/m(udrA;—a)*?; 1=3, 4. 

The frequencies v; and vs; belong to the same 
symmetry class and are the only fundamental 
frequencies of this class. Hence according to 
Teller’s theorem (3+, should be independent of 
the force constants. Upon making the addition 
and using the relations between \3 and ), we find 
(:+¢4= 3. This result may also be obtained from 
the formula }{=C/2A for the ZYX; molecule. 
In this latter system as the Z atom is changed 
into an X atom and the molecule assumes the 
tetrahedral form, the frequencies ve, v; and v¢, go 
over into the frequencies v3, ve and vs of YX, re- 
spectively.2* Therefore the fo.+f;+{.=C/2A of 
ZYX; becomes (3+f2+{,;=C/2A=4 for YX4. 
We must now consider-the following question. 
What is the internal angular momentum to be 
associated with the inactive frequency v2 of Y Xs? 
By drawing the displacements of the particles it 
is easy to show that each X particle does indeed 
possess an angular momentum, but that these 
angular momenta are all equal and all directed 
either towards, or directly away from the central 
Y atom. Consequently the total angular momen- 
tum vanishes and hence f2=0. Thus {3+{£,;=}3 in 
agreement with the result obtained by direct 
calculation. 

The formulae for the ¢’s are of great practical 
importance in the interpretation of the spectrum 
of methane. In the first place, they enable one to 


2 The formula for ¢; appears to be slightly different from 
that published by us earlier. The difference is due solely to 
a small change in the notation and to the fact that we there 
employed Dr. Rosenthal’s constants. 

23D. M. Dennison, Rev. Mod. Phys. 3, 304 (1931). 
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calculate the moment of inertia by summing the 
line spacings of the bands »; and ». Thus Avs; 
+ Avy =(2—£3—fs)h/42°A =3h/82°A. Secondly, 
they supply a fifth datum which together with the 
values of the four fundamental frequencies just 
allows us to determine the five potential con- 
stants of the molecule. The measurements of 
Cooley** gave the following values for the line 
spacings. Avy3=9.771 cm and Avyy=5.409 cm. 
Following through the indicated substitutions 
we find A =5.47X10~-*, ¢;=0.0345, ¢,=0.4655. 
The molecular dimensions are H —H = 1.81 K 10° 
and C—H=1.11X10- cm. 

The values of the fundamental frequencies 
v= 2914, ve=1520, v3=3014, and ¢,=1304 cm 
together with either ¢3 or ¢4 (the ¢’s are not inde- 
pendent since {=}) serve to fix the potential 
constants. We find 


a=7.670X10' dynes/em c= —0.018 e=+1.807 
b=1.616 d= —0.289 


or in Dr. Rosenthal’s notation 


A =7.670X10' dynes/cem C=0.341 E=0.313 


B=0.476 D=+1.278. 


These data do not determine the sign of e or 
of D, but the recent paper of Ginsburg and 
Barker on the spectrum of methyl deuteride 
shows conclusively that the positive sign is the 
one to be chosen. The work just cited is extremely 
important from our point of view since the values 
of the fundamental frequencies of CH;D together 
with the line spacing of the parallel bands yield 
independent values for the moment of inertia and 
the potential constants of methane. It turns out 
that A =5.298 K10-*°, which is about 3 percent 
lower than our figure. This agreement is quite as 
good as could be expected from a theory which 
neglects second-order perturbations. The poten- 
tial constants obtained by Ginsburg and Barker 
differ from ours by amounts around 5 percent. 
We regard a divergence of this order of magnitude 
as being inevitable for the following reason. Gins- 
burg and Barker in their calculations were forced 
to use the observed positions of the fundamental 
bands and not the actual normal frequencies 
themselves. The difference between these two is 
of the order of 10 to 50 cm™ and will introduce 


24 J. P. Cooley, Astrophys. J. 62, 73 (1925). 
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errors into the potential constants of around 5 
percent. Thus we do not believe that the con- 
stants derived from the spectrum of methyl 
deuteride are essentially more accurate than those 
quoted here. A second way to make this clear is to 
calculate the frequencies of CH;D by using our 
potential constants and the formulae of Dr. 
Rosenthal.** The differences between the calcu- 
lated and observed frequencies are 18, 39, 22, 17, 
7 and 5 cm™. But these are just of the order of 
magnitude of the shift in the position of a band 
caused by the existence of anharmonic terms in 
the potential. 

The spectrum of methyl deuteride presents 
one further point where we may apply our theory 
and that is the fine structure spacing of the per- 
pendicular bands v2, vy and vg. In the formula 
Av;=[(1-—¢,)/C—1/A ]h/42* the moments of 
inertia A and C are known. The internal angular 
momentum factors ¢; are completely fixed by the 
potential constants of methane and the masses of 
the particles. The actual computation of the ¢’s 
was carried out by means of the method used 
previously, but eventually became so involved 
that recourse was had to a numerical evaluation. 
The only check that we have on the calculation 
is that the sum of the ¢’s turns out to have its 
expected value C/2A=0.370. We obtained, 
f2=0.03, &:= —0.25, and ¢,=0.59 from which the 
line spacings Ave=2.3 Avy=5.2 and Avg= —3.3 
cm~ are readily deduced. 

The observations of Ginsburg and Barker are 
as follows. The frequency v2 lies at 2983.0 cm™ 
and has a strong central region of absorption 
indicating Ave~1 cm. This band, however, is in 
strong resonance interaction with the band 27, 
and consequently its line spacing will be greatly 
perturbed and no agreement with the theoretical 
spacing is to be expected. 

The band », occurs at 1477.1 cm™. Its line 
spacing is somewhat irregular but appears to lie 
between 5 and 6 cm™, in good agreement with 
our estimate of 5.2 cm. The frequency vg, at 
1156.3 cm consists of regularly spaced lines 
whose constant ranges from 3.8 to 4.5 cm™. Our 
value of 3.3 does not seem to be seriously in error. 


25 In the letter by Dennison and Johnston the frequencies 
of all the heavy methanes are computed and presented in a 
table which we omit in the present paper. 


§5. THE OVERTONES OF SPHERICALLY SYMMETRIC 
MOLECULES 


A treatment of the overtones of the YX, 
molecule may be made along the same lines as 
that used in the discussion of the overtones of 
ZY Xs, but it is somewhat more complicated in 
that we are here concerned with tetrahedral 
rather than threefold symmetry. We shall con- 
sider only two types of overtone, 2v;, i=3, 4, 
and v3+74. 

Dennison and Ingram” in their paper on the 
overtones nv; of methane, calculated the multi- 
plicity of the excited states. They showed that 
since the motion is similar to that of a three- 
dimensional isotropic oscillator, the wave func- 
tions may be expanded in terms of functions 
ViVi"? where m,+2+n3=n the quan- 
tum number of thestate. ¥ isoneof the Hermitian 
orthogonal functions and the rectangular coor- 
dinates x, y, 2 are proportional to the normal co- 
ordinates. (Because of the degeneracy the mo- 
tions corresponding to v3; and », each require 
three normal coordinates.) The function 
Vio VW via)™ is conveniently abbreviated by 
the symbol (9n3) and it was shown that a 
perturbing function having tetrahedral symmetry 
will connect (M3) with (m'ne'n3’) providing 
ny+mn,', No+ne’ and n3+n;3' are either all even 
or all odd integers. One further statement must 
be made before beginning the analysis. A transi- 
tion which is active must connect two wave 
functions for which one of the quantum numbers 
changes by one unit while the other two remain 
unchanged. Thus (,%%3)—>(m,—1, m2, ms) corre- 
sponds to a nonvanishing electric moment in the 
x direction. 

The normal state is characterized in first ap- 
proximation by (000) but the perturbation will 
add terms of the sort \(m,%en3) where , me, and 
nz are either all even or all odd integers. The first 
excited state is given (in first approximation) 
by the three functions (100), (010) and (001). 
Clearly, transitions to it from the normal state 
may occur resulting in the fundamental band »;. 
The internal angular momentum of this excited 
state is (;h/2m and it is easy to show that its 
component, say along the z axis, is (xy —yz){jh/27 


26D. M. Dennison and S. B. Ingram, Phys. Rev. 36, 
1451 (1930). 
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Thus the eigenvalues of (xy—y.z) for this state 
are +1, 0, —1. 

Dennison and Ingram obtained the result 
that the second excited state »=2 is broken up 
into three states. The first and second of these 
are single and doubly degenerate, respectively, 
and both are described by linear combinations of 
the wave functions (200), (020) and (002). No 
transitions may occur to either of these states 
since it is impossible to select a wave function 
(m\non3) entering the perturbed part of the normal 
state wave function which will combine actively 
with the functions (200), (020) or (002). The 
third part of the state n=2 is triply degenerate 
and uses the functions (011), (101) and (110). A 
transition may occur from the normal state to 
this state by means of the particular term (111) 
in the normal state wave function. Since the 
quantum number decreases rather than in- 
creases in going from the normal to the excited 
states, we would expect the ¢ of 2v; to be nega- 
tive. This may be verified by direct calculation. 
The magnitude of the ¢ depends upon the eigen- 
values of (xy—yaz) for the states (011), (101), 
and (110). A simple computation shows these 
to be +1, 0, and —1. Thus we obtain the result 
that the fine structure of the overtone 2v;, i1=3, 4 
has a spacing constant (1—{)h/4r?A_ where 
f=—%y. 

The excited state leading to the overtone 
nv3+myv, may be characterized by the symbol 
(n\n2n3)[mym2m3] where m+n2+n3=n and m 
+m3+m3;=m. The subscripts 1, 2, 3 refer to the 
xyz directions, respectively. The excited state of 
v3+v, is ninefold degenerate; the wave functions 
consisting of functions of the type (100)[100], 
(010)[001 ], etc. The first problem is to find how 
the tetrahedral symmetry will regroup the nine 
functions. This is easily done by writing out the 
secular determinant, using any perturbing func- 
tion AV’ of tetrahedral symmetry. It soon ap- 
pears that the only nonvanishing elements are 
those for which the changes in m+, m2+m2 
and n3+™m3; are either all even or all odd integers. 
The secular determinant now shows that the nine 
degenerate levels arrange themselves into two 
triple states, one double, and one single state. 
Tetrahedral symmetry will not further reduce 
the degeneracy. 

The first triple state is described by the wave 
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functions 


1/2 
1/2 


1{(100)[001 ]+(001)[ 100}, 
!{(010)[001 ]+(001)[010}} 


and 
1/2${(100)[010 ]+ (010)[100 ]}. 


The second triple state is given by three func- 
tions of the type 1/2!{(100)[001 ]—(001)[ 100 }} 
while the double and single states are given by 
certain linear combinations of the functions 
(100)[100], (010)[010] and (001)[001 }. To find 
the selection rules we must consider the perturbed 
part of normal state wave functions, the most 
significant terms of which for our purpose are of 


the type 
A{(100)[011 ]+(010)[101 ]+(001)[110}}. 


These show that transitions may take place to 
the first triple state but not to the second triple 
state nor to the double or single state. Moreover 
in the term giving rise to the allowed transition, 
the quantum number decreases rather than in- 
creases, showing that the ¢ will have the negative 
sign. 

To obtain the magnitude of the ¢ it is only 
necessary to calculate the eigenvalues of the 
angular momentum along some particular axis, 
say the z axis. The displacements used in com- 
puting the momentum are linear combinations of 
the normal coordinates with proportionality 
factors depending upon ¢; and ¢. Thus the angu- 
lar momentum of the system about the z axis in 
units of h/2rz is, 


(fsbxst fabxa) (Cs'yst Cabs) 
—(fsbvst fabys)(Cahest Fahey). 


The eigenvalues of this function for the triple 
state in question are easily found to be +3({3+£4) 
and 0. Thus the ¢ appropriate for the band 
v3t+vs is —3(f3+¢s) = —} and the line spacing of 
this overtone should be 54/167°A. 

These results may now be compared with the 
experimental data. Although Cooley observed 
the positions of a number of the overtone bands of 
methane, he succeeded in resolving the fine struc- 
ture of only one. He found a faint series of lines in 
the neighborhood of 3.5 with a line spacing of 
about 15.3 cm~. This series he assumed to be 
the negative branch of a band with its center at 
2824 cm and hence he identified it as ve+ 75. 
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Now the state ve has no angular momentum and 
therefore v2+ ; should have the same line spacing 
as vy, that is, 5.4 cm™, in direct disagreement with 
the observations. A plausible way out of this 
difficulty is to assume that the faint series of 
lines found by Cooley is really the positive 
branch of the overtone 275 with its center lying 
at about 2608 cm~'. Such an assumption is 
quite consistent with the observed data as far 
as it is now known. In that case the line spacing 
should be (1+ ¢,)h 427°A =14.8 cm™, which is in 
very good agreement with the experimental 
determination. 

The overtone bands of methane have also been 
studied by Moorhead.” He was able to resolve 
the fine structure of two bands; a band at 1.66 
which gave Av=10.4 cm™ and a band at 2.20u 
with a spacing of 10.7 cm™. The first of these 
Moorhead identified as 23, which according to 
our theory should have a spacing of (1+¢3)h/ 
47°A =10.5 cm™. The band at 2.20u which was 
also located by Cooley is probably v2+ 3. Its 
line spacing should be the same as that of the 
fundamental »; itself or about 9.8 cm™. This 
agreement, while less good, is still satisfactory. 

In a recent paper Steward and Nielsen* report 
on the infrared absorption spectrum of silane, 
SiH. The general features of the spectrum are 
very similar to those of methane with the excep- 
tion of the fact that the fine structure lines are 


27 G. J. Moorhead, Phys. Rev. 39, 83 (1932). 
28 W. B. Steward and H. H. Nielsen, Phys. Rev. 47, 828 
(1935). 


not single but show satellites. This is probably 
due to the influence of the higher order perturb- 
ing terms which have not been included in the 
present theory. Ignoring as far as possible the 
multiplet structure, they obtain for the line 
spacing of the fundamentals »; and », the values 
5.7 and 3.5 cm™, respectively. From these data 
it follows that ¢;=0.07 and ¢,=0.43. Steward 
and Nielsen also observed the positions of a 
number of the overtones of silane and succeeded 
in partially resolving the fine structure of two of 
them, namely, the bands at 3095 cm™ and 3153 
cm~!. These showed spacings of about 9.7 and 
4.7 cm™, respectively. The band at 3095 cm™ 
they identified as v3+»;. According to our theory 
it should have a spacing of 5h/162°A =7.7 cm™. 
The second band they consider to be »+ 7, in 
which case it would have the same spacing as 14, 
that is, 3.5 cm. The agreements in neither case 
are very satisfactory and certainly not as good as 
those found for methane. This is perhaps just 
what is to be expected since, if the higher order 
terms are sufficiently large to produce a splitting 
of the fine structure lines, they will no doubt ap- 
preciably affect the values of the line spacing. 

Note added in proof: Very recently A. H. Nielsen and 
H. H. Nielsen have succeeded in resolving the two overtone 
bands of methane »,+»4 and 3+ (Phys. Rev., this issue). 
They obtain the spacing constants 5.3 and 13.5 cm™, 
respectively. According to our theory the spacing of » +, 
should be just that of », namely, 5.4 cm™ while the spac- 
ing of vs+, is 54/16x°A =12.6 cm™. We wish to express 
our appreciation to the authors of this paper for allowing 
us to see and use their unpublished data. 
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The paper describes experiments made with Geiger- 
Miiller counters in which the cylinder potential is positive 
with respect to the wire instead of negative as ordinarily 
used. The counter operates under reduced efficiency with 
this reversed condition. When used as a double coincidence 
set the reduction in efficiency is greater than for a single 
counter. The paper also describes experiments confirming 
a view that the inefficiency of the reversed counter may be 


explained by supposing the counter action to be produced 
when a positive ion strikes the wire and liberates from it a 
secondary electron. Whenever such an electron is liberated 
a discharge takes place, but the probability of liberation is 
low. However, it seems reasonable that the probabilities 
of such liberation of secondary electrons as found by 
Penning and others stands in agreement with the efficiency 
observed in the counter experiments cited above. 
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i lipo Geiger-Miiller counter connected nor- 
mally with the wire at a positive potential 
with respect to the cylinder, has been used as 
a means of detecting various ionizing rays. 
It is well known that the counting action takes 
place with the potentials reversed, but it has 
been noticed that the rate of a reversed counter 
is much lower than with the potentials con- 
nected in the normal direction. For example, a 
double-coincidence! circuit, in which the record 
was made when the two counters discharged 
simultaneously, recorded two rays per minute 
with the potentials in the normal direction, 
while very few were recorded in a period of 
twelve hours with the potentials reversed. At 
first sight there is nothing remarkable in this 
result since one may contend that the counters 
operate by the release of an electron from the 
gas by an ionizing ray, with the subsequent 
transportation of the electron to the central 
wire, where under the large fields existing there 
it precipitates ionization by collision. With the 
potentials reversed, the released electron is 
carried away from the wire and the conditions 
for ionization do not develop. However, such a 
view discards the possibility that a positive ion 
when carried to the wire under a high potential 
difference may precipitate a discharge by second- 
ary actions involving the emission of electrons 
from the wire. Such possibilities have in fact 
been studied by W. E. Danforth and others. 


POSSIBLE INTERPRETATION OF THE RESULTS 


Two main possibilities suggest themselves as 
an interpretation of the foregoing facts. The 
first instance, the failure of the counters to 
operate when used as a double-coincidence set, is 
to be sought in some cause which results in a 
time delay between the instant of breakdown 
of the two counters, this time delay being greater 
than the recovery time of the circuit used for 
recording. Under such conditions the recording 
circuit would fail to record rays even though 
they had passed through both counter tubes. 
Single counter discharges would, however, be 
recorded by a circuit designed to record singles. 


1 The particular circuit used was that designed by G. L. 
Locher. The counters were argon filled, with their cylinders 
coated with copper oxide. 
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The time delay would presumably arise from the 
fact that the relatively low mobility of the 
positive ions would insure that the liberation of 
such ions from the gas by the cosmic ray at 
different distances from the wire, would result 
in an appreciable difference between the instant 
of arrival at the corresponding wires. 

A possible alternative interpretation would be 
that the counters are inefficient with the po- 
tentials reversed and do not respond to all 
radiation passing through them. 

In order to test this matter, two counters were 
arranged with independent counting circuits, so 
that the recovery time of the circuits was not 
involved, and coincidences were judged by the 
simultaneity of photographs of individual counts 
taken upon the same recording film. Two runs 
were made. In each the counters were 3.8 
centimeters apart, from center to center. In the 
first run, the applied potential of the counter 
reversed was just below the ‘plateau’ of the 
potential-vs.-counting-rate curve of the normal 
counter. In the second run, it was well along the 
“plateau.”’ The first experiment from a single 
counter gave 255 counts in an eight-hour run 
with the potentials reversed. The corresponding 
number of counts with the potentials in the 
normal direction was 2880. The value of the 
relative efficiency ¢, of the reversed counter with 
respect to the normal counter is 255/2880 
=0.088. Hence &, the relative efficiency for 
double coincidences is 7.7 X10-*. The number of 
coincidences observed during the eight-hour run 
was 10, with the potentials reversed. The number 
of coincidences with the potentials applied in 
the normal direction was 960. This ratio, 10~¢ 
to unity, agrees with the value 7.7107 for é 
well within the statistical limits of accuracy 
determined by the 10 counts cited above. 

In the second run the counting rates were 
different, as they should be, but the story of 
their ratios is as follows: With the potentials 
reversed, a run of 200 minutes gave 365 counts 
from a single counter. The corresponding number 
of counts for the same period with the potentials 
normal was 1600. The value of ¢ for this case 
was therefore 365/1600=0.23. The value of & is 
therefore 0.053. The number of reversed co- 
incidences for this period was 21, while the 
normal number with potentials applied normally 
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was 400. The ratio is 0.052, which again agrees 
with the value of e. 

In the foregoing experiments, the simple 
counter records gamma-rays as well as cosmic-ray 
particles, and to no negligible extent, while the 
coincidences between counters are concerned 
only with cosmic-ray particles. It is important 
to observe, however, that the foregoing con- 
clusions are not vitiated by participation of 
the gamma-rays. For while the two sets of data 
from the single counter are influenced by gamma- 
rays as well as by cosmic-ray particles, they 
give the true value of e, while the two experi- 
ments concerned with coincidences, though in- 
fluenced only by cosmic-ray particles, give the 
true value of e&. All that is assumed is that the 
counter possesses 100 percent efficiency for po- 
tential applied normally, and an electron such as 
would be emitted by a gamma-ray from the 
counter wall has the same characteristics for its 
subsequent activities as has that charged particle 
which is the representative of the cosmic ray. 

It is interesting to speculate as to the cause of 
the relative inefficiency of the reversed counter. 
Two possibilities suggest themselves. The first 
invoking only the electrons as the liberated 
entities responsible for the counter discharge. 
Since the electrons move away from the wire 
under the influence of the field, only those which 
are liberated very near the wire where the electric 
field is high can cause a discharge. In the case 
of a single counter the effect is therefore to limit 
the effective areas of the counter to a small 
region around the central wire, and so reduce the 
number of counts recorded. The double-coinci- 
dence counts, depending as they do upon the 
square of the area, will be much more affected 
by reversal than will the single counter; but 
even the single counter will be affected in pro- 
portion to the reduction of the effective area on 
reversal. As a matter of fact, the factor by which 
the single counting rate is reduced by reversal 
should be equal to the ratio of the effective 
counting area to the total counter area. We have 
called this factor « and should expect that the 


coincidence counting rate would be reduced in 
the ratio of & to unity. 

Considering the experimental values of «, the 
limiting fields where an electron would just be 
able to produce a discharge would occur at the 
distance of 0.044 centimeter from the central 
wire for the lower voltage and in the case of 
the higher voltage at 0.12 centimeter from the 
central wire. The respective fields at these points 
are 4905 and 2398 volts per centimeter. Since it 
would be expected on this hypothesis that the 
fields at the edge of the sensitive volume should 
be identical, the discrepancy of the above values 
militates against the adoption of this hypothesis. 
The discrepancy is particularly striking if one 
compares these fields with the fields required for 
counter operation in the normal manner. At the 
threshold voltage the field at the surface of the 
wire is 54,000 volts per centimeter. 

Following suggestions by T. H. Johnson and 
W. E. Danforth, a more probable interpretation 
of the inefficiency of the reversed counter can be 
made by supposing the counter action to be 
produced when a positive ion strikes the wire 
and liberates from it a secondary electron. 
Whenever an electron is so liberated a discharge 
results, but the probability of liberation is low. 
Penning? and others have investigated the proba- 
bility of emission of secondary electrons from 
tungsten by impact of positive ions of several 
gases. The number of electrons per ion varies 
from 0.03 to 0.05 for energies of incident ions 
below 50 volts and increases nearly linearly 
with energy up to 0.42 for 1000-volt ions. 

With the variations of the phenomenon at- 
tributable to adsorb gaseous layers, it seems 
probable that this range of probabilities stands 
in agreement with the efficiency found in the 
above counter experiments. 

The writer wishes to express his appreciation 
of the help given by Professor W. F. G. Swann 
in formulating the final method of attack adopted 
in the problem. 


2? Penning, K. Akad. Wet. Amst. Proc. 33, 841 (1930). 
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Since the capillary ion source described by Tuve, Dahl 
and Van Atta was a decided advance over previous types, 
it seemed desirable to investigate its possibilities further. 
The essential modification in the present design lies in the 
elimination of the auxiliary electrode for drawing ions out 
of the arc. This simplifies construction and greatly reduces 
the difficulties usually experienced in focusing ion beams. 
Currents as high as 4.2 milliamperes have been converged 
into a beam 12 mm wide with a focusing voltage of 10 kv. 
Possible explanations are presented for the observed 
dependence of the total ion current on arc current density, 
condition of the walls, and design of the outlet. Relative 


INTRODUCTION 


ECENT developments in nuclear physics 

have served to emphasize the demand for 
high current sources of protons, deuterons, and 
helium ions. Further tests on the type of pure 
proton source described by Lamar and Luhr! 
have confirmed the conclusion reached by them 
that extreme cleanness of the various metal 
parts is necessary for the production of pure 
proton beams in this type of source. In their 
results there remained also much to be desired in 
the direction of increased ion current for a given 
gas transport. The capillary type ion source 
described by Tuve and his collaborators? gave 
ion currents, for a given gas transport, far in 
excess of those reported for the type of source 
mentioned above, though the proton purity was 
less. The present investigation was undertaken, 
therefore, in an effort to evaluate the various 
factors involved in the operation of the capillary 
discharge and to see in what directions one might 
hope for improvements. This paper is a presenta- 
tion of the results of experiments carried on in 
this laboratory during the past year and reported 
first at the Washington meeting of the Physical 
Society. It is believed that our work, together 
with that of Dr. Tuve and his collaborators, 


1 Lamar and Luhr, Phys. Rev. 44, 947 (1933); Phys. Rev. 
46, 87 (1934). 

? Tuve, Dahl and Van Atta, Phys. Rev. 48, 315 (1935); 
Tuve, Dahl and Hafstad, Phys. Rev. 48, 241 (1935). 


ease of starting and stability of the arc are discussed for 
various designs of tubes. When the arc was operated in 
hydrogen, the mass spectrograph showed a preponderance 
of diatomic and triatomic ions with proton fractions rang- 
ing from 5 to 20 percent of the total ion current. Ions other 
than those of hydrogen were not present in fractions 
exceeding 1 percent of the total. Evidence is presented for 
the breaking up of the heavy ions into protons and neutral 
particles on collision at high speed with gas molecules. 
Upper and lower limits are given for the efficiencies of these 


processes. 


constitutes a reasonably complete picture of the 
capillary type ion source. 


APPARATUS AND EXPERIMENTAL PROCEDURE 


A number of discharge tubes were tested, dif- 
fering in detail one from another, but each having 
the same essential elements. A typical one is 
shown in Fig. 1. The metal envelope which 
enclosed the discharge was made in three pieces; 
a central stainless steel block, and two stainless 
steel tubes press fitted and soft soldered into the 
block. The discharge space consisted of two 
cylindrical enclosures (in every case 3.18 centi- 
meters inside diameter) connected by a smaller 
hole or capillary. There was, in the side of the 
capillary, a tapered outlet through which ions, 
electrons, and neutral gas molecules might escape 
from the discharge. One of the enclosures con- 
tained an anode and the other an oxide-coated 
filament cathode. The stems for cathode and 
anode were supported by ground glass to metal 
joints on the ends of the steel tubes. These joints 
were made vacuum tight by means of picein w 
applied externally. The connection between the 
outlet side of the discharge tube and the mini- 
ature high voltage tube shown in the figure was 
also accomplished by means of picein wax. 
Hydrogen was introduced into the discharge tube 
and was pumped out of the high voltage tube 
with a pumping system having a net speed for 
air of 14 liters per second. 
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Fic. 1. Diagram of experimental tube. 


In the first tube tested the ends of the two 
enclosures, instead of being conical as shown in 
the figure, were flat and the capillary, instead of 
being concentric with the two cylindrical en- 
closures, was on the bottom. In other words, the 
lowest generating line of the capillary was a con- 
tinuation of the lowest generating lines of the 
two cylindrical enclosures. This design was 
adopted in order to minimize the distance 
between the outlet and the first high voltage 
electrode thus giving the most favorable arrange- 
ment for the initial focusing of the ion beam. 
Since no difficulty was experienced in initial 
focusing, and since the lack of symmetry in the 
position of the capillary made the starting of the 
discharge difficult, that design was abandoned in 
favor of the one shown. 

In order to start the discharge the metal walls 
are connected to the anode through a resistance 
of from 20 to 30 thousand ohms. A low current 
discharge then starts in the cathode enclosure, 
the walls of the enclosure acting as anode. 
Electrons from this discharge are shot through 
the capillary and precipitate a high current dis- 
charge to the anode through the capillary. A 
symmetrical field about the end of the capillary 
is most favorable for shooting electrons through 
the capillary and thus for starting the main 
discharge. The symmetrical field also tends 
toward greater stability of the arc for, if any 
transients tend to extinguish the arc, it immedi- 
ately starts again. As regards convenience in tube 
design, it seems safe to say therefore that wide 
variations are permissible provided the walls at 
the ends of the capillary are geometrically sym- 


metrical and, of course, provided cathode and 
anode are not too far removed from the ends of 
the capillary. . 

Various different filaments were tried but the 
one finally adopted consisted of a 20-mil molyb- 
denum wire about 6 inches long wrapped with 
8-mil nickel. The resulting ‘piano wire’’ was 
wound on a 3/32-inch mandrel to form a helix. 
The emitting surface was a mixture of barium 
and strontium oxides. The molybdenum served 
as a heating element and the nickel as a binder 
for the oxides. 

Several 1/4-inch holes (not shown in the 
figure) were drilled through the stainless steel 
block for water cooling and water jackets were 
provided for each of the two stainless steel tubes. 

A number of different capillaries were tried 
ranging in diameter from 1.5 to 4.0 mm and in 
length from 4 to 20 mm. The difficulty in starting 
a discharge increased with increasing length and 
with decreasing diameter of the capillary. In 
some instances the discharge started as soon as 
the e.m.f. was applied and in others it was 
necessary to apply a “leak tester” to the glass 
on the anode side. In most of the tests an e.m.f. 
of 230 volts was used simply for convenience 
although arcs have been maintained easily with 
an e.m.f. of 115 volts. The arc drop ranged from 
60 to 140 volts. In general, the arc drop is least 
for the capillary having the greatest diameter and 
smallest length. For any one discharge tube the 
arc drop is a minimum at some one pressure and 
rises gradually as the pressure is changed on 
either side of the minimum. 

Measurements of the number of ions leaving 
the discharge per second through the outlet were 
complicated by space charge near the outlet and 
by the emission of secondary electrons from the 
electrodes in the high voltage tube. Ten kilo- 
volts, applied to the first electrode, were usually 
sufficient to overcome the space charge near the 
outlet and thus give saturation ion current. With 
this voltage most of the ions went through the 
first electrode and were collected by the second. 
The resulting secondary electrons, emitted by the 
second electrode, were prevented from con- 
tributing to the measured ion current by means 
of a small retarding potential (45 to 350 volts) 
maintained between the two high voltage elec- 
trodes. 
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The ion current was measured for each dis- 
charge tube as a function of the variables in the 
arc proper. These were the arc current, the gas 
pressure, the potential of the steel walls with 
respect to the anode, and the previous history of 
the tube. A slit (5 mm long and 1.5 mm wide) in 
the bottom of the second high voltage electrode 
allowed a small fraction of the ions to enter a 
mass spectrograph (not shown in the figure). It 
was possible, therefore, to determine the nature 
of the ions in a sample of the beam as a function 
of the variables mentioned above. 


RESULTS AND CONCLUSIONS 


The current of ions leaving the discharge 
through the outlet depends upon the number of 
ions produced per second in the discharge and 
upon the electric field in the vicinity of the outlet. 
These two variables and the variables on which 
they in turn depend will now be considered. 

The number of ions produced per second 
depends upon the current density and potential 
drop in the capillary and upon the concentration 
and nature of the gas present. Experimentally it 
was found that, with the metal walls floating or 
connected to the anode through a high resistance, 
the ion current (leaving the discharge through 
the outlet) for constant arc current varied over 
fairly wide limits, in most cases by a factor of 
two. When an arc was first started, the ion 
current was always low and increased with con- 
tinued operation. Since there was no accompany- 
ing change in the relative numbers of ions of 
different masses collected in the mass spectro- 
graph or of are drop, it was concluded that the 
increase in current resulted from an increase in 
the current density in the capillary rather than 
from any change in the nature of the gas or drop 
along the capillary. The following is a possible 
explanation of the change. There are two possible 
paths through which current can flow from anode 
to cathode: the capillary, and the metal walls. 
The fraction of the arc current which follows the 
metal path (and thus produces no ions) is greater 
or less, depending upon the number of secondary 
electrons emitted by the walls of the anode en- 
closure. Continued operation of the discharge 
tends to remove surface films from the walls, thus 
reducing the number of secondaries emitted and 





COMPTON 


AND 


forcing the arc current to take the path through 
the capillary. If the walls are made negative with 
respect to the cathode, electrons from the cathode 
are prevented from reaching the walls and thus 
most of the arc current must go through the 
capillary. This procedure does not result in any 
great increase in ion current for a given power 
input since the current resulting from the emis- 
sion of secondaries from the walls of the anode 
enclosure must be supplied by the metal wall 
circuit. Reductions in the fraction of the are 
current taking the metal path can be accom- 
plished by selecting, for tube construction, a 
metal which does not emit many secondaries and 
by reducing to a minimum the area of metal 
exposed to the discharge on the anode side of the 
capillary. 

Let us consider now the electric field in the 
vicinity of the outlet and its effect upon the ion 
current. In general a metal surface in contact 
with an ionized gas and negative with respect to 
it is covered with a positive ion space-charge 
sheath. The thickness of the sheath increases 
with increasing negative voltage and with de- 
creasing ion density. In the capillary, then, the 
region of ionized gas in the center is surrounded 
by a space charge sheath separating this region 
from the walls. The thickness of the sheath 
depends upon the ion density in the capillary and 
upon the potential of the central region with 
respect to the walls. Consider now ions from the 
discharge which leave the capillary through the 
outlet in the side of the capillary. If the thickness 
of the sheath is large in comparison with the 
diameter of the outlet, the outlet introduces no 
distortion in the sheath. Thus ions drifting into 
the sheath from immediately over the outlet 
suffer accelerations normal to the sheath and 
leave the discharge through the outlet. As the 
sheath thickness decreases, a depression develops 
immediately over the outlet. Ions which then 
drift into the sheath are accelerated in a diverg- 
ing field and thus not all of them get through the 
outlet. If the sheath is extremely thin it bends 
around the edges of the outlet thus allowing a 
part of the discharge itself to stick out through 
the outlet, the ions being removed from this part 
of the discharge by the field in the high voltage 


3 Langmuir and Mott-Smith, G. E. Rev. 37, 449, 538, 
616, 762, 810 (1924). 
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Fic. 2. Ion current vs. potential of the metal walls with 
respect to the cathode. Arc current 1 ampere, pressure 
0.095 mm, capillary diameter 4 mm, and outlet diameter 
0.83 mm. 











tube. The fraction of the area of the outlet 
available for this “sticking out’’ of the discharge 
increases, therefore, with decreasing sheath 
thickness and thus, in this region, the ion current 
increases with decreasing sheath thickness for 
constant ion density. The correctness of this 
picture is indicated by the results presented in 
Fig. 2. In this figure the ion current is plotted as 
a function of the negative potential of the metal 
walls with the cathode as a fixed reference, other 
measurable variables remaining constant. As one 
would expect from the above picture, for high 
negative voltages, as the negative voltage and 
thus the sheath thickness increase, the ion 
current increases. In the intermediate region the 
ion current is low. As the sheath thickness 
decreases after having passed the minimum, the 
ion current again increases indicating a ‘‘sticking 
out” of part of the discharge through the outlet. 
It was not possible to extend this curve in either 
direction in search of a leveling off, since further 
increase in negative voltage put the arc out, and 
since any decrease made the walls positive with 
respect to the cathode, thus allowing electrons to 
reach the walls and complicating the interpre- 
tation of the curve. Maximum ion currents are 
obtained, therefore, either from a number of 
extremely small outlets or from one large one. 
Fig. 3 is a plot of ion current leaving the dis- 
charge with arc current as the independent vari- 
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Fic. 3. Ion current vs. arc current with metal walls con- 


nected to anode through 20,000 ohms. Pressure 0.15 mm, 
capillary diameter 3.2 mm, and outlet diameter 1.2 mm. 











able. If the effective area of outlet were constant, 
one might expect an approximate proportionality 
between ion current and arc current. But since 
the sheath thickness decreases with increasing 
ion density and thus with increasing arc current, 
the effective area of outlet is not constant. Since, 
in this case, the outlet was relatively large (1.18 
mm in diameter) it is believed that the sheath 
was small by comparison and thus the effect of 
changing sheath thickness on outlet area was 
small, and thus the lack of any great curvature 
is taken as good evidence for a proportional 
dependence of ion current on arc current for 
constant outlet area. 

A number of tests were made to determine the 
dependence of ion current and arc stability on 
gas pressure. Over the range from 0.035 to 1.0 
mm no marked dependence of ion current on gas 
pressure was observed although it increased 
slightly with increasing pressure. In pure hy- 
drogen, however, for pressures below 0.06 mm 
the arc shows a marked tendency to go out. A 
pressure of 0.1 mm is suggested, therefore, as a 
sufficiently high pressure to insure reliable opera- 
tion. As far as the operation of the arc is con- 
cerned, it can be made stable at low pressures by 
introducing traces of impurities into the hy- 
drogen although this is undesirable since under 
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these conditions an appreciable fraction of the 
ions are those of the impurity. 

The design of capillary arc shown in Fig. 1 
presents a different focusing problem from that 
usually encountered since there is no auxiliary 
electrode for drawing the ions out of the arc. In 
this problem there are two limiting cases: one in 
which the outlet diameter is large in comparison 
with the sheath thickness so that a part of the 
discharge actually sticks out through the outlet, 
and the other in which the outlet is small so that 
the ions are shot out by the potential drop in the 
sheath. In either case there is very little net 
space charge in the immediate vicinity of the 
outlet, for in the first case electrons drifting out 
of the discharge with the ions tend to neutralize 
space charge, and in the second not only does 
the smallness of the outlet reduce the number of 
ions but the sheath drop furnishes them with an 
initial velocity. Studies of focusing were greatly 
simplified by the fact that slight traces of gas in 
the high voltage tube rendered the beam suf- 
ficiently luminous to be clearly visible in a 
darkened room. With low voltage on the first 
high voltage electrode, the beam appeared to 
come from a rather large area indicating con- 
siderable spreading of the ions after leaving the 
outlet. As the voltage was increased, not only did 
this area decrease, but the ions continued to 
converge in their flight down the tube. With an 
ion current of 4.2 ma and a voltage of 10 kv the 
beam reached a minimum width of 12 mm at a 
distance of about one centimeter below the first 
high voltage electrode. Since this beam width 
was sufficiently small to allow the ions to pass 
through the first high voltage electrode, and 
since the work of Dr. Tuve and his collaborators? 
on further focusing was already so complete, 
these studies were not continued. It should be 
added, however, that the apparent beam edge 
was always extremely sharp. That this corre- 
sponded with the physical facts was shown by the 
sharp drop observed in the current collected by 
the first electrode when the beam appeared just 
to clear its edges. There are, therefore, certain 
advantages to be gained by allowing the ions to 
leave the discharge with negligible velocities and 
impart all of the energy in a region of rela- 
tively high vacuum. Under these conditions, one 
is not faced with the problem of refocusing an 
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and direc- 
which gas 


initial beam, heterogeneous in energy 
tion, which comes from a region in 
scattering is possible. 

Studies of totai ion current were, 
complicated by the passage of part of the arc 
current through the metal and by the dependence 
of the effective area of the outlet on sheath 
thickness, as discussed above. The indications 
were, however, that the ion density, and thus the 
ion current for constant effective area of outlet, 
were approximately proportional to the current 
density in the capillary. With a 3.2-mm capillary 
and a 1.2-mm outlet, the ion current for a 1.92 
ampere arc was 4.2 ma. Thus one may expect 
about 16 ma per square centimeter effective area 
of outlet for a current density of 1.0 ampere per 
square centimeter in the capillary. In the 3.2-mm 
capillary, the current density for a 1.0-ampere 
arc is about 12.5 amperes per square centimeter. 
Currents as high as 5 amperes have been passed 
through such a capillary, giving a current density 
of 62.5 amperes per square centimeter and thus 
an ion current density of 1000 ma per square 
centimeter. For the sake of comparison, the 
source described in the first reference gave 8 ma 
per square centimeter for a 1.0-ampere arc. 

In Fig. 4 a fairly typical mass spectrograph 
curve is presented merely to show the positions 
and shapes of the various peaks. Although a 
number of complete curves were taken, in the 
interest of speed only the heights of the various 
peaks were recorded in most of the tests. In most 
of the tests taken with the mass spectrograph, 
the two high voltage electrodes and the mass 
spectrograph were 1160 volts negative with 
respect to the steel walls of the discharge tube. 
As can be seen, the diatomic and triatomic ions 
are in the majority, the relative magnitudes of 
the two peaks depending on the pressure, while 
the primary proton peaks are disappointingly 
low. With pure hydrogen, the fraction of the ion 
current carried by primary protons ranged from 
5 to 20 percent, and was usually about 10 
percent. No marked dependence, of primary 
proton percentage on the arc variables, was ob- 
served over the ranges of currents and pressures 
used although there was some slight indication 
of an increase with measured pressure. There was 
some evidence that the addition of helium to the 
hydrogen tended to increase the proton per- 
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F1G. 4. Typical mass spectrograph curve. Pressure 0.13 mm. 


centage, but the evidence was not conclusive 
and since the effect, if present at all, was small, 
only a few tests were made. 

It was hoped that the mechanism employed by 
Lamar and Luhr' for producing protons might 
be of use in this type of ion source. Briefly, if both 
ions and neutral particles have a low accom- 
modation coefficient, and if the walls are main- 
tained highly negative, ions striking the walls 
rebound as neutrals retaining a large fraction of 
the energy gained in the negative potential drop. 
This energy going into the gas tends to increase 
its temperature. If the neutral particles also have 
a low accommodation coefficient large tempera- 
ture differences are possible between gas and 
walls. With the gas at a high temperature, the 
hydrogen exists in the atomic rather than the 
molecular state and the ions are therefore 
protons. Tests on the capillary type source and 
further tests on the type described in the refer- 
ence above confirmed the conclusion that ex- 
tremely clean surfaces are required for the 
existence of low accommodation coefficients. 
This surface cleanness, although possible, re- 
quires care in technique which is hardly practical 
in a high voltage installation. 

Some hope for increased proton percentages is 
offered by the fact that heavy ions at high speeds 
break up, on collision with gas molecules, into 
protons and neutral particles. That this occurs 
is indicated by the presence of peaks in Fig. 4 
at 0.4 and 0.5 ampere magnet current. If these 


are protons, they appear with one-half and one- 
third energy as would be expected if they were 
formed by the break-up of H2*+ and H;"*, respec- 
tively. This conclusion is further substantiated 
by the fact that the relative heights of these two 
peaks vary in the same way as the relative 
heights of the H2*+ and H;* peaks. The peak at 
0.8 ampere results from the break-up of H;* ions 
to form Hg. It is possible, on knowing the relative 
heights of the various peaks, the pressure in the 
discharge tube, and the various dimensions of the 
apparatus, to get some idea of the efficiency of 
these break-up mechanisms, the only uncertainty 
being the angular distribution of the protons 
after break-up. A lower limit to the efficiency of 
break-up can be set by assuming no angular 
scattering. With this assumption, every heavy 
ion, which would have entered the spectrograph 
had it not been broken up, produces a proton 
which does enter the spectrograph. Knowing, 
then, the number of protons per heavy ion, the 
length of path in the discharge tube at constant 
potential, and the pressure in the discharge tube, 
one can calculate from these data the average 
number of collisions per centimeter path at 
normal pressure and temperature which result in 
break-up. This calculation gives, as a lower limit, 
for the H,* five and for the H;* two collisions 
per centimeter path at one millimeter pressure 
and room temperature resulting in break-up. As 
an extreme upper limit, one might assume that 
on break-up the protons are scattered equally in 
all directions. Knowing the solid angle subtended 
by the slit system of the spectrograph at any 
point along the path of the ion beam, one is able 
to perform the necessary integrations and solve 
for the efficiency necessary to account for the 
observed number of protons. This assumption 
leads to about 250 times as many collisions per 
centimeter path as the former one. A more 
reasonable assumption is that the protons are 
scattered as though they had collided as elastic 
spheres with spheres of equal mass. Since the 
resulting distribution gives four times as many 
protons scattered in the forward direction as the 
equal angular distribution, the ratio between 
upper and lower limits is now reduced to about 
60. For 1160-volt ions, therefore, the number of 
collisions per centimeter path at normal pressure 
and temperature resulting in break-up lies some- 
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892 LAMAR, 
where between 5 and 300 for H2* ions and 
between 2 and 120 for H;* ions. These tests were 
not repeated for other voltages. However, the 
absence of marked broadening of the primary 
proton peak for high pressures, and negative 
voltages of the capillary walls as high as 180 
volts, indicated that the efficiency of break-up 
must be extremely small for voltages less than 
180. 

The possibility of using the break-up mechan- 
ism for producing protons depends, of course, 
upon the relative importance of scattering in 
removing ions from the beam as it passes through 
the gas. In order to test this an additional elec- 
trode was introduced 4.0 centimeters below the 
outlet to the tube. It consisted of a cylinder 1.0 
centimeter in diameter and 4.0 centimeters long 
with a plate extending from the upper edge of the 
cylinder radially to the glass walls of the dis- 
charge tube thus sealing this region from the 
region of high vacuum. Previous tests, without 
the additional electrode, having shown little 
change in total ion current with pressure, it was 
assumed that, with the additional electrode any 
decrease in ion current with increasing pressure 
resulted from scattering in the region between 
outlet and additional electrode. The measure- 
ments indicated an approximate exponential 
dependence of ion current on pressure as would 
be expected from loss due to scattering. The loss 
for a gas path equivalent to 1 mm at 1 mm 
pressure was 3/4 of the current in the beam. With 
this gas path the fraction of primary protons 
increased to 30 percent. These were definitely 
primary protons from the discharge since they 
possessed the full energy and their high per- 
centage merely showed preferential scattering of 
the heavier ions. In all of the tests with the addi- 
tional electrode, the maximum current was 
realized with a potential of about 3 kv on the 
electrode, showing for lower voltages an under- 
focusing and for higher an overfocusing of the 
ion beam. 

Since the angular distribution of the protons 
after formation from the break-up of heavier ions 
is not known it is not possible at the present 
stage to give an accurate estimate of the ultimate 
limit to the total proton current. However, a few 
remarks as to the directions in which to proceed 
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in designing a “‘break-up’’ region seem to be in 
order. The design of such a region presents two 
problems: that of breaking up the ions and of 
focusing the resulting protons. The first of these 
problems is relatively simple. One needs, simply, 
to give the ions an energy of 1 kv or more and to 
have sufficient gas present to insure the desired 
number of break-up collisions. In order to focus 
the protons the field in the break-up region 
should be strongly converging so that scattered 
ions may be drawn back into the ion beam. The 
region should be wide in comparison with its 
length so as to provide the field with ample 
opportunity to act on the ions before they reach 
a wall and are lost. In the source described by 
Tuve and his collaborators? some of the lines of 
force from the auxiliary electrode pass through a 
constriction, separating this electrode from the 
discharge, and up into the capillary. This ar- 
rangement provides a converging field in a region 
of reasonable gas pressure, but the voltage drop 
in this leakage field is necessarily low since high 
voltages here extinguish the arc. Since the efh- 
ciency for break-up with low voltages is small, 
few protons are produced in this region. It seems 
more advisable, therefore, to get the ions com- 
pletely out of the discharge, possibly by one of 
the mechanisms outlined above, and then provide 
them with the necessary energy in the region of 
gas. The outlet to the break-up region should be 
a hole in a diaphragm rather than a long tube 
since the angle spread of the beam leaving this 
region is probably large. A discussion of the 
relative merits and demerits of diaphragms as 
compared with tubes is given in the second 
reference. Although the “break-up”? mechanism 
for forming protons has with it the undesirable 
feature of gas scattering which makes focusing 
difficult, it does not make focusing impossible and 
seems to offer the most promising means of 
obtaining high proton yields. 

In conclusion, the authors wish to express 
their indebtedness to Mr. Walter Kallenbach, 
head of the Physics Shop, whose advice in design 
and whose cooperation in execution made 
possible the testing of a large number of tubes. 
We wish also to thank Dr. C. M. Van Atta for 
his excellent suggestions and advice during the 
preparation of the manuscript. 
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Polymorphism, Principally of the Elements, up to 50,000 kg/cm’ 
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A new technique is described by which pressures of 
50,000 kg/cm? and more may be applied to solids, and the 
parameters of any transitions measured. A systematic 
examination has been made for polymorphism of many of 
the elements in the new pressure range. New modifications 


are found for Bi, Hg, Tl, Te, Ga, and Ip, and the transition 
parameters measured. A beginning has been made of a 
systematic study of polymorphism of compounds, and 
results obtained for KCI, KBr, and KI, which assume at 
about 20,000 kg/cm? presumably the CsCl type of structure. 





THE NEw TECHNIQUE 


ITHERTO nearly all my high pressure 

measurements have been restricted to a 
range of about 20,000 kg/cm*, and by far the 
larger part have reached only 12,000. The limit 
has been set by the bursting strength of the 
cylinders. By subjecting the cylinders to an 
initial stretching considerably beyond the elastic 
limit it is possible to put them into such a state 
of internal stress that the best grades of steel 
will withstand pressures of the order of 40,000 
kg/cm?, but the benefit of any such treatment 
rapidly tapers off with time, so that 20,000 is as 
high as I have found it profitable to go, and even 
then rupture is frequent. The practical limit of 
the pistons with which pressure is produced is 
considerably higher, around 40,000 kg/cm’. 

It is evident, and in fact a number of persons 
have independently made the suggestion to me, 
that it would be ideally possible to reach any 
pressure whatever by making a nest of pressure 
apparatuses, one within the other, each inner 
apparatus being exposed to a uniform hydrostatic 
pressure all overexerted by the next outer ap- 
paratus. The extreme mechanical difficulty of 
carrying out this idea, even in an apparatus of 
only two stages, has hitherto deterred me from 
any such attempt at extending the pressure 
range. It is quite conceivable, however, that it 
would be possible to accomplish the desired end 
by subjecting the internal apparatus to a less 
complete support than would be afforded by a 
hydrostatic pressure all over its external surface. 
Such, has in fact, proved to be easy, and in the 
following results are given for a large number of 
routine experiments made with such an apparatus 
up to 50,000 kg/cm’. 

The cylindrical vessel exposed to internal 


hydrostatic pressure may be supported with suf- 
ficient effectiveness by applying an external 
pressure to the curved surface only of the 
cylinder. A very simple calculation by elasticity 
theory shows the effectiveness of such support. 
Consider, for example, two cylinders both ini- 
tially of the same internal diameter, one with 
infinitely thick walls, and the other with an 
external diameter four times the internal diam- 
eter. To the first cylinder is applied a certain 
pressure, and to the second cylinder the same 
pressure and at the same time an external 
pressure equal to one-half the internal pressure. 
Calculation by the conventional formulas, as for 
example in Love’s Elasticity, shows that the 
radial stretch at the inner surface of the second 
cylinder is only 0.28 as great as that of the first. 
A factor so advantageous as this cannot be gained 
by cold working. If the ratio of external to in- 
ternal pressure could be automatically main- 
tained in some way as the internal pressure in- 
creases, the second cylinder could obviously be 
carried to a much higher pressure without rup- 
ture than could the first cylinder with its in- 
finitely thick walls. 

A simple scheme by which the external pres- 
sure may be automatically maintained at a pre- 
determined fraction of the internal pressure is 
illustrated in Fig. 1. The external walls of the 
cylinder are given a conical shape like a stopper, 
and this steel ‘stopper’ is pushed into an 
external sleeve by the same force with which the 
piston is driven into the cylinder, compressing 
the material within it. The friction on the walls 
of the cone is of course a complicating factor 
that must be allowed for empirically, but by 
using carefully polished surfaces, and a graphite 
and vaseline lubricant the purpose is accom- 
plished with about the proportions shown, and 
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Fic. 1. Fic. 2. 


Fic. 1. Illustrates the new arrangement by which an 
external pressure is automatically applied to keep pace 
with the internal pressure. The diagram is approximately 
to scale; the diameter of the ‘‘stopper”’ at the upper end is 


2.5 cm. 
Fic. 2. The double ended cylinder and pistons by which 
the transitions of bismuth were measured. 


the stopper comes back out of the sleeve without 
sticking on release of pressure. 

As already stated, with this arrangement many 
routine measurements have been made to 50,000 
kg/cm?. The weak part of the apparatus is now 
the piston, 50,000 being beyond the reach of any 
steel that I have been able to find. Fortunately 
carboloy, the recently developed cemented 
carbide‘of tungsten and cobalt, turns out to have 
a much higher compressive strength than the 
hardest steel, and 50,000 kg/cm? or even con- 
siderably more may be easily reached with it; I 
am exceedingly indebted to Dr. Zay Jeffries of 
the Carboloy Company of America, through 
whose kind offices I was presented with sufficient 
carboloy in the form of rods about 6 mm in 
diameter to make 100 pistons 6 mm long. It is 
usually possible to use a single piston for half a 
dozen excursions to 50,000 kg/cm? before it 
breaks. 

Full details of the apparatus are deferred to 
another paper, only a few of the more important 
points will be mentioned here. The simple ap- 
paratus sketched in Fig. 1 is evidently capable of 
yielding quantitative information about poly- 
morphic transitions by merely determining the 
position of the piston as a function of pressure, 
which is exactly the method used in my former 
work at lower pressures. It is probable that the 
apparatus can also yield information about com- 
pressibility, although this second use has not 
been seriously attempted yet. There will obvi- 
ously be complications because of hysteresis and 
friction, which, however, I believe can be over- 
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come. In any event, the accuracy cannot be as 
high as in my previous work, in which pressure is 
truly hydrostatic and is transmitted by a true 
liquid. But such experiments are no longer 
possible at ordinary temperatures at pressures 
as high as this, for nearly everything freezes to 
the solid. The substance which would longest 
resist freezing is doubtless helium; Simon's! ex- 
trapolation indicates that even it will freeze at 
room temperature at a pressure in the neighbor- 
hood of 100,000 kg/cm?, a pressure which I have 
actually reached under special conditions on one 
occasion. It has been a suprise, however, to find 
that the apparatus is capable of accuracy not 
much lower than that of the former apparatus, 
as has been established on several occasions 
when measurements of the same transition have 
been made with the two types of apparatus. In 
any event, the field is sufficiently new, and the 
results to be expected of sufficient interest to 
justify the use of methods not of the highest 
accuracy. 

One reason that the pressure fails to be truly 
hydrostatic arises from the friction exerted on 
the walls by the material undergoing compression 
as it moves down into smaller volume. If the 
ratio of length to diameter is high, enormously 
high friction may be developed in this way, so 
that the pressure at the lower end may be only 
a small fraction of the pressure immediately 
under the piston. In this apparatus the dimen- 
sions were such, however, 6 mm diameter and 
length of the order of 9 mm, that the friction 
could not be an important fraction of the mean 
hydrostatic pressure. The friction is less when 
soft materials are used, because of their low 
plastic limit, and it is this kind of material that 
has been used in most of these investigations. 
Harder materials can be experimented with by 
embedding them in a sheath of some softer 
material, such as lead or even copper. Such a 
sheath becomes a positive necessity for many 
materials, particularly above 100°. In my early 
work I had found that mercury amalgamates 
and breaks the cylinder at a few thousand kilo- 
grams pressure. Action quite similar to this takes 
place with many of the low melting metals at 
these new higher pressures. Sodium, for example, 


1F, Simon, M. Ruhemann and W. A. M. Edwards, 
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cannot be used directly in contact with steel at 
any temperature, but alloys with it and ruptures 
it. Tin, bismuth and others act similarly. Also 
such substances as sulfur attack chemically the 
steel, which must be protected from direct 
contact. 

In addition to the friction of the material 
under compression, there is also considerable 
friction at the closure under the piston, which is 
an arrangement utilizing a conical ring as shown. 
The effect of both sorts of friction was eliminated 
from the transition pressure by taking mean 
values obtained with increasing and decreasing 
pressure. One serious result of the very consider- 
able friction, however, was loss of sensitivity; a 
number of substances were found to have 
probable transitions by a shearing method, de- 
scribed in a previous paper, but showed nothing 
large enough to be measurable by this direct 
volumetric method. 

Temperature control was effected through a 
temperature bath, which consisted of a pan 
soldered to the conical supporting sleeve. Solid 
carbon dioxide temperatures were easily reached 
in this way, and by means of an electric heater in 
oil, temperatures up to the softening point of the 
steel were reached. 200°C was about the limit for 
measurements to 50,000; at higher temperatures 
the cylinder softens perceptibly, as does also the 
carboloy piston, which seems to be notably 
stronger at the lower temperatures. 

The detailed results now follow. First are given 
those elements for which it has been possible to 
carry through a full determination of all the 
thermodynamic parameters, namely, transition 
pressures and accompanying change of volume 
as a function of temperature. Next those ele- 
ments are enumerated for which negative results 
were found, with a list repeated for the sake of 
completeness from a previous paper of those 
elements which the shearing method has shown 
to be probably polymorphic. Finally are given 
the transition parameters for the potassium 
halides, with a mention of other exploratory work 
for a few simple compounds. 


DETAILED RESULTS FOR THE ELEMENTS 


Bismuth 


This metal was found to be polymorphic at 
high pressures over a year ago, and part of the 
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results have been already summarized in print.? 
The first work with bismuth was done with a 
different type of apparatus from that described 
above, not capable of nearly such high pressures, 
but of sufficient usefulness to justify a brief 
description. In general idea it is not dissimilar 
from a method used by Tammann.* The method 
is sufficiently well indicated in Fig. 2. The bis- 
muth is contained in a short cylinder, and com- 
pressed from both ends by two pistons pushed 
together by a hydraulic press. By the use of two 
pistons friction is very materially reduced, and 
the sharpness of the readings increased. Mechan- 
ical stability is increased by the large ends shown 
on the pistons. The cylinder is made of ‘‘Solar’’ 
steel; this is a water hardening silicon-manganese 
steel of very materially higher tensile charac- 
teristics than the chrome-vanadium steel used 
hitherto in my work. In fact, it was the avail- 
ability of this new grade of steel that gave the 
impetus to all this attempt to reach higher 
pressures. The solar steel cylinder of Fig. 2 was 
treated by an initial stretching exactly as my 
former cylinders. Temperature was varied by 
means of a heating coil wound over the outside 
of the cylinder. With this arrangement measure- 
ments could be made up to 30,000 kg/cm? or 
perhaps more. 

With this first apparatus two new high 
pressure forms of bismuth were discovered, and 
measurements were made on the transition lines 
I-II and II-III between 79° and 182°C. At 197° 
the pressure of the melting of I was located, 
15,700 kg/cm’, increasing. pressure from the 
domain of the solid to that of the liquid. The 
pressure found in this way agreed well with the 
pressure extrapolated from former measurements 
up to 12,000. The melting pressure had not been 
held for longer than two minutes, however, when 
the cylinder broke, amalgamated by the liquid 
bismuth; the amalgamation was very obvious in 
the fracture. 

With the “‘stopper”’ apparatus an exploration 
was made for other modifications at still higher 
pressures. The I—II and II-III transitions were 
first located at 134° and the change of volume 
determined. At somewhat less than 30,000 at 


2 P. W. Bridgman, Phys. Rev. 47, 427 (1935). 
3G. Tammann, Zeits. f. anorg. u. allgem. Chemie 199, 
209 (1931). 
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this temperature the cylinder broke, amal- 
gamated by the solid bismuth. The changes of 
volume and also the transition pressure so found 
agreed well enough with the values obtained 
with the first apparatus. A new cylinder was now 
used and the bismuth enclosed in a copper sheath. 
This was successfully carried to 50,000 at 20°. 
At the upper end of the pressure range there 
was considerable plastic yield of the cylinder, so 
that a small transition might have escaped 
notice. On release of pressure, a very small dis- 
continuity in the volume of not more than 5 
percent of the volume change of the other transi- 
tions was found near 40,000. No special sig- 
nificance was attached to this at the time. But 
on making shearing measurements with the 
apparatus described in the previous paper, there 
appeared unmistakable evidence that this discon- 
tinuity must have been really a transition. The 
experiment was then repeated with two different 
set-ups of the “‘stopper’’ apparatus at —79° and 
0°C. On each occasion the break was found with 
decreasing pressure and there seems no doubt 
of the reality of the transition. The “region of 
indifference” of this transition unfortunately 
appears to be rather wide, the transition 
occurring with increasing pressure so near to 
50,000 that it is obscured by volume creep in the 
cylinder. It was therefore not possible to deter- 
mine the average of the transition pressure with 
increasing and decreasing pressure, so that the 
true transition pressure is restricted to a much 
broader region of uncertainty than is usually 
possible. The width of the band of indifference 
did not seem to be importantly less at 20° than 
at —79°, and for this reason the attempt was not 
made to get more accurate values by increasing 
the temperature, a procedure which would 
ordinarily be effective, but here would have been 
of doubtful value because of the rapid increase 
of plastic creep with higher temperatures. The 
volume change was so small that it could be 
determined with only small accuracy; the actual 
value found at —79° was a little more than one- 
half that at 20°. 

At solid CO, temperature the two transition 
points I-II and II-III are so close together that 
it was not possible to separate them. The total 
change of volume, I-III, found at this point was 
0.0073, which agreed exactly with the sum of the 
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values for I-II and II-III extrapolated from 
readings at higher temperatures. In view of the 
breaking of the cylinder when melting of I 
took place, no attempt was made to obtain 
points on the melting curves of II or III. The 
parameters at the triple points given in Table I 
were obtained by calculation. The parameters for 
the melting of II and III are very rough, since 
they involve extrapolations and the difference of 
nearly equal quantities. The parameters of the 
transition III—IV given in the table are also very 
rough. 

The transition parameters are collected in 
Table I, and the phase diagram is shown in Fig. 3 
and the changes of volume in Fig. 4. It will be 
seen that there is a rough analogy between the 
phase diagram of bismuth and that of water, 
which possibly may be of significance. 


Mercury 


The transition was first found at CO. tem- 
perature at about 14,000 kg/cm*®. The mercury 


TABLE I. Transition parameters of bismuth. 








PRES- TeM- 





























SURE PERATURE dr Ar LATENT HEAT 
(kg ‘cm?) (°C) dp (cm*/g) (kgem/g) (gcal./g) 
Lieuip-I 
1 271.0 —0.00342 0.00345 —549 — 12.87 
5,000 251.9 418 386 485 11.37 
10,000 228.8 506 419 416 9.76 
15,00Q 200 61 440 346 8.11 
I-II 
18,000 170 —.196 — .00480 108 2.53 
19,000 150 —.196 485 104 2.44 
24,100 50 —.196 475 78 1.83 
29,200 —50 —.196 455 $1 1.20 
II-III 
22,800 170 —.298 — .00300 445 1.04 
23,600 150 —.298 .00310 44 1.03 
26,900 50 —.298 00310 33 77 
30,300 —50 —.298 .00275 20 47 
IlI-IV 
43,000 —100 +.3(?) —.0003(?) —.2(?) — .0047(?) 
43,500 + 50 3(? .0004(?) -4(?) — .0094(?) 
APPROXIMATE TRIPLE POINTS 
L-I-Il 
17.300 183 L-I —.0066 .0045 —310 —7.27 
I-II —.196 —.0047 110 2.58 
L-II +.00046 —.0002 — 200 —4.69 
L-II-IlI , 
22,400 185 L-II +.0004(?) -—.0002(?) -—200(?) —4.69(?) 
II-III —.298 —.0029 45 1.05 
L-III +.0091(?) —.0031(?) —155(?) —3.64(?) 
I-II-III 
32,300 -—110 I-II —.196 .0043 36 .84 
II-III —.298 .0025 14 .33 
I-III —.221 .0068 50 1.17 
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TABLE II. Transition parameters of mercury. 





250" 


200" ae 

150° | a a ce On 
a) ow 

100° | | 










































































nn “TT, 

50" + i | | 

| Y\] | 
RN tthe | 
\ : 

-507 

h let 
-100 | | <---|---t- ane 
10,000 20,000 30,000 40,000 50000 

Bismuth 


Fic. 3. The phase diagram of bismuth, pressure in 
kg/cm? as abscissa against temperature in degrees Centi- 
grade as ordinate. The arrows on the transition line III-IV 
mark the limits_of indifference. 
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Fic. 4. The volume changes for the various transitions 
of bismuth, temperature in degrees Centigrade as abscissa 
against volume change in cm*/g as ordinate. 


was frozen directly into the apparatus without 
protecting sheath. Determinations at this tem- 
perature were made with increasing pressure, 
decreasing pressure, and then with increasing 
pressure again. Temperature was then raised to 
0°C and the pressure to about 36,000 to maintain 
the mercury solid, and the transition located with 
decreasing pressure. The run was terminated by 
leak of liquid mercury past the piston due to 
melting, in spite of the fact that the pressure in 
the apparatus was well above the melting 
pressure at this temperature. The change of 
volume of the transition was moderately large, 
corresponding to the total volume decrement 
produced by an increase of pressure of about 
2000 kg, but nevertheless it was not sufficiently 
large. In conjunction with the super- and sub- 











TEm- 

PRESSURE PERATURE dr Av LATENT HEAT 
(kg/cm?) (°C) dp (cm*/g) (kgcm/g) (gcal./g) 
10,000 — 109 0.0077 0.00104 22.1 0.518 
15,000 — 73 67 .0007 1 21.1 495 
20,000 — 43 57 51 20.6 .483 
25,000 — 17 48 41 21.8 511 
30,000 + 5 40 36 25.0 .586 
35,000 + 23 32 32 29.6 694 








pressing which the transition will support, it was 
not found possible to land exactly on the equi- 
librium pressure, but the transition, when it 
started, always ran to completion. The apparatus 
was then set up again, and the measurements at 
—79° repeated, with good checks for both 
pressure and change of volume. Temperature was 
then raised to —43°, where it was maintained 
constant by continual manual manipulation by 
two assistants. After determination of the transi- 
tion parameters from above and below, temper- 
ature was raised to +22°, and again two deter- 
minations made. No attempt was made for 
higher pressures because there was slight leak 
due to melting at +22°. At room temperature 
no other transition was found up to 45,000. A 
notable feature was the very markedly lower 
compressibility of the high pressure modification. 
The transition parameters are given in Table II, 
and the experimental points are shown in Fig. 5. 

Extrapolation of the transition curve to atmos- 
pheric pressure suggests that the transition 
should run in the neighborhood of liquid-air 
temperature. Such a transition has not been 
reported, but should now be searched for. There 
are measurements of the x-ray structure down to 
liquid-air temperature; the best recent work is in 
accord about the structure and finds no change 
with temperature. There was an early deter- 
mination that gave a discordant structure, but 
apparently even this cannot be explained in 
terms of a transition. In all these x-ray determi- 
nations the mercury was, of course, frozen at a 
higher temperature and then cooled to liquid-air 
temperature; under these conditions the transi- 
tions may have been suppressed, even if the 
temperature of the transition had been slightly 
exceeded, because of viscous resistance to the 
transition at low temperatures. If the mercury is 
cooled sufficiently beyond this temperature, 
however, the transition should run. There are 
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Fic. 5. The phase diagram and the volume change of the 
transition of mercury. Pressure in kg/cm? as abscissa, 
against temperature in degrees Centigrade as left-hand 
ordinate scale, and volume change in cm*/g as right- 
hand ordinate scale. The pressure-temperature points are 
shown as open circles, and the volume changes as filled-in 
circles. 


measurements of the specific heat and the elec- 
trical resistance of mercury all the way from 
hydrogen temperatures. For the specific heat 
there are measurements by Simon‘ and by Pol- 
litzer.5 The measurements of Simon are very 
much the smoother; if these are plotted against 
temperature a very distinct break in the tangent 
will be found between 53° and 68°K. Simon does 
not comment on this, but it seems to me too 
pronounced to be spurious, and may perhaps be 
evidence of the transition. The resistance of 
mercury has been carefully measured at Leiden ;° 
there is a gap in the temperature range from 
—156° to —183°. If the deviations from line- 
arity of resistance as a function of temperature 
are plotted, a break will be found in this tem- 
perature interval, which again appears to be 
beyond experimental error. However, I do not 
regard this as definite as the specific heat 
evidence. 

There is another consideration; it is possible 
that the transition curve found above does not 
extend all the way to atmospheric pressure, but 
it might be terminated by a triple point with 
still another modification somewhere below 
—79°, and with a low temperature transition 
line running steeply enough not to cut the tem- 


4F. Simon. Ann. d. Physik 68, 241 (1922). 

5 Pollitzer, Zeits. f. Elekchemie u. angew. physik 
Chemie 19, 513 (1913). 
- Onnes and Holst, Comm. Phys. Lab. Leiden, No. 142a, 
1914. 
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perature axis at atmospheric pressure. Such a 
state of affairs would demand the existence of 
yet another transition line, with a transition 
pressure at room temperature somewhere be- 
tween 35,000 and the melting pressure, 12,000. 
Such a transition could not be observed at room 
temperature with the 50,000 apparatus because 
of leak, and certainly has not been observed with 
my other pressure apparatus. But at lower tem- 
peratures such a transition should be observable 
with the former apparatus. Accordingly a special 
examination was made for a new transition at 
— 79° with the regular 12,000 apparatus. Pressure 
was now transmitted with a true liquid, and the 
indications were very much more sensitive than 
with the 50,000 apparatus. No transition was 
found up to 10,000 kg/cm’. It therefore seems 
to me highly probable that the transition curve 
discussed above extrapolates to atmospheric 
pressure, and if the transition is not entirely 
suppressed by viscosity it will be found at 
atmospheric pressure not far from liquid-air 
temperature. 


Thallium 


The material was some which I had inherited 
from the late Professor T. W. Richards, and had 
been highly purified. Four different set-ups were 
used. In the first two the thallium was allowed 
to come directly in contact with the steel. The 
first run gave successful readings at 0° and 35°, 
reaching at 0° a maximum pressure of 51,000. 
There is a high pressure transition; the band of 
indifference is fairly wide, and although it was 
possible to locate the beginning of the transition 
with increasing pressure it was not possible to 
get a good value for Av with increasing pressure, 
because it is obscured by plastic yield in the 
cylinder. The second set-up gave good readings 
at CO, temperature, and again Av only with 
decreasing pressure. At 0° the run was terminated 
by the thallium blowing out past the piston. The 
third set-up gave good readings at 0° and 75°, 
the transition having become so sharp at 75° as 
to give good values for Av both with increasing 
and decreasing pressure. At 180° the cylinder 
broke, amalgamated by the thallium, at a 
pressure which later proved to be almost exactly 
a transition pressure. It has been my rather 
common experience that rupture due to the 
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Fic. 6. The phase diagram of thallium, pressure in kg/ 
cm? as abscissa against temperature in degrees Centigrade 
as ordinate. 


amalgamation effect is quite likely to take place 
while a transition is in progress; perhaps there is 
greater mobility during the transition, or perhaps 
the new phase has an intrinsically greater pene- 
trating effect. The fourth and last set-up was 
with the thallium completely surrounded with a 
copper capsule to protect the steel. With this 
set-up a point was successfully determined at 
170°, which proved, however, on plotting to have 
both too low a pressure and too low a change of 
volume to fit in with the other points. Thallium 
is known to be polymorphic at atmospheric 
pressure, and the transition line has been fol- 
lowed by Werner’ up to 3000 kg/cm*. This 
transition is of the abnormal ice type, tem- 
perature falling with increasing pressure. My 
point at 170° fits almost exactly on a linear pro- 
longation of Werner’s curve, and there can be 
little doubt but that this is the transition I-II. 


TABLE III. Transition parameters of thallium. 

















Tem- 

PRESSURE PERATURE dr Ar LATENT HEAT 
(kg/cm?) (°C) dp (cm3/g) (kgcem/g) (gcal./g) 
I-II 

1 227 0.00004 «fF 
—0.00187;5 
39,000 153 .00024 —41.7 —0.98 
II-III 
43,000 — 100 — .00039 + 1.1 + .026 
— .0632 
39,000 +153 —.00053 + 2.7 + .063 
I-III* 
39,000 153 + .00242 — .00029 —39.0 — 91 








* The parameters for I-III are at the triple point I-III and are 
calculated. 
t+ Werner's value. 


7M. Werner, Zeits. f. anorg. Chemie 83, 275 (1913). 
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Fic. 7. The volume changes for the transitions of thal- 
lium, temperature in degrees Centigrade as abscissa against 
volume change in cm*/g as ordinate. The point at 227° 
for the transition of I-II is the value of Werner. 


The phase diagram demanded by these con- 
siderations is shown in Fig. 6, and the changes of 
volume in Fig. 7. The II-III line, obtained by 
calculation from the thermodynamic parameters 
which are given in Table III, is shown dotted. 
At 170° the I-III transition was not found up to 
45,000; this maximum at this temperature was 
set by plastic yield in the cylinder. It will be 
seen that this is consistent with the calculated 
I-III line. 

Thallium was also examined by the shearing 
stress method at room temperature. The shearing 
strength curve has marked changes of direction 
in it, and there can be no doubt of the existence 
of other forms. It is not easy to determine from 
the shearing curve alone, however, whether the 
proper interpretation of it is one or two transi- 
tions. If the results of this volumetric analysis 
are accepted as meaning that there cannot be 
more than one transition, then the shearing 
experiments give 25,000 for the mean pressure of 
the transition. This is much lower than the value 
found above, and indicates an unusually large 
displacement of transition pressure by shearing 
stress. ' 


Tellurium 


The material was obtained from Kahlbaum; 
it was formed into a single crystal by my method 
of slow solidification, which should have effected 
an additional purification. This substance has a 
transition at room temperature in the neighbor- 
hood of 45,000 with a very large volume change, 
and at lower pressures probably another transi- 
tion with volume change so much smaller that 
there is not complete certainty about it. Runs 


_ were made with three different set-ups. The first 


showed no transition to 52,000 at CO, tempera- 
ture. A transition with unusually large volume 
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change was then located at 41,000 at 145°, but 
the volume change was not measured. The 
volume change was so large that it was easy to 
obtain pressure readings on the two-phase 
system. At 90° the experiment was terminated 
at 47,000 by fracture of the piston before the 
transition had put in an appearance. The second 
set-up gave at 171° an equilibrium point at 
41,000. By cautious manipulation, what ap- 
peared to be a satisfactory value for the change 
of volume was also obtained at this temperature. 
The difficulty is that because of indifference and 
friction the transition runs to completion only at 
pressures very near the maximum attainable 
with the apparatus. At 142° a good value was 
obtained for the equilibrium pressure, but the 
indifference and friction were here so much 
greater that the piston was fractured in trying to 
force the reaction to run to completion in order 
to determine Av. All the behavior at high pressure 
apparently indicates that we are here dealing 
with a transition much like that between ice I 
and III. That is, there is very little latent heat, 
and the transition may run with almost explosive 
rapidity at high temperatures, but there is an 
enormous temperature coefficient of transition 
velocity, so that at lower temperatures the transi- 
tion will not run perceptibly under practical con- 
ditions. 

After the measurements just described, explo- 
rations were made with the shearing apparatus 
at room temperature; this gave large and striking 
breaks, the only interpretation of which seemed 
to be two transitions. It is highly probable that 
the transition at the higher pressure, near 50,000, 
is the transition described above, the viscous 
resistance at room temperature being overcome 
by the shearing forces, as is not unusual. The 
transition at lower pressures seemed new. The 
volumetric apparatus was accordingly set up 
again, this time with the tellurium enclosed in a 
tin sheath, not to protect the steel, which was 
not necessary because there seems to be no 
chemical attack, but to reduce friction, which is 
unusually great with tellurium and makes the 
results somewhat irregular. At room temperature 
a doubtful transition with volume change near 
the limit of the detectible was found with in- 
creasing and decreasing pressure at a mean 
pressure of 36,000. At 85°, doubtful results were 
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Fic. 8. The phase diagram of tellurium, pressure in 
kg/cm? as abscissa against temperature in degrees Centi- 
grade as ordinate. 


again obtained, the readings during increasing 
pressure being obscured by an effect which later 
proved to be due to slight leak of the tin. The 
results at 100° of the original set-up were now 
reexamined, and a slight break found at a 
pressure which extrapolates linearly from the two 
just found. I therefore think that the location of 
this transition is established sufficiently well to 
justify putting it in the phase diagram, but it 
must be remembered that it is not entirely free 
from doubt. The Av values are too uncertain to 
more than estimate; they cannot be much larger 
than 0.0001 cm?/g. 

The transition parameters are indicated with 
sufficient accuracy in Fig. 8. In addition to the 
pressure-temperature values shown there, Av 
II-III at 171° is 0.0058 cm*/g. 


Gallium 


A high pressure transition would be expected 
here for the same reason as for bismuth, 
namely, because the melting curve is of the ice 
type. The transition was at once found at CO» 
temperature at about 20,000; the change of 
volume is large, readings are easily made with 
the two phases present, and repetition gave good 
checks. Two other points were then obtained at 
—41° and —11°. No other transition was found 
to 46,000. The transition line is of the ice type 
and intersects the extrapolated melting line only 
slightly beyond the limit of my previous measure- 
ments.* It was tempting to try for points on the 
melting line of II; this was not possible with the 
50,000 apparatus, but the pressures are so low 
as to be within the reach of previous types of 
apparatus. A special cylinder was therefore con- 


8 P. W. Bridgman, Proc. Am. Acad. 56, 104 (1921). 
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structed, much like that used in my measure- 
ments on the alkali metals to 20,000, only 
somewhat smaller. I was fortunate to obtain 
from Dr. W. M. Craig a supply of 16.5 grams of 
highly purified gallium. With this, measurements 
were made with the new apparatus by the con- 
ventional piston displacement method, pressure 
being transmitted by isopentane, and measured 
by a manganin pressure gauge in the liquid, so 
that all frictional effects were eliminated. 

A large number of readings were made with 
the special apparatus. At first there appeared to 
be certain very puzzling inconsistencies. The 
explanation was presently found to be that there 
is still another high pressure modification, II’, 
very much like II in its thermodynamic proper- 
ties, but totally unstable with respect to it. 
Points were obtained on both transition lines 
I-II and I-II’, but only on the one melting line 
L-II’. The points on the melting line were ob- 
tained by the procedure of raising temperature at 
constant volume. The reason for this was that 
the preliminary, and as it later proved, incorrect 
values for the change of volume on melting led 
to the expectation that this melting line would 
be nearly horizontal, and in this case the ordinary 
procedure would have been useless. The pro- 
cedure of changing temperature gives much less 
accurate values than the other. The values so 
obtained checked within the expected accuracy 
with the values calculated from the data for the 
other transitions, so that it did not seem worth 
while to repeat for the improvement in accuracy 
that could doubtless have been obtained by the 
other method. Although no points were obtained 
on the L-II line because of some caprice with 
which the phases appeared, the existence of a 
melting line lying higher than L—II’ was checked 
on one occasion by carrying the solid phase II 
a degree or more above the melting temperature 
of II’ without melting. In view of the fact that a 
solid practically never superheats, this consti- 
tutes almost a certain check of the explanation 
in terms of a third absolutely unstable phase. 

Advantage was taken of the opportunity to 
measure the changes of volume on the melting 
line of I, thus making possible a calculation of 
its thermodynamic parameters. This had not 
been previously possible because the quantity of 
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gallium then available was too small. Over the 
first 6000 kg the new determinations of the 
melting temperature of I as a function of pressure 
agree almost exactly with the former values, but 
above 6000 the new temperatures are somewhat 
lower, the maximum difference being 1° at 12,000. 

The results are shown in Figs. 9 and 10, and 
the parameters of the transition are collected in 
Table IV. The value for Av at atmospheric 
pressure shown in Fig. 10 is the value of Boyer, 
0.00529 cm*/g; it is seen to lie almost exactly on 
the curve with my high pressure values. No 
parameters are given in the table explicitly for 
the melting of II and II’ except at the triple 
points; experimental accuracy did not justify the 
attempt to estimate the variation of these 
parameters along the melting lines. The initial 
calculated slopes of these lines are indicated with 
sufficient accuracy in the figure. 

Gallium is only the second example I have 
encountered in which a totally unstable phase 
appears at high pressures. The other example is 
water (both H,O and D,.Q). Perhaps it is sig- 
nificant that in both these cases the ordinary 
melting curve is of the ice type. 


TABLE IV. Transition parameters of gallium. 








TeM- 


























PRESSURE PERATURE dr Ai LATENT HEAT 
(kg/cm?) (°C) dp (cm'/g) (keem/g) (gcal./g) 
L 
1 29.85 —0.00205 0.00529 —782 — 18.3 
2,000 25.7 210 555 790 18.5 
4,000 21.5 215 585 803 18.8 
6,000 17.1 221 625 822 19.2 
8,000 12.6 234 665 813 19.0 
10,000 7.7 252 705 785 18.4 
12,000 2.5 270 740 755 17.7 
I-ll 
12,230 0 — 0085 — .0098 +314 7.4 
14,580 —20 — .0O88 98 281 6.8 
16,800 —40 — 0095 98 240 5.6 
18,800 —60 — .0108 OR 193 4.5 
20,500 —0 — .0130 98 145 3.4 
I-II’ 
12,800 0 — 0098 +314 7.4 
— 0.0852 
13,390 —- 5.0 
TRIPLE Points 
L-I-Il 
12,050 2.4 L-I — .00270 00740 —755 —17.7 
I-Il — .00852 — .0098 +317 + 7.4 
L-Il + .0O151 — .0024 —438 — 10.3 
L-I-II’ 
12,750 04 L-I — .00277 .0075 —740 —17.3 
I-III’ — .0O852 — 0098 +314 + 7.4 
L-Il’ + .00137 — .0023 —426 -— 99 
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Fic. 9. The phase diagram of gallium, pressure in 
kg/cm? as abscissa against temperature in degrees Centi- 
grade as ordinate. 
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_ Fic. 10. The volume changes for the transitions of gal- 
lium, temperature in degrees Centigrade as abscissa 
against volume change in cm*/g as ordinate. The filled-in 
circles are for the transition I-II’. 


Iodine 


This was found to have a transition at CO» 
temperature with small but unmistakable change 
of volume. The pressure limits of the transition 
are unusually wide; the transition runs to com- 
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pletion at once with increasing pressure at 20,600 
kg/cm’, and on decreasing pressure the transition 
runs at once to completion on reaching 10,800. 
The run was repeated with complete check both 
for the pressure limits and for the volume change, 
which is between 0.004 and 0.005 cm'/g. A run 
was then made at 40°, and no perceptible differ- 
ence was found in the pressure limits or in the 
change of volume. No other transition was found 
up to 53,000 at either temperature. A second 
set-up was then made, and readings made on a 
time schedule twice as long as on the first set-up, 
two minutes between readings instead of one, in 
the endeavor to narrow the pressure range. This 
was moderately successful, the pressure range 
now being 6000 instead of 9800, and the mean 
pressure 12,500 instead of 15,700. The volume 
change was essentially the same as before. It is 
therefore probable that the transition line of 
iodine runs nearly vertically, that is, with small 
latent heat, at a pressure not far from 13,000, 
and a volume change not far from 0.0045 cm*/g. 
The shearing apparatus gave both on the first 
trial and on repetition an unmistakable discon- 
tinuity with increasing pressure at about 20,000, 
which is doubtless the transition just described. 
But a second set-up showed practically nothing 
with increasing pressure and only a minor break 
with decreasing pressure. This emphasizes the 
fact that the shearing measurements give only 
qualitative indications of transitions. 


SUMMARY OF EXPLORATIONS FOR THE ELEMENTS 
GIVING NEGATIVE RESULTS 


There follows now a brief summary of these 
elements for which negative results were obtained 
by the volumetric method. It must be again 
emphasized that the sensitiveness is not as great 
as is possible at low pressures with the former 
apparatus, in which the pressure is transmitted 
by a true liquid. The greatest difficulty arises 
from friction; this results in irregularity of the 
readings, the motion of the piston having a 
tendency to take place in jumps. These jumps 
simulate small transitions. What is more, some- 
times the jumps repeat themselves both with 
increasing and decreasing pressure. In these cases 
it is particularly difficult to be sure that small 
transitions have not been missed. However, I do ° 
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not believe that on the average any transition 
with volume change much larger than 0.001 
cm*/g would have escaped detection. 

The arrangement in the following is by columns 
in the periodic table. 

Lithium. Nothing to 46,000 at CO» tempera- 
ture. 

Potassium. Nothing to 52,000 at CO, tem- 
perature. 

Caesium. Probably nothing to 52,000 at 0°. 
This is more doubtful; there is some evidence for 
a transition at 25,000. Discontinuities were found 
four times both with increasing and decreasing 
pressure, but the uncertainty is particularly large 
here because of the very great compressibility of 
caesium, and I am inclined to think that the dis- 
continuities were the result of friction. 

Gold. Nothing to 51,500 at room temperature. 
The compressibility of gold stands strikingly out 
of line with that of copper and silver, and it 
therefore seemed to me not unlikely that a 
transition might be found. 

Magnesium. Nothing to 51,000 at CO, tem- 
perature nor to 39,000 at 150°C. 

Calcium. Nothing to 51,000 at CO, tempera- 
ture. 

Cadmium. This is known to have two new 
modifications at room temperature with very 
small volume change.’ The transitions are very 
hard to detect, and up to now have been found 
only with single crystal material; the pressure of 
the two transitions is nearly the same, and in the 
neighborhood of 5000. Runs with the new volu- 
metric apparatus were made at CO, temperature. 
It is fairly certain that there was a transition with 
increasing pressure near 20,000. The same result 
was found on repetition, but the transition could 
not be found with decreasing pressure. It is not 
unlikely that this is the transition already found 
at room temperature, made highly viscous by the 
low temperature. If this is the explanation it 
means that the reverse transition ran so close to 
atmospheric pressure that it could not be disen- 
tangled from the various creep effects always 
present at the zero on release of pressure. No 
other transition was found to 44,000 at CO» 
temperature. It would probably pay to examine 
this metal further at other temperatures. 


9 P. W. Bridgman, Proc. Am. Acad. 60, 346 (1925). 


Barium. Nothing to 51,000 at CO», tempera- 
ture nor +150°. 

Indium. Nothing to 50,000 at room tempera- 
ture. The crystalline structure of indium is 
tetragonal, rather unusual in a metal; a transition 
would not have been surprising. 

Carbon, graphite. Nothing to 52,000 at COs: 
temperature nor 155°. 

Germanium. Nothing to 52,000 at room tem- 
perature; the friction was unusually high. 

Tin. Probably nothing to 51,000 at CO, and 
room temperature; the shearing experiments 
suggest a transition. 

Lead. Nothing to 52,000 at CO, temperature 
nor at 150°. 

Phosphorus, dense variety of black phosphorus. 
At CO, temperature nothing was found to 52,000; 
the cylinder broke at the maximum. The second 
set-up gave nothing at 150° to 36,000, where the 
cylinder broke. In both cases the break was 
evidently due to chemical attack. 

Arsenic. Two independent set-ups were used. 
The first gave nothing to 50,000 at 140° and the 
second nothing to 52,000 at 20° and 100°. A 
transition was expected because of the similarity 
of the crystal structure to that of bismuth. 

Antimony. Nothing was found to 52,000 at 
room temperature or at 125°. A transition was 
expected because of the similarity of bismuth and 
antimony. However, it was found in the shearing 
experiments that there is very probably a 
transition near 50,000, which may have been 
suppressed in these volumetric experiments by 
viscosity. 

Sulfur, amorphous modification. This showed 
nothing to 40,000 at room temperature, where the 
cylinder broke, doubtless by chemical attack, 
since the sulfur was directly in contact with the 
steel. 

Selenium. The amorphous variety was exam- 
ined at CO, temperature to 44,000. There were 
two rather large jumps during increasing pressure 
between 25,000 and 30,000 in the direction of 
increasing density, and one would therefore 
suspect an irreversible transition to the metallic 
form. However, a direct determination of the 
density after the run showed no permanent 
change. The matter should be examined further. 

The metallic form showed nothing to 46,000 
at CO temperature nor at +155°. 
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OTHER ELEMENTS PROBABLY POLYMORPHIC AT 
H1GH PRESSURE 


There is repeated here for the sake of com- 
pleteness a list of elements which the method of 
shearing described in a previous paper has shown 
to be probably polymorphic at high pressure, 
namely: lithium, calcium, vanadium, manganese, 
zinc, strontium, yttrium, antimony, barium, 
lanthanum, cerium, praseodymium, erbium, and 
thorium. 

Some of these elements appear in the list 
which have given negative results by the volu- 
metric method of this paper. It is probable that 
in such cases the volume change may be too small 
to detect. The other elements which have given 
positive indications by the shearing method have 
not yet been investigated by the volumetric 
method. 


RESULTS FOR COMPOUNDS 


We now pass from the elements to consider a 
few simple compounds. The systematic inves- 
tigation of polymorphism among the compounds 
will evidently be very time consuming; I am 
planning to take up this investigation in the near 
future. Most of my activity in this preliminary 
investigation was confined to the alkali halides. 
One or another of these salts was investigated a 
good many times in the course of the develop- 
ment of the apparatus. One would suspect that 
by the application of sufficient pressure the 
potassium, sodium, and lithium salts could be 
made to undergo transition from the simple 
cubic NaCl structure to the body-centered CsCl 
structure, since this transition has already been 
found in the rubidium halides at moderate 
pressures. The expected transition was found 
almost at once with the potassium salts at 
pressures in the neighborhood of 20,000, but 
apparently pressures higher than 50,000 will be 
required for the transitions of the sodium and 
lithium salts. . 

KCl. Clear single crystal material was used. 
Measurements have been made on this substance 
with five different set-ups, and in spite of this 
the data are not as complete as for the other two 
potassium halides. Three cylinders ruptured 
during the measurements; the rupture of one was 
almost certainly due to the chemical action of 
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the metallic sodium in which the NaCl was 
embedded in order to make pressure more nearly 
hydrostatic. The other ruptures I think were 
probably just bad luck, due to damage received 
in previous exposures. One cylinder ruptured in 
which the KCI was protected from contact with 
the steel with a sheath of lead, and in the other 
no protection was used. Four of the five set-ups 
were made with the same type of cylinder as that 
used for the first measurements with bismuth; 
the fifth was an early form of the small (6 mm 
piston) apparatus, from which the stopper shaped 
cylinder later evolved. 

The experimental results are shown in Fig. 11. 

KBr. Clear single crystal material was used. 
Two independent sets of measurements were 
made, at six temperatures altogether, with the 
same cylinder, of the double piston bismuth type, 
with no accident. The results are shown in Fig. 
12. The regularity of the results is, I think, sur- 
prising, in view of the very great difference 
between the readings obtained with increasing 
and decreasing pressure, due both to the friction 
of the pistons, and the intrinsic width of the 
region of indifference itself. This latter, judging 
from comparison with other transitions, must be 
very much the more important factor at low 
temperatures. At room temperature there was a 
difference in the transition pressures of 13,600 
kg/cm? when approached from above and below. 
This decreases irregularly with increasing tem- 
perature; at 235° it had dropped to about 3000. 

KI. In addition to the first exploratory runs, 
when the existence of transitions in this series 
of compounds was first found, runs were made 
at four temperatures with the double piston 
bismuth apparatus, and at room temperature 
with the same early form of the 6 mm piston 
apparatus as used with KCl. The results are 
shown in Fig. 13. The pressure difference between 
increasing and decreasing readings is about the 
same as for KBr, being 13,400 at room tem- 
perature, and dropping off with increasing tem- 
perature, but perhaps not as rapidly as for KBr. 
The great width of the band of indifference 
explains why the transition was not found before 
at room temperature in an exploration up to 
20,000.” 

Comparison of the results for the six rubidium 


1” P. W. Bridgman, Zeits. {. Krist. 67, 371 (1928). 
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Fic. 11. The transition parameters of KCl. 
Temperature in degrees Centigrade as abscissa 
against pressure in kg/cm? as left-hand ordi- 
nate, with the corresponding points shown as 
open circles, and volume change in cm'/g as 
right-hand ordinate, with the corresponding 
points as filled-in circles. Notice that the 
usual order of the pressure-temperature co- 
ordinates is here reversed. 


and potassium salts shows some interesting re- 
semblances. At temperatures below 200° the 
order of the transition pressures is the same in the 
potassium and the rubidium series, the transition 
pressure increasing from the iodide to bromide to 
chloride. This is the reverse of the order in the 
series NH,Cl, NH,Br, NH,I. At 0° the pressure 
difference between iodide and chloride is not 
dissimilar in the potassium and rubidium series. 
Above 200° the transition curves for potassium 
iodide and bromide cross; in the rubidium series 
there will probably also be crossing, but at higher 
temperatures. 

The important thermodynamic parameters of 
the transitions at 130° are collected in Table V. 


TABLE V. Transition parameters of the potassium halides 











at 130°C. 
FRAC- CHANGE 
TIONAL OF 

DECREASE INTERNAL 

PRESSURE AV OF dr LATENT HEAT ENERGY 
SaLt (kg/cm?) (cm’/g) VOLUME dp (kgem/g) (g cal./g)(kg cm/g) 
KCI 19,900 0.055 0.108 +0.200 —84 —1.97 + 1006 
KBr 19,300 0325 .102 —0.188 +52 +1.22 680 
kl 17,850 026; .087 +0.430 —18 — .42 472 








The density of the ordinary modification at the 
pressure of transition was obtained by extra- 
polating the second degree formulas which have 
been found for the compressibility in the range 
up to 12,000.° The corresponding fractional 


changes of volume at the transition in the rubi- 


uJ. C. Slater, Proc. Am. Acad. 61, 135 (1926). 
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Fic. 12. The transition parameters of KBr. 
Temperature in degrees Centigrade as ab- 
scissa, against pressure in kg/cm? as left-hand 
ordinate, with the corresponding points shown 
as open circles, and volume change in cm*/g as 
right-hand ordinate, with the corresponding 
points as filled-in circles. Notice that the 
usual order of the pressure-temperature co- 
ordinates is reversed. 


dium series, from chloride to iodide was: 14.6, 
11.3 and 10.7 percent. The order is thus seen 
to be the same in the rubidium and _ potas- 
sium series. The numerical magnitudes are some- 
what less in the potassium series; this is to be ex- 
pected because of the much higher pressure of the 
transition. 

The latent heats are seen to be only a small 
part of the total energy change during the 
transition. The order of energy change is the 
same in the potassium and rubidium series; the 
absolute magnitudes are several-fold greater in 
the potassium series, as is also the range from 
chloride to iodide. 

Other compounds. This covers the compounds 
the transitions of which have been measured up 
to the present with the 50,000 apparatus. In 


19,000 





20 


0° ! 200 300 
KI 


Fic. 13. The transition parameters of KI. Temperature 
in degrees Centigrade as abscissa, against pressure in 
kg/cm? as left-hand ordinate, with the corresponding points 
as open circles, and volume change in cm*/g as right-hand 
ordinate, with the corresponding points as filled-in circles. 
Notice that the usual order of the pressure-temperature 
coordinates is reversed. 
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addition, a number of explorations have been 
made with negative results. Naturally the 
greatest amount of effort was spent in trying to 
find the transition in the sodium series. NaCl 
gave nothing at room temperature with two 
different set-ups, the first to 32,000 and the 
second to 50,000. Nal is the member of the series 
that would be expected to show the transition 
most readily; this gave negative results at room 
temperature with four different set-ups, to 
46,000, 49,000, 48,000, and 57,000, respectively. 
However, it has already been mentioned in the 
previous paper that the shearing curve of Nal 
begins to bend over at 50,000, and I believe the 
transition is not far away. 

Another interesting salt of the alkali halides 
series is CsF; this crystallizes with the NaCl 
structure, although all the other caesium salts 
are body-centered cubic. One might expect to be 
able to force the transition to the body-centered 
structure by high pressure. Two separate set-ups 
with CsF, for which I am again indebted to Dr. 
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R. W. G. Wyckoff, gave negative results, the 
first to 42,000 at room temperature, and the 
second to 56,000 at 137°. It would thus appear 
that the fluorides are qualitatively different from 
the other halides. In view of the failure to obtain 
the transition with CsF, it seemed hardly worth 
while to explore the fluorides of any of the other 
alkali metals. 

In addition to these, sugar gave negative 
results to 40,000 at room temperature, and BaS 
to 39,000 at room temperature. PbTe, on the 
other hand, probably has a small transition at a 
mean temperature of 18,000, but this is not 
certain; there is no further transition up to 
39,000. The explorations mentioned in this 
paragraph were made with a preliminary form of 
apparatus, only one-third as sensitive as the final 
form. 

I am indebted to the Rumford Fund of the 
American Academy of Arts and Sciences for 
financial assistance in purchasing some of the 
supplies. 
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Quantitative information as to the nature of intranuclear 
forces has been obtained by considering the hydrogen and 
helium isotopes, using a Hamiltonian of the exchange type 
based on the neutron-proton model with interactions be- 
tween all the particles. The interactions between like 
particles are taken to be the same for protons and for 
neutrons except for the small Coulomb force between 
protons which accounts for at least half of the difference 
in the binding energies of He* and H‘*. It is not inconsistent 
with existing experimental evidence’ to assume the range 
of such forces to be about the same as the range of neutron- 
proton forces. By making these assumptions and by using 
simple attractive potentials, it is found that the forces 
between like particles are given by a potential with depth 
not greater than 30 mc? and with an effective radius of 
action which must be less than 2.8 10~'* cm to be con- 
sistent with experimental data for the binding energies. 
A more accurate determination of the constants based on 
assuming a close correlation of the three- and four-body 
problems with “equivalent” two-body problems gives the 
depth of the proton-neutron potential as 74 mc? and of the 


* Now at University of Wisconsin. 
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like-particle potentials as 26 mc? and a range of 2.2 10°" 
n for the forces. Reasons for believing in the reliability 
of Nhe “equivalent” two-body method for determining 
imate binding energies for three- and four-body 
problents are discussed. It is shown furthermore that if the 
rators multiplying the neutron-proton inter- 
action potenNals are assumed to be linear combinations of 
the Majorana ahd Heisenberg types, it is possible to explain 
the large scattering of neutrons on protons at low velocities. 
The proportion of the Heisenberg to the Majorana operator 
in the linear combination necessary to obtain the correct 
scattering is about one-fifth. It is then shown that in the 
three- and four-body problems such an operator can to a 
good approximation be represented by an equivalent 
Majorana operator and that because of this fact the results 
obtained for nuclear energies and for the magnitude and 
range of the forces (in which Majorana operators were 
used) remain the same except that the depth of the poten- 
tial for like particles is increased from 26 to 41 mc. The 
latter value yields a scattering intensity for protons in 
hydrogen which agrees qualitatively with the experimental 
results. 


exchange 
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SECTION I. INTRODUCTION 


E propose to use the binding energies of the 

hydrogen and helium isotopes to fix the 
parameters in a simple exchange type nuclear 
Hamiltonian based on the neutron-proton model 
with interactions between all pairs of particles. 
Several authors have found evidence for large 
interactions between like particles in the odd- 
even properties of nuclei.! In a previous investiga- 
tion? the same problem was studied, but subject 
to the assumption, which we do not retain, that 
the interaction between like particles could be 
treated as a small perturbation in comparison 
with the neutron-proton interaction. 

The Hamiltonian has the form 


H=— 32 s— 2S (rein j) Pox —='K (rviv,) P’ viv; 
— DL (rein) P' xixj + (27)? ‘82/1 /reix; (1) 


with the units m,c? (510,000 ev) for energy and 
h/2xc(m.m,)'” (8.97 X10- cm) for length. There 
exist several possibilities for the exchange opera- 
tors P and P’ depending on the general character 
of the forces between the nuclear particles. If 
these forces are first attractive as two particles 
approach and then strongly repulsive upon closer 
approach, the exchange operators may be re- 
placed by unity. The inner repulsive regions pre- 
vent the collapse of heavy nuclei which would 
otherwise occur when the forces are of the 
ordinary (Wigner) type. Because of the mathe- 
matical difficulties connected with the correct 
treatment of the repulsive regions, this model has 
not yet been adequately investigated. 

If the inner repulsive regions are lacking, the 
burden of keeping the nuclear volume and bind- 
ing energy proportional to the total number of 
particles is put on the exchange operators. These 
may have the form’ 


Py, =(1—g) Pi" +2Pi", (2) 


a linear combination of the Majorana and Heis- 


1K. Guggenheim, J. de physique 5, 253 (1934); L. A. 
Young, Phys. Rev. 47, 972 (1935); H. Bethe, unpublished 
lectures. 

2 E. Feenberg, Phys. Rev. 47, 850 (1935), denoted by FI 
in text; Phys. Rev. 47, 857 (1935), denoted by FII in text. 
For discussion of related problems see L. H. Thomas, 
Phys. Rev. 47, 903 (1935) and W. V. Houston, Phys. Rev. 
47, 942 (1935). 

3 P;;™ interchanges the space coordinates i and j in a 
function to which it is applied; P;;” interchanges both space 
and spin coordinates. 


enberg exchange operators, and‘ 
Pi) =3(1+2P i") = — foi o;. (3) 


Because the wave functions for the hydrogen and 
helium isotopes are symmetrical in the space co- 
ordinates of like particles, we may replace P’ by 
unity for our purposes without loss of generality. 
It is known from the calculations in FI that the 
effective width and depth of the neutron-proton 
potential well required to account for the binding 
energies of the deuteron and the alpha-particle 
depend only slightly on whether the forces are of 
the Wigner or of the Majorana type. Thus, if g 
is small in comparison with unity, it is possible 
to compare directly the exchange type nuclear 
model for the hydrogen and helium isotopes with 
a model based on a Hamiltonian containing only 
ordinary forces. For this reason a calculation 
with the exchange model has a certain generality 
in that it gives at once quantitative information 
about the alternative model with attractive forces 
of the Wigner type and also yields a good qualita- 
tive picture of the magnitude of the attractive 
potential wells when these are associated with 
inner repulsive regions. 


Section II. POTENTIALS 


To obtain a manageable problem we set K(r) 
equal to L(r) and suppose that K(r) and J(r) 
are proportional. The first assumption finds some 
support in the fact that at least half of the small 
difference between the binding energies of H* and 
He’ can be accounted for as resulting from the 
Coulomb interaction between the protons in He®. 
The second assumption is arbitrary, but does not 
conflict with the experimental evidence as yet 
available. The definition of the Hamiltonian is 
completed by the assumption 


J(r)=A,,e-*", K(r)=L(r)=A,e*". (4) 


The experimental energies have the values® 
E(H2)2 —4.0, E(H5)2—15.8, E(He*)= —54.0. 
From the deuteron equation, HY=Ey with 


* In Eq. (3) o; represents the Pauli spin operator for the 
ith particle. The equivalence of the two expressions for 
P;;' can be deduced from a relation derived by Dirac in his 
Quantum Mechanics, first edition, p. 215, Eq. (35). See also 
Van Vleck, Phys. Rev. 48, 367 (1935) for a discussion of 
the operator P;;’. 

> Chadwick and Goldhaber, Nature 134, 237 (1934); 
Oliphant, Kempton and Rutherford, Proc. Roy. Soc. A149, 
406 (1935). 
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TABLE I. A,g(@) for EU?) = —4.0. 


4.35 4.60 4.85 








Ayr/a 345 3.66 3.88 4.11 

1/a'? O12 O15 O18 0.21 0.24 0.27 0.30 
Ayr/a 5.10 5.36 5.71 6.16 665 7.17 7.70 
l/e’? 0.33) (0.36 «60.40 «20.45 0.50 0.55 0.60 
E=-—4.0, A,, is obtained as a function of a by 


means of a series of numerical integrations: 
Table I makes possible the ready determination 
of the lowest eigenvalue of any two-particle 
Schrédinger equation with an arbitrary error 
function type of potential® and will be used for 
that purpose in Section IV. The function A,,(a@) 
does not depend on the value assigned to g. 
This is true because the wave function y for the 
deuteron is symmetrical in both the space and 
spin coordinates of the two particles and hence 
in this case 


((1—g)P”+gP"}y=y. (5) 


When, however, the spins of unlike particles are 
not completely parallel, the properties of the 
model do depend on g. In Sections III and IV 
the discussion is based on the simplifying assump- 
tion that the forces between unlike particles are 
of the pure Majorana type (g=0). The general 
case in which g differs from zero is considered in 


Section V. 

SeEcTION III. THE SIMPLE VARIATIONAL METHOD 
With the Gaussian type of wave function 
go(1 . 2, 3) = Nem (v!2) (ria? +113? )—(u/2) r23? (H), (6) 

go(1, 2; 3, 4) 


= 24. vec? -4-res2 )— a2 2 
= Nem (v/2) (ri3?+ria?+r23?+r24?)— (u/2) (ria +r3¢?) (He), (7) 


the average value of the Hamiltonian operator 
(omitting the Coulomb terms and setting g=0) 
can be written 


E°(H*) =(2+p)ao—2A,,{p(4—p)/(1+2p) }*” 
x (a/(o+1))??—A,(po/(po+1))*? (8) 
for the three-body problem and 
E*(He*) = (3/2)(2+p)ac—4A,,.(2p— p*)*” 
X (¢/(a+1))*?—-2A,(po/(po+1))*? (9) 
for the alpha-particle. These differ from the cor- 


‘FI, Eq. (4). 
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responding expressions in FI (Eqs. (18) and (19)) 
in the presence of the term proportional to A, and 
in a more suitable choice of the variation param- 
eters.’ The contribution from the Coulomb inter- 
action to the energy has the form 


AE=1(2apo)"? 


for both He*® and He‘. The values of p and o 
which minimize (8) and (9) must be found by 
trial. However the dependence on ? is slight, 
the best value generally being between 0.9 and 
1.0. 

Let A,(a, H*) and A,(a, He*) represent the 
functions defined by the conditions that the ac- 
curate lowest eigenvalues of the three- and four- 
body problems have the experimental values 
E(H’) = —15.8 and E(He*) = — 54.0, respectively. 
Eqs. (8) and (9) can be used to obtain upper 
bounds on A,(a, H*) and A,(a, He*). Consider 
first Eq. (9). There exists a function A,°(a, He‘) 
such that if A, is replaced by A,°(a, He*) in (9) 
the minimum value of E°He*)+A£ with re- 
spect to p and o coincides with the experimental 
energy. By the minimum property of the varia- 
tion method, we have 


A,(a, Het) =A," (a, He‘). (11) 


We find that A,(a, He‘) is practically constant 
with the value 30 mc? on the range 10=a=50. 
In the same way we find that A,°(a, H*) varies 
from the value 30 mc? at a=10 to 100 mc? at 
a=50. Furthermore the difference A,°(a, H*) 
—A,(a, H*) is much greater than the correspond- 
ing difference A,°(a, He*) —A,(a, He‘). This can 
be seen from the fact that the binding energy in 
the three-body problem is a small difference be- 
tween large potential and kinetic energy terms. 
Hence the energy is sensitive to errors in the 
wave function. The same statement is true for 
the alpha-particle, but in a much less extreme 
degree. Thus the correct value of A, is less than 
30 mc? and the correct value of a@ is greater than 
10 (1/a"!? <2.810- cm). 


(10) 


SecTion IV. THe “EQUIVALENT” Two-Bopy 
METHOD 


To improve upon these inequalities we have 
recourse to the method of the “equivalent’”’ two- 


7In Eq. (8) 4ac=Sv+y, 2acp=v+2y; in Eq. (9) 2ae 


=3v+y, acp=v+u. 
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body problem described and used in FI. The 
method is not rigorous, but its general correctness 
is supported by the considerations advanced in 
FI and by the results of attempts to improve on 
Eqs. (8) and (9) with more complicated varia- 
tional calculations (see appendix). The procedure 
involves first replacing p by unity in both (8) 
and (9) in order to reduce the two different poten- 
tial energy terms to the same form. In FI because 
there was only one kind of potential energy term, 
this simplification was unnecessary. The result- 
ing simplified expressions have exactly the form 
taken by the corresponding expression for the 
variational problem associated with the equation 


{d?/dr?+E’'+Be-*"} o=0 (12) 


when the approximate wave function g=e-”"”? 
is used to evaluate the energy integral. The formal 
analogy is made exact by the identifications 


B=2a sas B=3a 
B=2A,,+A, B=2(2A,,+A,) 


He‘. (13) 


These identifications establish a correspondence 
between the three and four-body problems and 
“equivalent” two-body problems which can be 
solved accurately. To proceed further we insert 
into (12) the experimental value for E’ and, using 
Table I, compute B and hence A, as a function 
of a for both H* and the alpha-particle. We 
designate these functions by the symbols 
A,'(a, H*) and A’(a, He*). Small corrections to 
make up for the inaccuracy resulting from setting 
p=1 are easily determined by means of Eqs. 
(8) and (9). Results are exhibited in Tables II 
and III. Table II includes a column giving the 
binding energy of H®* as a function of a when 
A,(a)=A,'(a, He*). At the point of intersection 
(A,’(a, H*)=A,’(a, He*)) ahas the value 17 (1/a'”” 
=2.17X10-" cm), A,,=74.3 mc? and A,=26 
mc. The figures given by Young! can be com- 
pared with the a=10 column in Table II by 
noting that 3(2A,,+A,) corresponds to the 
depth of the rectangular potential function used 
in Young’s calculation. There is satisfactory 
qualitative agreement. 

All the results of this and the preceding sections 
can be extended without difficulty to other forms 
of the attractive potentials, in particular, to the 
rectangular potential well and the straight ex- 


TABLE II. The “equivalent” two-body method for the deter- 
mination of A,'(a, He*). 








a Ayn A, —E’(He*)C.E.(He*)* —E’(H®) 





10 49.7 29 55.4 1.4 

20 84.4 24 55.6 1.6 15.2 
30 117.1 19 55.8 1.8 13.7 
40 148.8 12 56.0 2.0 12.8 








TABLE III. The ‘‘equivalent"’ two-body method for the deter- 
mination of A'(a, H*). 


—E’(H’) C.E.(He*)* 











@ Ay A,’ 

10 49.7 22 15.8 1.2 
20 84.4 27 15.8 1.4 
30 117.1 31 15.8 1.5 
40 148.8 33 15.8 1.7 








*C. E. in Tables Il and I1I—Coulomb energy. 


ponential function. The interesting generaliza- 
tion of taking different radii of action for the 
forces between like and between unlike particles 
can also be treated by the same methods, but 
would involve a large amount of numerical 
computation. 


SEcTION V. THE MIXED EXCHANGE OPERATOR 


The experimental cross section for the elastic 
scattering of slow neutrons in hydrogen has 
about the value* 30 X 10-*4 cm’, which exceeds by 
a factor of eight the theoretical cross section® 
computed from the neutron-proton potential de- 
termined in the preceding section. One of the 
writers has attempted to explain this descrepency 
by supposing that the potential depends on the 
state of the system and on the velocity of the 
colliding particles in such a way that the inter- 
action potential for scattering is smaller than 
for binding and decreases with increasing relative 
velocity.’ This hypothesis, if correct, should ap- 
ply equally well to the interaction of two protons. 
Thus we should expect the anomalous proton- 
proton interaction potential required to explain 
the anomalous scattering of protons in hydrogen" 
to be smaller than K(r). In fact K(r) is very much 
too small to explain White's results.” Thus the 


8’ Dunning, Pegram, Fink and Mitchell, Phys. Rev. 47, 
970 (1935); Bjerge and Westcott, Proc. Roy. Soc. A150, 
790 (1935). 

®FII, Eq. (5). 

FIT. 

1M. G. White, Phys. Rev. 47, 573 (1935). 

2 R. D. Present, to be submitted to the Physical Review. 
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hypothesis advanced in FII must be abandoned. 

Wigner" has suggested that the very large ex- 
perimental neutron-proton scattering cross sec- 
tion indicates a dependence of the interaction 
potential on spin orientation. Eq. (2) represents 
the mathematical formulation of Wigner’s sug- 
gestion.'* In the scattering problem based on the 
exchange operator of Eq. (2), there occur four 
distinct wave functions: 


¥1(12) = o(ri2)5(mi, 3) 5(me, 3), (14) 
¥2(12) = ga(ri2)5(mi, 3)6(me, — 3) 
+ yn(ri2)6(m1, —3)6(me, 3), (15) 


and two additional functions obtained by revers- 
ing the directions of the spin coordinates in y 
and ye. It is readily found that ¢ and ga+¢> are 
solutions of the equation 


(Vi2?+AW+J (riz) Pio™ } W(ri2) = 9, (16) 
while gu — ¢» satisfies the equation 
{Vi2k? + W+(1 —22)S (riz) Pio} W(ri2) =0. (17) 


An appropriate set of boundary conditions is 
provided by the following statement: the func- 
tions ¢g and ¢, have the form of an incident plane 
wave plus a scattered wave: the incident wave 
does not occur in g» which contains only a scat- 
tered wave. The scattering intensity obtained 
from ye is simply the sum of the intensities as- 
sociated with yg, and yg, because there is no inter- 
ference between scattered waves associated with 
orthogonal spin wave functions. From this fact 
can be deduced that the statistical weights of the 
scattering intensities obtained from the potential 
fields J(r) and (1—2g)(Jr) are in the ratio three 
to one. The scattering cross section for slow neu- 
trons can be made as large as desired by choosing 
g to bring the lowest eigenvalue of Eq. (17) 
sufficiently near the top of the potential hole, or, 
equally well, by taking g so large that Eq. (17) 
just fails to have a discrete state. For the effective 


radius determined in the preceding section these . 


conditions are met with g~0.2. 
The effect of the mixed operator on the three- 
and four-body problems may be studied most 


18 Unpublished; mentioned by Professor Bethe at the 
Washington Conference on Nuclear Physics. 
4 Van Vleck, Phys. Rev. 48, 367 (1935), note 20. 
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simply when g is small by means of a first order 
perturbation calculation. When g vanishes, the 
alpha-particle wave function can be written in 
the form 


¥(1, 2; 3, 4) =9(x1, x2; X3, X4) 


« S(m, me2)S(mz3, m4) (18) 
with 
S(m,, me) =1/2"!?| 5(my, 3)6(me, — 3) 

—65(me, 3)6(m1, —3)}. (19) 


The average value of the mixed neutron-proton 
potential is simply 


4> fo fva.2i34 
Spin 


X {(1—g)Pis”+¢Pis¥!}¥(1, 2;3,4)dr, (20) 
which reduces to 
4(1-—g, 2) f- . +f ole Xe} X3, X4) 
X J (riz) Piz” (x1, Xo; X3, x4)dr (21) 


when the summation over the spin coordinates is 
carried out. Thus in this problem for small 
values of g the mixed operator can be replaced by 
an “equivalent’’ pure Majorana operator: 


(1—g)P"™+gP"#—(1-—2/2)P”. (22) 


The same result is obtained from the analogous 
treatment of the three-body problem. 

't is clear from Eq. (13) that the calculations 
in Section IV determine essentially the quantity 
2A,,+A, and not A, directly. Consequently 


2Aret+A,~(2—g)Arst+Ay, (23) 


A,~A,+gA), (24) 


and 


if A, is written for the amplitude of the interac- 
tion between like particles when g differs from 


zero. With 1/a'?~2.2XK10°" cm, A,~26, 
A,,~74 and g~0.2, we get 
A,~4lmce?. (25) 


Present” has computed the anomalous scatter- 
ing of protons in hydrogen to be expected from 











INTRANUCLEAR FORCES 


the potential defined by Eq. (25) and finds rough 
qualitative agreement with White’s measure- 
ments."' To summarize: Wigner’s suggestion that 
the neutron-proton interaction potential depend 
on spin orientation does not appreciably change 
the binding energy problem, but does enable us to 
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understand the remarkable experimental results 
on the scattering of slow neutrons and fast pro- 
tons in hydrogen. It should be mentioned that 
the theoretical cross section for elastic scattering 
of fast neutrons in hydrogen is somewhat larger 
than the rather uncertain experimental value. 


APPENDIX 


COMPARISON OF IMPROVED VARIATIONAL CALCULATIONS FOR TWo- AND THREE-Bopy PrRoBLEMS 


Further justification for the use of the method of the 
“equivalent” two-body problem for calculating energies 
for three-body problems is obtained when the results of 
attempts to improve the variational calculations in the 
two problems are compared. The “equivalent” two-body 
problem is determined in both the Wigner and Majorana 
cases by identifying certain constants in the expressions 
for the energy obtained when simple Gaussian-type wave 
functions are used. This identification is possible because 
the energy expressions in this “zeroth” approximation 
have the same form. 

We have used three different methods for improving the 
variational calculations beyond the zeroth stage. In the 
first method the original Gaussian functions are multiplied 
by polynomials in the various distances. For the two-body 
problem we then have for the wave function (taking only 
quadratic terms) 


¢(r) =e7 2) {ag+ayrr?}. (26) 


A similar wave function for the three-body problem is given 
by the expression 


¢(1, 2, 3) = 7 (/2)(r12? +713") H/2) naa” | Cyt Cy(r122 +113?) 
+ Cyre;*}. 


The calculations were made by using the Majorana 
operator for the three-body problem (H*) including the 
neutron-neutron attraction term. The parameter p was 
set equal to unity in order to simplify the calculation of 
the matrix elements. With a given 8, B was determined to 
give E’ = —16 for the exact eigenvalue of the “equivalent” 
problem. The corresponding values of A,x, A», @ for the 
three-body problem are taken from Eq. (13). In the three- 
body problem, the inclusion of the quadratic terms gives an 
improvement in the energy of only a small fraction of an 
mc? unit over the value obtained with the simple Gaussian 
function. The fact that the more flexible three-parameter 
function is no better than the two-parameter Gaussian 
function simply means that the quadratic terms serve only 
to correct vy and uw. When these parameters are both varied, 
the quadratic terms lose their importance. Correspondingly, 
for the ‘“‘equivalent” two-body problem, the quadratic 
term gives no appreciable improvement. This result is of 
value because it brings out the closeness of the correlation 
between the three-particle problem and the associated 
“equivalent” problem. 

In general, computations with Majorana interaction 


(27) 


terms in the Hamiltonian are more involved than those with 
the ordinary Wigner type of potential. Moreover, as we 
show in the next paragraph, there is a Wigner problem 
closely related to every Majorana problem. For the pur- 
pose of showing the close correlation of three-body prob- 
lems to their “equivalent” two-body problems it is just as 
instructive and very much simpler to use Wigner interac- 
tion potentials. 

Let us consider the following closely related three-body 
problems. Call problem (a) a problem in which we have 
only neutron-proton potential terms of the Wigner type in 
the Hamiltonian: 


(28) 


— Aemans® — Aan’, 


In problem (b) take only neutron-proton potential terms 
of the Majorana type: 


—Ae-@n’ PM — Ae-ans’ PM, (29) 


Let problem (c) be the same as (b) except for the addition 
of a neutron-neutron potential term to the Hamiltonian: 


—(A—C)e-arie* PM — (A —C)e-@71s" PM —2Ce-arn®, (30) 


The close relation between problems (a) and (b) was shown 
in FI; (a) gives slightly lower energies than (b). On the 
other hand, problems (b) and (c) are both associated with 
the same ‘‘equivalent”’ two-body problem provided the 
parameter ? is set equal to unity since then for both B = 2A. 
Protilem (c) gives slightly lower energies than problem (b) 
since in passing from the ‘“‘equivalent” problem for (c) 
to that for (b) the parameter p is fixed at unity. We thus 
see that problem (a) and problem (c) are very closely 
related. For the rest of the discussion we shall use the 
simple Wigner-type Hamiltonian of problem (a) since the 
results so obtained are quite general and because to do so 
greatly simplifies the computations. 

We sought to improve the variational calculations by 
using wave functions which are the sum of functions of the 
Gaussian type but with different v's and y's. If we take 


¥=Yothyo’" (31) 
then appropriate wave functions for the two-body problem 
are obtained by substituting for yo 


vo = Ne~rl2), (32) 


v=Bo, 


and for yo’ a similar expression having v’=8e’. The 
matrix elements then are 
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(a/1/o0')2= {400'/(a+0')?}34=y, 
(33) 
(¢@/H/o’)2=y{38e0'/(o+0') —B((a+0’)/(a+o0'+2))*?}. 


In our wave function for the three-body problem we use for 
yo the function given by Eq. (6), with y=ae(1+mn)/(1+2n) 
n=y/v, and for yo’ a similar expression. The matrix 
elements are simplified if we specialize the function by 
taking n’=y'/v’ =n, thus reducing the number of param- 
eters by one. This simplification makes it possible to 
put the matrix elements in the same form as those of the 
two-body problem. In fact we have 


(a/1/o"); = {4a0'/(a+0')?} 9? =¥?, 
(34) 
2+3n+n? oo’ -\* 
(o/H/a'),=3"| 3a £2 _p4 (ste) | 
1+2n oa+0’ ao+o'+2 { 


The value of » which minimizes the zeroth approximation 
(¢/H/o); to the three-body energy is very nearly }, 
independent of A, a and ga. As in the case of the polynomial 
method, B and 8 are determined so that E’ for the two- 
body problem has the value —16. A and a@ are such as to 
establish an “‘equivalence”’ between the two problems in 


the zeroth order and are given by the relations 
B= (28/15)a, B=2A. 


Although A is the amplitude of the neutron-proton 
interaction, it is not directly connected with the A,, 
considered in Section IV since here we have a Wigner 
Hamiltonian and no interaction between like particles. 
In this problem A plays the same role as A,x+}A, in the 
actual problem. Taking the same value of ¢ and  through- 
out (¢, best value for the zeroth approximation, m = 4) and 
varying o’ we obtain about 19 and 17 percent improvement 
in the two- and three-body problems, respectively. If 
a, o’ and » are all varied the improvement is 51 percent 
in the three-body problem and 86 percent in the two-body 
problem. A part of this difference must result from the 
assumption n’=m which restricts considerably the flexi- 
bility of the three-body function. Moreover the three-body 
wave function contains only even powers of the distances 
12, 713, 23 and hence lacks an essential part of the true 
solution. For these two reasons it is thought that the 
difference in the improvements found in the two problems 
is without significance. 

The Hassé method™ was the third method used to 
improve the wave functions and the one which gave the 
greatest improvement in the case of the three-body prob- 
lem. The Hassé method has the great advantage that it 
automatically introduces into the wave function essential 
terms which may be lacking in the zeroth approximation. 
Suppose that in either problem H is the Hamiltonian and 
Yo is the Gaussian wave function; we take as our new 
wave function 


1° H. R. Hassé, Proc. Camb. Pil. Soc. 26, 542 (1930). 


(35) 
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¥ =(1+AH)Yo, (36) 


where \ is a parameter to be varied. In the two-body 
problem y is a two-parameter function while in the three- 
body problem it is a three-parameter function. 

The matrix elements of H, H?, H* and unity are neces- 
sary to the calculation. For the “equivalent” problem 
they are (with B=2A, 8 =(28/15)a) 


(0/1/0) =1, 

(0/H/0) =2.8a0 —2A(o/(o+-1))*2, 

(0/H?/0) = 13.07 (ac)? —0.4a@A (¢ / (a +1))*2(280 + 56) 
+4A?(o/(o+2))*?, 

(0/H*/0) =85.37(aec)*—5.23a°A(a/(o+1))*?? 

xX (150?+520+52) 

+11.2a@A*(o/(o+2))*°*(36+124+4/c) 
—8A%(o/(o+3))*?2, 


(37) 


The matrix elements of unity and H are given correctly by 
Eq. (37) for both problems. Those of H? and H® in the 
three-particle problem are (with = }) 


(0/H?2/0) = 10.67 (ac)? —0.4aA (0 /(o +1))*/2(280 +43) 
+2A*[(o/(o+2))9!? + (0?/(o?+20+15/16))*7], 
(0/H*/0) =51.77(ae)*—a®A (a/(o+1))7? 
X (6407+ 184.80 + 143.3) 
+0.2a@A?[(o?/ (2+ 20+15/16))5/2(840+ 258+ 165/c) 
+(a/(o+2))5/*(840+258+60/c) ] 
—2A%[(o/(o+3))*2+3(62/(o2+3e+15/8))*/?]. 


We find about 65 percent improvement in the energy for 
the two-body problem with this method and about 60 
percent improvement for the three-body problem. It is 
interesting to note that the matrix elements have some- 
what the same form in the two problems except for different 
numerical factors and certain correction terms. The 
numerical values of the matrix elements (0/H?/0) and 
(0/H3/0) are quite different in the two problems because of 
these differences. The fact that the corrections obtained are 
a few percent smaller for the three-body problem than for 
the two-body problem is thought to be of no significance 
because one of the parameters which appear in the three- 
body problem was fixed at a probable value and not varied 
further. 

These investigations establish the inequality 


E(H*) < E°(H*) +0.6{ £’(H*) — EH *)} 


and lend strong support to the view that E(H*) does not 
differ very much from E’(H4). 


(38) 


(39) 


This paper was begun at Harvard and com- 
pleted at the summer session of the University of 
Michigan. We wish to express our appreciation 
of the courtesy and hospitality shown by the 
Michigan Physics Department. 
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Simple approximate calculations of the ranges and 
strengths of the interactions of nuclear particles lead to the 
following conclusions: The interactions are of an exchange 
nature with the range approximately 2.8<10~ cm. The 
neutron-neutron and proton-proton interactions are ap- 
proximately equal. The neutron-proton interaction proba- 


bly depends upon spin orientation. For parallel spins the 
order of magnitude of the interaction is 40 me, for anti- 
parallel spins 25 mc*. It is further concluded that a di- 
neutron or di-proton may be dynamically stable, but 
unstable with respect to 8-decay. 





[* a recent publication,' the writer gave 
arguments for the existence of neutron- 
neutron and proton-proton (to be denoted here- 
after as y—v and r—7) exchange interactions 
comparable in range and strength with the »y—7 
interactions. It was found that the empirical 
data on binding energies of heavy nuclei could 
be explained fairly well by taking all these 
possible interactions of the same range and 
strength. In terms of a potential ‘hole’ the 
range was found to be just &/mc?=2.810-" 
cm, and the strength 40 mc*. It is the purpose 
of the present note to give very simple approxi- 
mate calculations of these quantities and to 
indicate how the assumption of equal inter- 
actions needs to be modified. 


A. RANGE OF INTERACTIONS 


From experimental data on nuclear radii we 
can obtain considerable information as to the 
range of the interactions among the nuclear 
particles. For if, as is certainly the case, the 
interactions are of an exchange nature the forces 
holding a-particles and paired-neutron groups 
together will depend on wave function over- 
lapping in much the same manner as valence 
binding in molecules. In such molecules the 
interatomic distances are related simply to the 
radii of the atoms concerned (d;;~7r;+7;). Thus 
in a nucleus we should expect each a-group and 
v’-group to occupy its own volume, so that the 
nuclear volume, V, is of the order A Vo. This is 
in agreement with experiment.? Further, Vo 


1L. A. Young, Phys. Rev. 47, 972 (1935). In this paper 
a factor 4 is omitted in the proposed nuclear potential 


energy. 

3 For this reason we give below numerical calculations 
of range of interactions for only one nucleus. Any other 
nucleus could have been used. 


cannot be far from the volume of the sphere of 
action of the y—z, v—v, and m—7a exchange 
forces. This means that we can use experi- 
mentally determined nuclear sizes to get at the 
radius of action of these forces. 

The lead nucleus has a cross section for the 
scattering of fast neutrons’ ¢=5.7X10-** cm’, 
corresponding to a collision radius of 1.35 10-" 
cm. This collision radius we may expect to be of 
the order R,+r, where R, is the nuclear radius 
measured to the outermost layer of »v*-groups 
and r is the radius of action of the exchange 
forces. There are 60 groups of paired neutrons 
in the Pb nucleus, so R,~(60)'v, whence 


r~2.7X10-" cm. 


Similarly the radius of the same nucleus extra- 
polated from a-decay data is ~7.810-" cm. 
This radius probably measures the distance from 
the center of the nucleus to somewhere within 
the outermost layer of a-particles, i.e., 7.8 10-" 
cm ~R,—ir. There are 40 a-particles in this 
nucleus, so R, ~ (40)!7, whence r~ 2.6 10-" cm. 
The agreement of these two values for r does not 
mean that either is very accurate. Actually by 
such calculations we can only say that r lies 
between 2.4X10-" cm and 3.010~-" cm. The 
fact that nuclear radii from neutron scattering 
are, in general, larger than radii from a-decay 
data makes it seem reasonable to believe that 
extra neutrons in nuclei are in »*-groups which 
form something of the nature of a “‘mono- 
molecular”’ layer over the core of the nucleus 
which is made up of a-particles. 


8J. R. Dunning, G. B. Pegram, G. A. Fink and D. P. 
Mitchell, Phys. Rev. 48, 265 (1935). 
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B. ORDER OF MAGNITUDE OF y—v INTERACTION 


The more accurate knowledge of 
weights (particularly of the lightest elements) 
given by transmutation data‘ makes it possible 
to gain some idea of the neutron-neutron inter- 
action in a very direct way. Bethe gives 


Ox; oa Org = 1 .0040. 


atomic 


Cis a Cie = 1.0032 > 


These differences correspond to the addition of 
a single neutron to a nucleus made up of a- 
particles. Subtracting these from the neutron 
mass gives the binding energy between the 
unpaired neutron and the rest of the nucleus. 
From the above this binding energy is ~8 mc?.° 
If we now take mass differences for nuclei which 
differ by one v?-group we can hope to get at the 
v—v interaction. Unfortunately, Bethe did not 
have such data to include in his study of masses. 
Determinations have been made with the mass 
spectrograph, however, and yield for Nex». — Nevo 
and Cl3;—Cls,; the value 1.9980. This corre- 
sponds to a mass decrease of 36 mc? from the 
mass of two free neutrons. Of this ~16 mc? can 
be accounted for as due to the binding of the 
two neutrons to the core of the nucleus and 
20 mc? remains to be attributed to »y—v binding 
within the v?-group itself. This does not mean 
that the free di-neutron is stable with a binding 
energy ~ 20 mc’. It is probably just stable. Such 
a great difference between the binding energy of 
a free di-neutron and a v?-group in a nucleus 
can be understood if one remembers that the 
motion of the two particles in a free di-neutron 
would be of a very constrained type (the center 
of gravity must remain at rest or move uni- 
formly) with a reduced mass, 3M. Assuming 
that the radius of interaction for the »—v forces 
is just €/mc? the energy quantization condition 
for a free di-neutron is 


tan (0.312[v—w }!) = —[v/w—1} 


where v is the depth of the potential hole and w 
the binding energy, both in units mc’. The effect 
of binding the di-neutron to a nucleus would be 


4H. Bethe, Phys. Rev. 47, 633 (1935). 

5 This is the average of two values which differ by 
almost 20 percent; so from these data alone one cannot 
say that the binding of a single neutron is independent 
of the number of a-particles in the nucleus. Mass differ- 
ences for heavier nuclei, however, give considerable weight 
to the belief that this is approximately true. 


LLOYD A. 





YOUNG 


to increase the effective reduced mass from }\/ 
to practically 7, so that the quantization condi- 
tion would become 


tan (0.312[2(v—w) }!) = —[v/w—1 }}. 


If, in this equation we give w the value 20, we 
can solve for v. We find v~45. Thus, if the 
v—v interaction be represented by a potential 
hole of radius e&/mc*? the depth is of the order 
45 mc. Since v*-groups have zero spin this is 
the y—v interaction for antiparallel spins. Since 
such paired spins are the usual thing it is hard 
to find a way to determine how two neutrons 
(or protons) interact when their spins are parallel. 
The above calculation is quite rough, but in the 
next section we give a comparison of the »—vp 
and v—7 interactions which should be fairly 
reliable. 


C. COMPARISON OF y—v AND y—f 
INTERACTIONS 


The nuclei Lis, Bio, and Nis consist of one 
(vr)-group outside of a-groups. The next nucleus 
of this sequence, Fis does not exist, but Os, 
having a v?-group, is stable. This may be con- 
sidered as evidence that the »y—7 interaction is 
slightly smaller than the v—7 interaction, al- 
though v?-groups finally become more stable due 
to Coulomb interactions. The Coulomb inter- 
action of the proton of the (yz)-group in Fis 
may be estimated to be ~6 mc? and this means 
that the effective internal binding of the (vz)- 
group is ~7 mc? greater than that of the v*-group. 
This will be reflected in terms of the interactions, 
so that we may say that the »—» potential hole 
is about 7 mc? shallower than the »y—7 potential. 


D. THE FILLED POTENTIAL HOLE 


The periodic system of elements ends with the 
natural radioactive families which are unstable 
with respect to a@ and §-decay. If we visualize 
such a nucleus in terms of neutrons and protons 
moving in a spherical potential hole, then this 
instability is due to the fact that the energy 
levels of this hole are filled up to a point so near 
the top that such transformations are energeti- 
cally possible. Up to the present, evidence points 
to the fact that the effective potential for 
neutrons and protons differs only by the Coulomb 
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interactions which affect the latter. Thus, if the 
effective depth of the neutron potential is #, the 
effective depth for a proton is @—J where J is 
the average Coulomb energy of a proton within 
the nucleus. This is given by J~(1/2)Zé/R, 
where R is the radius of the nucleus or of the 
potential hole. Now according to the Thomas- 
Fermi theory, the number, WN, of particles 
occupying states with momenta between O and 
P is given by 


» —f ” 8rV -) 
7 h 0 P p= 3 h 





Since we are going to consider the potential hole 
for the heaviest nucleus, we write 


Z max = (84 V/3)(P,/h)', 


(A —Z) max = (84 V/3)(P,/h)$ 
with 
P,=(2M(#—Zeé/2R)}!} and P,=(2M®)!. 
The uranium nucleus is the heaviest, so if we 
substitute 92 for Zmax and 146 for (A —Z) max 
we have two equations to solve for @ and V 
(or R). The solution is 


R~8.1X10-" cm; &~ 140 mce?. 


It is indeed satisfactory that such a simple 
approximate calculation gives such reasonable 
results. It should be emphasized that # is not 
the magnitude of the »—v, »y—z, or r—7z inter- 
actions, although these may be deduced from ®. 
The value of R which we have just obtained is 
somewhat too small. It should be nearer to 
9X 10-8 cm. Utilizing this value and only one of 
the above equations we find @~110 mc?. From 
this value of J the heavy particle interaction 
strength may be deduced in the following 
manner. In the uranium nucleus 184 particles 
are in a-groups and 54 neutrons are in v*-groups. 
Each particle in an a-group feels the full force 
of interaction with 3 others while each particle 
in a v?-group feels only one particle. Therefore 
the average nuclear particle feels about 2.6 other 
particles, and the average strength of these 
interactions is (110/2.6) mc? or ~40 mc’. 
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E. NEUTRON-PROTON SCATTERING 


The very large neutron-proton scattering cross 
section observed’ in slow collisions cannot be 
understood on the basis of the neutron-proton 
interaction which explains the binding of the 
deuteron. This probably means that this inter- 
action depends on the relative orientation of the 
particles’ spins. An interaction of strength ~40 
mc? and range e/mc? leads to the observed mass 
defect of the deuteron, but gives a scattering 
cross section ~10~** cm’, which depends very 
little on the energy of collision. In the deuteron 
the neutron and proton have their spins parallel, 
while in collisions both the parallel and anti- 
parallel situations occur. To explain the cross 
section we may therefore take a different inter- 
action for the latter case. An interaction of 
strength ~25 mc*® and range é&/mc® will be 
capable of explaining the collision data. 

We are therefore led to the following conclu- 
sions concerning the interactions among nuclear 
particles: 

1. They are of an exchange nature with a range ~2.8 
X 10-8 cm. 

2. The y—vand x—7z interactions are approximately equal. 

3. The »—7 interaction probably depends on spin ori- 
entation. 


4. The strengths of the interactions are of the order of 
magnitude given below. 


v—-T spins 40 mc? 

y—-r spins a— 25 mc? (22-28) 
v-—voraw—x spins a—|| 33 mc? (30-36) 
v-vora7—7 spins || ? 


It is possible for a di-neutron or di-proton to be dy- 
namically stable, but unstable with respect to £- 
decay. 


wm 


It should be emphasized that the above results 
are provisional and will certainly have to be 
modified somewhat. In the first place, the inter- 
actions are expressed in terms of an effective 
“potential hole.”” This will have to be replaced 
by some continuous interaction function or 
functions. Also, the strengths given above de- 
pend on the range, which may need to be 
modified slightly. 
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Prompt publication of brief reports of important 
discoveries in physics may be secured by addressing 
them to this department. Closing dates for this 
department are, for the first issue of the month, the 


On the Predissociation in the Sulphur Bands 


The phenomenon of predissociation was observed for the 
first time by V. Henri! in the sulphur bands and studied 
later by several authors. The principal feature consists in 
the fact that the emission band system abruptly stops at a 
definite value v’, namely, v’=8, and also that the last 
observed band (8,0) is much shorter than the others. The 
upper limit for the dissociation energy which may be 
derived is D=4.41 ev; G. Herzberg? has shown, by con- 
sidering the character of predissociation, that the right 
value of the heat of dissociation, is in this case, probably 
not much lower than this limit. 

It is, however, impossible to determine the right value 
of D from spectroscopic data if the abrupt break of rota- 
tional structure has not been observed at least in two 
successive vibrational levels; the overlapping in the ultra- 
violet bands makes such a study difficult. We have tried 
to solve the problem by analyzing the bands (6,28), 
(6,29), (6,30), (7,29), (7,30), (7,31), (8,30), (8,31) and 
(8,32). The value of B30” being lower than that of 
Bo” and much nearer to that of Bs’ the overlapping in the 
red region of the band system will not be so important. 
The dispersion used by us up to the present was not suf- 
ficient for a rotational analysis; however, it gave good 
evidence for the abrupt termination of the rotational 
structure at 35 cm~ from the head in the v’ =8 progression 
and at 92 cm in the v’=7 progression. The bands of the 
v’=6 progression were almost too broad and the high 
rotational lines were blended by other bands. 

In the (8,0) band at 2829A we have observed the ter- 
mination of rotational structure at 150 cm@! from the 
head; according to Naudé and Christy’ this corresponds 
to the rotational line K=37. Considering the widths of 
the bands (8; 30,31,32) and (7; 29,30,31), and still using 
the data of Naudé and Christy we may conclude that 
there is a predissociation in the v’=7 level at about K =60. 

If we apply the reasoning established by G. Herzberg? 
in the case of P2 it is easy to show that we have to deal 
with the J, case of predissociation, i.e., the perturbing 
curve producing predissociation has a minimum; the heat 
of dissociation which may be derived in this case with 
great precision is equal to 4.41+0.02 ev. 

Recent investigations by Badger,‘ Van Dijk and Lameris® 
seem to show that Naudé and Christy’s values of B’ and 
B” are too large. If this is correct, the values derived by 
us for K should be modified, but the results concerning the 
shape of the perturbing curve and the energy of dissoci- 
ation should be the same. Our observations give an approx- 
imative determination of a=(B,’’—Bpo”’)/v; we find 


a=0.002. This value may be noticeably in error but the 
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twentieth of the preceding month, for the second 
issue, the fifth of the month. The Board of Editors 
does not hold itself responsible for the opinions 
expressed by the correspondents. 


value a=0.0007 given by Naudé and Christy is certainly 
much too low. 
A high dispersion analysis is in progress. 
Boris ROSEN 
Maurice DésiRANt 
JuLes DucHESNE 
Department of Astrophysics, 
University of Liége, Belgium, 
September, 1935. 


1V. Henri, Structure des Molécules (Paris, 1925), p. 93. 

2G. Herzberg, Ann. d. Physik 15, 677 (1932); note on p. 701. 

3S. M. Naudé and A. Christy, Phys. Rev. 37, 490 (1931). 

4R. M. Badger, Phys. Rev. 46, 1025 (1934); see also A. Christy, ibid. 
47, 251 (1935). 

5 E. W. Van Dijk and A. J. Lameris, Physica 2, 785 (1935). 





After Effect of Aluminum Bombarded by Electrons 


Using a high voltage Lenard tube, J. J. Livingood and 
A. H. Snell! have searched for the radioactivity induced by 
electron bombardment, and having been unable to observe 
it they have attributed the fact to the following four cases: 
(a) The range of the radiation is less than 4 cm of air, or 
the periods are either (b) a small fraction of a second or 
else (c) many hours, or (d) the effect does not exist at the 
potential used. (In their experiment the voltage applied 
was 845 kilovolts.) The present investigation has been 
done with 250-300 kilovolts tube voltage and 1 milli- 
ampere tube current. The high speed electrons produced 
were bombarded upon an aluminum foil placed in high 
vacuum. After the bombardment, the tube potential was 
removed and the after effect of aluminum was examined 
with a Geiger counter. A certain kind of radiation was then 
observed which decayed with half-value periods of about 7 
seconds, 40 seconds and 10 minutes. This radiation was 
determined later to be an electron radiation, and in the 
case of 40-second period, the maximum energy was 2.5 
kilo-electron volts. These results are not inconsistent with 
those of Livingood and Snell whose investigation has not 
been extended to the range of present investigation as 
mentioned above. The similar effect has also been observed 
in the case of many elements including Ni, Cu and Ag. 
It is under investigation whether these are nuclear phenom- 
ena or extra-nuclear ones. This investigation has been 
started by the suggestions of Dr. Y. Nishina, and the 
author’s thanks are due to him and to Professor S. Nishi- 
kawa for their kind advice. 

MASAMICHI TANAKA 

Research Laboratory, 

Tokyo Electric Co., 
Kawasaki, Japan. 
October 20, 1935. 
1 Livingood and Snell, Phys. Rev. 48, 485 (1935). 
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Inhomogeneous Fields for Mass Spectrography 


The use of a nonuniform electric field perpendicular to 
a homogeneous magnetic field provides the possibility of 
tremendously increasing the dispersion in the mass spectral 
analysis of ions. A detailed discussion of the case where the 
fields are uniform has already been given' wherein it is 
demonstrated that the trajectories followed will be cycloids. 
When, in particular, an electric field E,° superposed on a 
magnetic field H, are chosen, the parametric equations for 
the ion paths are 


x =ao+ (eE,°/mw*)+(A/mw) sin wt, 
¥ =Bo— (eE,°/mw)t+(A/mw) cos wt, 


where w=el7,/mc and ao, Bo, A are constants of integration. 
Depending on the initial conditions, the cycloids will be 
curate, prolate or common; and will be the latter when 
éo=cE,°/H,., Yo=20=0. Then, the cusps will lie among a 
line parallel to the y axis; Ay./y.=Am/m; further, the 
addition of the electric field does not disturb the direc- 
tional refocusing characteristic of ion paths in a uniform 
magnetic field. 

Let E, now be a function of y. In a qualitative way one 
can readily see what the resulting effects would be. Suppose 
that E, increases with increasing negative values of y. 
Since the paths followed by the heavier ions will be in a 
region of greater E, than those of the lighter ones, the 
former will spread out more along the y axis, and the dis- 
persion will increase progressively with m in a manner 
depending on the form of E,(y). Rigorous solutions of the 
equations of motion are difficult to obtain in this case; a 
method of successive approximations is applicable. 

Let the zeroth approximation be expressed by Eqs. (1), 
and choose the linearly varying field 


E,=E,+ey, E,=E,+ex. (2) 


Substitute in (2) the results of (1) obtaining E, and £, as 
functions of the time only, as a first approximation. The 
equations can then be rigorously solved. One finds 


xX=a,+e(E,°+8o)/muw?+ (1/mw) (Ki sin wt —K2 cos wt) 
+(e/mw)(E, °+eao)t, (3) 
y=B, +e(E,°+ cao) /mua*+ (1/mw)(K, cos wt — K» sin wt) 
— (1/mw) {e(E,°+¢Bo)t — (e*eE.°/2mw)t* 
+(eeA /mw*) sin wt}, 


(1) 


where a, 6:1, K, and Ke are, once more, constants of in- 
tegration. It is a simple matter to solve for the extremal 
values of x and y as given in (3). One finds that the greater 
the number of periods taken for these modified cycloids, 
the more nearly correct becomes the relation ¥n~m?; a 
large increase in the dispersion has, therefore, been made 
possible. 

It does not appear desirable at this point to go into 
details as far as the trajectories given by (3) are concerned 
until higher approximations have been investigated and 
other functions giving rise to greater dispersion con- 
sidered. It is to be pointed out, however, that in order to 
obtain the full benefit of the large dispersion brought about 
by the nonuniform field, sharpness in the definition of the 
traces will probably be sacrificed. Hence it would be of 
great interest to study the refocusing characteristics of such 
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a field combination. From an experimental point of view, 
the necessity of increasing the pole-gap separation as well 
as the precise determination of the nonhomogeneous 
electric field may offer some difficulties.? 
S. H. BAUER 
Gates Chemical Laboratory, 
California Institute of Technology, 
October 31, 1935. 


'L. Page, Phys. Rev. 33, 553 (1929); 42, 101 (1932). 
2 F. W. Aston, Mass Spectra and Isotopes (Longmans, Green and Co. 
1933), p. 77; A. E. Shaw, Phys. Rev. 44, 1006 (1933). 





A New Experimental Method for Study of the Upper 
Atmosphere 


Calculations and laboratory tests indicate that by pro- 
jecting a modulated searchlight beam and measuring the 
light scattered from sections of its high altitude path at an 
observing station some kilometers distant, by using a large 
mirror with a photo-cell and synchronized amplifier, 
information can be obtained regarding the molecular 
density and certain other characteristics of the atmosphere 
in the almost unexplored region from 30 to 70 km, and 
possibly for even greater heights. A brief description of the 
method and its possibilities is given in the December, 1935 
issue of the Journal of Terrestrial Magnetism. This notice is 
published to reach the attention of workers in numerous 
other fields who may be interested. 

M. A. TUVE 
E. A. JOHNSON 
O. R. WuLF 
Department of Terrestrial Magnetism 
Carnegie Institution of Washington (M.A.T., E.A.]J.), 
Bureau of Chemistry and Soils 
U.S. Department of Agriculture (O.R.W.), 
Washington, D. C., 
November 6, 1935. 





Doubly Excited States of Helium—A Correction 


In a recent article, Wilson! noticed a discrepancy 
between his value for (2p)?*P obtained by Hartree’s 
method, and that obtained by one of us by the variational 
method.? This discrepancy is unfortunately due to an 
algebraic error in the latter calculation. The correct value 
obtained by minimizing the integral for (2p)?'P is 
(2p)? *P = —1.3976, and the normalized wave function for 
3P is 
(8/9) Vio(1 ax®!*r; exp (—ar:) Vi: (2)8*'*r2 exp (—Brz) 

— Yio(2)a*/*r2 exp (—arz) ¥i1(1)8*/*r, exp (—ar,)I, 


where a=0.84, and 8=0.83. The values for (2p)? 'D and 
(2p)? 1S obtained with the variational parameters so deter- 
mined are (2p)? 1D = —1.326, (2p)? \S= —1.217. Thus the 
value (2p)? ?P = —1.3976 from the variational method is 
quite close to Wilson's value — 1.4018. 
Ta-You Wu 
S. T. Ma 
National University of Peking, 
October 24, 1935. 


1Wm. Wilson, Phys. Rev. 48, 536 (1935). 
2Ta-You Wu, Phys. Rev. 46, 239 (1934). 
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The Value of the Electronic Charge 


Kellstrém's' recent determination of the viscosity of air 
seems to furnish the final piece of evidence required to 
establish the essential correctness of the so-called grating 
value of the electronic charge e. As is well known, there has 
been for a number of years a sharp discrepancy between the 
value of e given by oil-drop work and that deduced from the 
absolute wave-lengths of x-rays, as determined by means of 
a grating. In recalculating the oil-drop data, in 1929,? I 
adopted Millikan’s assumed value of the viscosity of air 
(nos = 1822.7 X 10-7), and also his assumed probable error 
(0.05 percent), since I did not feel competent to express an 
expert opinion in that field. This value had been obtained 
by the method of rotating cylinders, a method that has 
generally been considered more reliable than the method of 
capillary tubes. I thus obtained e=(4.768+0.005) x 107'° 
e.s.u., as the final result of Millikan’s oil-drop work. 

Shiba* was perhaps the first to challenge the correctness 
of Millikan’s assumption regarding the viscosity of air. He 
presented a table of values of » which indicated clearly that 
the capillary tube method gives in general a definitely 
higher result than the rotating cylinder method. Shiba him- 
self adopted 1831.2 x 1077 as the best average value. Now, 
however, Kellstriém, using the method of rotating cylinders, 
obtains (1834.8+3.0) 10-7, an even higher value. The 
use of this last value with Millikan’s data leads to e 
= (4.816+0.013) K 107" e.s.u. 

Meanwhile, as Bearden‘ has recently shown, all meas- 
urements on the absolute wave-lengths of x-rays aie 
remarkably consistent. These results can be used to deter- 
mine a value of e only if one assumes a geometrically 
perfect crystal (in practice, calcite), and certain observed 
constants (density, etc.) for this crystal. With such an 
assumption, and with the latest values of the constants, one 
obtains, in agreement with Bearden,‘ e = (4.8036 +0.0005) 
xX 10-'° e.s.u. Hence there is no longer any outstanding 
discrepancy between the values of e determined by these 
two distinctly different methods. 

This fact, however, does not settle the problem of the 
values of the three interrelated constants e, e/m, and h. 
As a result of a number of recent investigations on e/m, 
it seems more than probable that its true value lies between 
1.757 and 1.758 (107 abs. e.m.u.). Let us adopt 1.7575 
107 as the best average. To get a value of the Planck 
constant h, one may now use Bohr’s formula for the Ryd- 
berg constant, a formula that is still believed by theoretical 
physicists to be correct to a high degree of accuracy. With 
e=4.8036 and e/m =1.7575, one obtains h =6.6286 x 10727 
erg - sec. 

About a year ago’ I presented a diagram showing all 
important experimental results involving e, e/m, and h. 
This diagram indicated only too clearly the impossibility 
of assigning any set of values to these three constants 
that would even reasonably satisfy all of the experimental 
results. This situation still remains essentially unchanged. 
The set of values that I gave then was e=4.768, e/m 
= 1.7574, h=6.547. A few of the experimental results are 
equally well satisfied by either set of values, but in general 
the two sets predict very different experimental results. 
Thus if the new set (e =4.8036, etc.) is correct, the experi- 
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mental results for the limit of the continuous x-ray spec- 
trum, as a function of voltage, are all in error by amounts 
up to 20 times the assigned probable error. These various 
results, however, are quite consistent, and the method is 
usually considered especially reliable. Similarly the average 
value of h/e, from determinations of ionization and 
resonance potentials, is in error by some five times its 
apparent probable error. The determinations of h/e from 
the photoelectric effect and from the value of the radiation 
constant C2, are equally in error. 

It appears to me that the most desired experiment in this 
field, at the present time, is a really reliable determination 
of h/e. Since the two sets of values of e and h that have been 
discussed give values of h/e differing by 0.58 percent, it 
should be possible to distinguish between them. 

It may be noted, in closing, that the new set of values of 
e, h, and e/m, leads to 1/a = 137.06. If one retains the same 
value of e/m (1.7575), but shifts e from 4.8036 to 4.810, 
1/a becomes exactly 137, and h becomes 6.6433. 

RayMonpD T. BirGE 

University of California, 

November 15, 1935. 

1G. Kellstrém, Nature 136, 682 (1935). 

2R. T. Birge, Rev. Mod. Phys. 1, 1 (1929). 

3K. Shiba, Inst. Phys. and Chem. Research, Tokyo, Sci. Papers 19, 
97 (1932). 


4J. A. Bearden, Phys. Rev. 47, 883 (1935). 
5 R. T. Birge, Science 79, 438(A) (1934). 





Note on Majorana’s Exchange Energy 


The exchange energy introduced by 
Yis,J(rij)Pi; = V, where 7, 7 refer to protons and neutrons, 
respectlvely, and P;; is the exchange operator which inter- 
changes the space coordinates x;', x2", x3' of the proton 7 
with the space coordinates £)/, 7, &3/ of the neutron j, and 
ri; is the distance between 7 and j. The interaction energy 
operator V is Hermitian and it commutes with the three 
components of the total momentum operator 


(h/i)(2;0/dx,' + 2,0/dE,) =G,. 


Majorana_ is 


It is satisfactory in these respects. There is, however, an 
undesirable feature of this interaction energy which is due 
to the difference of the mass u of the neutron and the mass 
m of the proton. It shows itself in an improper behavior of 
the center of mass and a lack of invariance to Galilean 
transformations. The coordinates of the center of mass 





X,=(2mx,'+ Zyué,’)/(Zim+ Zu) have the following equa- 
tions of motion 
dx, G; a(m—p) y . ° 
= — +———_ 2 SI (rj) (Ef —x,*) Pj. 
ei OM en Oe ee 


Here M is the total mass. If m= this is the usual relation 
between the velocity of the center of mass and the mo- 
mentum. Because of the extra term in m—y there is an 
additional tremblatory component in the motion of X, 
somewhat analogous to Schridinger’s ‘‘Zitterbewegung.” 

The lack of invariance to Galilean transformations can 
be seen by considering the deuteron. The solutions cor- 
responding to an energy E and a momentum G are obtained 
by the substitution y= ¢(y) exp {(¢/2h)2(x,+£,)G,} with 
y=&—x which gives 
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The operator on the left side of this equation is Hermitian 
and E—(G?2/8)(1/m+1/) has real eigenwerte. It should 
be noted that the substitution y= x(y) exp {(7/A)z, 
X (mx,+put,)G,/(m+u)} leads toa non-Hermitian operator 
for E—G?/2(m+u). There is thus no freedom in the choice 
of the substitution and the energy values are determined 
by Eq. (1). If instead of using the Majorana operator one 
used ordinary potential forces, the term in 4.—m would be 
absent and (G?/8)(1/m+1/u) would be changed into 
G*/2(m+ yu). These two effects do not compensate in 
general because the change in E due to the term in u«—m 
depends on ¢ and hence on J. It is thus impossible to 
interpret the eigenwerte of E as the sum of G?/2(m+y) 
and an internal energy independent of G. 

These anomalies are clearly due to the fact that every 
exchange of the positions of a neutron and a proton 
changes the position of the center of mass. The anomalies 
can be avoided by modifying the operator. A simple change 
leading to the desired result is to use instead of P;; an 
operator P;;’ defined by 





=-——. . . 
soo, xt+ (t?—x"*), ), 


M 
utm 
This operator P’ leaves unchanged the center of mass of 
the neutron-proton pair. It is Hermitian, commutes with 
the total momentum, and with the center of mass. The 
wave equation modified by changing P into P’ leads to 
results invariant under Galilean transformation. In ordi- 
nary calculations of mass defects and similar problems it 
is practically indistinguishable from P on account of the 
smallness of (u—m)/(u+m). 

It is not probable however that such a complicated inter- 
action energy can have the significance of a fundamental 
law of nature. It would be more satisfactory to have a 
general theory in which the Majorana force is the first 
term of an expansion in (u—m)/(u+m). 

G. BREIT 
E. WIGNER 
The Institute for Advanced Study, 
Princeton University, 
November 14, 1935. 





Proton-Proton Forces in Anomalous Scattering and in 
Nuclear Binding 


Recent experiments by White! on the scattering of 
600-750 kv protons in hydrogen show large deviations from 
the theoretical cross sections calculated by quantum 
mechanics and indicate a departure from the Coulomb law 
at a distance of about 5X10-" cm. Using the Feenberg- 
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Knipp nuclear model, discussed elsewhere in this issue, we 
have calculated the scattering cross sections on the assump- 
tion that the interaction between two protons includes, in 
addition to the usual Coulomb repulsion, an attractive 
interaction —A,xe~*®” with a=2.210-" cm and Agr 
adjusted to fit the experimental scattering at 45°. The value 
of A,;x so determined (45 mc?*), is in very satisfactory 
agreement with the value obtained by Feenberg and Knipp 
(41 mc*), when the binding problems are treated with a 
mixed Majorana and Heisenberg potential for neutron- 
proton interaction. These calculations have been repeated 
with a rectangular potential hole under varied conditions. 

When the calculations are extended to other angles of 
scattering we encounter the gravest difficulties. It is 
readily seen that when the phase shift Ko, which repre- 
sents the effect of the modified field on the s scattering, is 
adjusted to fit the observed scattering at 45°, it is impossi- 
ble to get agreement with the scattering at other angles. 
This result is entirely independent of the type of potential 
assumed for interaction at small distances. Both the 
theoretical and the observed scattering show a minimum 
in the neighborhood of 30° in contrast to the Mott- 
Coulomb scattering which decreases monotonically with 
increasing angle up to 45°. However, the theoretical 
scattering at 30°, when a fit is made at 45°, is over six 
times as large as the experimental. White’s results are 
undoubtedly subject to large statistical fluctuations, as 
may be seen from irregularities in the data and from the 
fact that at 30° the observations correspond actually to 
about four proton collisions. Nevertheless the disagreement 
seems too considerable to dismiss cavalierly in this way. 
One alternative that remained open was to investigate the 
possible influence of p scattering on the angular distribu- 
tion. If this were to be of importance an immediate conse- 
quence would be the breakdown of the Coulomb law at the 
unbelievably large distance of 210~-" cm. Examination 
reveals that the modification in the scattering produced by 
a distortion of the p wave is zero at 45° and much too small 
at 30° to reduce the scattering to White’s value. Here one 
adjusts the phase shift Ko for s scattering to give a fit at 45°, 
and then the phase shift K, for p scattering is varied to 
minimize the scattering at 30°. If a smaller value of Ko is 
used so as to obtain a theoretical scattering at 45° about 
0.8 of the experimental, it is still impossible to fit the 
scattering at 30° by a suitable choice of K,. It is to be 
emphasized that these conclusions are completely inde- 
pendent of the nuclear model. They are true even if the 
potential is not centrally symmetric or depends in some 
way on the relative momentum or spin orientation. More- 
over, they must remain true even if the wave mechanics 
breaks down for small distances, since only the asymptotic 
form of the wave function is involved (in a region where 
wave mechanics is known to be accurate) and any future 
theory must reproduce the same expression. Distortion of 
the higher order harmonics would be entirely negligible, 
for the de Broglie wave-length of the fastest protons 
(7X10-" cm) is so large compared with the range below 
5X 10-" cm where departures from the Coulomb law occur 
that the wave function must have there an appreciably 
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constant amplitude. Thus it appears impossible in any 
way to account for the experimental results. Further 
experimental work is urgently needed to confirm or pre- 
clude the possibility that these discrepancies may be due 
to statistical fluctuations or systematic error. 

A discussion of the calculations with the nuclear model 
will shortly appear in a more extended form in this journal. 
This work was originally suggested to me by Dr. Eugene 
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Feenberg, to whom I am grateful for much friendly and 
valuable discussion. I also wish to thank Dr. M. G. White 
for communicating to me unpublished data. 
R. D. PRESENT 
Purdue University, 
Lafayette, Indiana, 
November 13, 1935. 


1M. G. White, Phys. Rev. 47, 573 (1935). 
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